Pathogen detection using next generation sequencing

ABSTRACT

Embodiments are directed to systems and methods for pathogen detection using next-generation sequencing (NGS) analysis of a sample. Embodiments may apply alignment algorithms (e.g., SNAP and/or RAPSearch alignment algorithms) to align individual sequence reads from a sample in a next-generation sequencing (NGS) dataset against reference genome entries in a classified reference genome database. Embodiments of the present invention may include classifying, filtering, and displaying results to a clinician that can then quickly and easily obtain the results of the sequencing to identify a pathogen or other genetic material in a sample that is being tested. A negative sample and a corresponding database can be used to remove contaminants from a list of candidate pathogens. Thus, embodiments are directed to a system that is configured to filter the results of a sequencing alignment and classify a sample quickly.

CROSS-REFERENCES TO RELATED APPLICATIONS

The present application is a continuation of International Patent Application No. PCT/US2016/052912 filed Sep. 21, 2016, which claims the benefit of and priority to U.S. Provisional Application No. 62/221,574, filed Sep. 21, 2015, the entire contents of which are herein incorporated by reference in their entirety for all purposes.

GOVERNMENT LICENSE RIGHTS

This invention was made with government support under grant no. HL105704 awarded by the National Institutes of Health. The government has certain rights in the invention.

The Sequence Listing written in file SequenceListing_081906-220810US-1078352.txt created on Mar. 7, 2018, 765 bytes, machine format IBM-PC, MS-Windows operating system, in accordance with 37 C.F.R. §§ 1.821-1.825, is hereby incorporated by reference in its entirety for all purposes.

BACKGROUND

Infectious diseases affect the lives and health of millions of patients annually. Failure to obtain a laboratory-confirmed diagnosis for many acute infectious diseases directly contributes to poor patient outcomes and a high cost burden to the health care system. Key areas of unmet clinical need include neurological infections (encephalitis and meningitis), pulmonary infections (e.g., pneumonia), blood infections, and sepsis.

Traditional diagnostic methods, including culture, antigen detection, and nucleic acid amplification, are limited in scope in cases for which there is little clue regarding the identity of the causative agent. DNA sequencing is the process of determining the precise order of nucleotides within a DNA molecule. It includes any method or technology that is used to determine the order of the four bases—adenine, guanine, cytosine, and thymine—in a strand of DNA. DNA sequencing can be helpful in identifying potential causes of disease in patients. For example, alignment processes may be used to identify matching portions of a sample sequence with a reference database of classified reference sequences.

However, DNA sequencing of a sample and the alignment of such sequences for identification of a potential source of disease, typically includes long processing times due to the large amount of information to compare and process in order to identify a matching sample sequence and reference sequence. Additionally, because of the vast amount of data within sequencing data, sequence alignment results can return large numbers of false positives where portions of sequence reads may appear to match to portions of a reference genome that is not in fact present in the biological sample. As such, many alignment results from a DNA sequence alignment process may not be accurate and raw results of such alignments are not useful as-is in a clinical environment because an expert and/or clinician must analyze the results and manually interpret the returned alignment results to interpret the results of the sequencing reads.

Thus, challenges in the field include developing accurate pipelines (or processes) that can quickly analyze millions of reads that include millions of data points that come out of a DNA sequence system, as well as interpret the data so that it is clinically useful to laboratory scientists and/or a physician. Accordingly, there is a need for systems that are capable of quickly and efficiently identifying and interpreting next generation sequencing data for detection of potential causes of disease and/or any other potential applications of DNA sequence alignment information.

Embodiments of the present invention solve these and other problems individually and collectively.

BRIEF SUMMARY

Embodiments are directed to systems and methods for pathogen detection using next-generation sequencing (NGS) analysis of a sample. Embodiments may apply alignment algorithms (e.g., SNAP and/or RAPSearch alignment algorithms) to align individual sequence reads from a sample in a next-generation sequencing (NGS) dataset against reference genome entries in a classified reference genome database. Various embodiments can filter, classify, and display results to a clinician to identify a pathogen or other genetic material in a sample that is being tested. Embodiments can provide various systems that are configured to filter the results of a sequencing alignment and classify a sample quickly and accurately.

As an example, two alignment techniques (one being faster than the other) can be used together to speed up alignment, without sacrificing accuracy. An initial alignment technique can identify which reference genomes in a database match to which sequence reads. For a matching reference genome, an optimally-aligning sequence read can be identified. For the optimally-aligning sequence read, a different alignment technique can be applied, and it can be determined whether any of the new alignment scores to other reference genomes exceed the optimal alignment score for the matching reference genome. If an new alignment score exceeds the optimal alignment score for the optimally-aligning sequence read, the matching reference genome can be removed from a set of matching reference genomes. The set of matching reference genomes can then be output.

As another example, sequence reads can be assigned to a particular classification level, so as to provide accuracy identified of a particular pathogen. Sequence reads can be identified that match to two or more matching reference genomes of the classified reference genomes with at least the minimum alignment threshold. For such a sequence read, a taxonomy identifier can be assigned from classification information to each of the two or more of the classified reference genomes. The taxonomy identifier can include at least two levels of classification. The assigned taxonomy identifier of each of the two or more classified reference genomes can be compared at each of the at least two levels of classification, with levels that do not match being removed. The lowest level shared between the two or more reference genomes can be assigned to the sequence read. Updated alignment results can be provided and include a number of corresponding sequence reads for each of the plurality of taxonomy identifiers.

As another example, background contaminants can be identified and removed from a set of potential (candidate) pathogens that are clinically-relevant. A negative control sample can be used to identify sequence reads from potential contaminating organisms. A ratio can be taken of a first amount of sequence reads from a test sample that align to a matching genome and a second amount of sequence reads from the negative control sample that align to the matching genome. The ratio can be compared to a threshold to identify a set of one or more matching reference genomes have the ratio exceed the threshold. An output can identify a set of one or more matching reference genomes as potential pathogens in the test sample.

Other embodiments are directed to systems and computer readable media associated with methods described herein.

A better understanding of the nature and advantages of the present invention may be gained with reference to the following detailed description and the accompanying drawings. Reference to the remaining portions of the specification, including the drawings and claims, will realize other features and advantages of the present invention. Further features and advantages of the present invention, as well as the structure and operation of various embodiments of the present invention, are described in detail below with respect to the accompanying drawings. In the drawings, like reference numbers can indicate identical or functionally similar elements.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a block diagram of a system 100 for sample analysis, identification, and classification according to embodiments of the present invention.

FIG. 2 shows an exemplary method of analyzing and matching sequence reads of a sample to one or more classified reference genome databases according to embodiments of the present invention.

FIG. 3 shows an exemplary method of analyzing and matching sequence reads of a sample to one or more classified reference genome databases and identifying biological material corresponding to sequence reads of the sample according to embodiments of the present invention.

FIG. 4 shows an exemplary method of applying a taxonomic classification algorithm to the results of an alignment process to one or more classified reference genome databases for sequence reads of a sample according to embodiments of the present invention.

FIG. 5 shows an exemplary method of filtering the results of an alignment process to one or more classified reference genome databases for sequence reads of a sample according to embodiments of the present invention.

FIG. 6 shows an exemplary table including a summary of the results of an exemplary sequencing process for three different sample sequence reads according to embodiments of the present invention.

FIG. 7 shows an exemplary result table including matching reference genomes to three different sample sequence reads as a result of an exemplary sequencing process without filtering or taxonomic classification processes being applied to the results according to embodiments of the present invention.

FIG. 8 shows an exemplary result table including matching reference genomes to three different sample sequence reads as a result of an exemplary sequencing process with a filtering process being applied to the results but no taxonomic classification process being applied to the results according to embodiments of the present invention.

FIG. 9 shows an exemplary result table including matching reference genomes to three different sample sequence reads as a result of an exemplary sequencing process with a taxonomic classification process being applied to the results but no filtering process being applied to the results according to embodiments of the present invention.

FIG. 10 shows an exemplary result table including matching reference genomes to three different sample sequence reads as a result of an exemplary sequencing process with both a filtering process and a taxonomic classification process being applied to the results according to embodiments of the present invention.

FIG. 11 is a flowchart of a method for identifying pathogens in a test sample of biological material according to embodiments of the present invention.

FIG. 12 shows table of sequencing statistics according to embodiments of the present invention.

FIG. 13 shows a table of the number of reads from DNA libraries and RNA libraries aligning to viral sequences according to embodiments of the present invention.

FIG. 14 shows an exemplary computer system.

An Appendix includes Supplementary Tables 1-4.

DETAILED DESCRIPTION

Embodiments can provide processes for rapid analysis of next generation sequencing (NGS) data for pathogen detection. For example, embodiments may be used in a broad, comprehensive pathogen diagnostic for infectious diseases by analyzing sequencing results against many reference genomes, as a metagenomics analysis. The use of unbiased metagenomic next-generation sequencing (mNGS) can provide for detection of all potential pathogens in a single assay. An advantage is the ability to detect all viruses, bacteria, fungi, and parasites in a single, standardized universal test directly from diverse clinical sample types such as cerebrospinal fluid (CSF), bronchoalveolar lavage (BAL), and plasma, thereby maximizing the potential impact on patients with acute, life-threatening infections by early and more accurate diagnosis.

Embodiments may apply alignment algorithms (e.g., SNAP and/or RAPSearch alignment algorithms) to align individual sequence reads from a sample in a next-generation sequencing (NGS) dataset against reference genome entries in a classified reference genome database. For example, as a result of sequencing of a sample, the system may obtain a list of results of reference sequences that align with the sequencing reads of a sample. Individual sequences in the reference genome database (also referred to as a GenBank) are referred to as reference genome sequences and may be identified by genome identifiers (GI). The genome identifiers include reference identifiers used to identify reference sequence genomes stored in the GenBank database. A best GI match can be assigned to each read, and the taxonomy assigned to each GI according to the classified reference genome in the GenBank can be assigned to each read.

Such unbiased analysis (e.g., using the large number of reference genomes in the GenBank) can make accurate detection of specific pathogens difficult as many human samples include flora or background colonization organisms. For instance, a nasal swab has a lot of bacteria in it because that's going to be bacteria that sort of colonize your respiratory tract. Thus, when the system makes a detection, it is difficult to know whether the system has detected a pathogen or a colonizer. Further, because the analysis does not bias or target any individual pathogen, the system could be detecting, for instance, a bacterial contamination of enzyme preps that are used in making the sequencing libraries. So effectively the system may identify laboratory reagent contamination instead of matching to a particular infectious disease.

Accordingly, while the sequencing analysis and matching is sensitive, the analysis may return false-positive matches (e.g., a human read will align to a viral GI due to database misannotations). And, reference genome databases may include many GI entries that are poorly annotated (e.g., “uncultured eukaryote” instead of a particular identification) or, even worse, incorrectly annotated (e.g., http://www.ncbi.nlm.nih.gov/nuccore/KC506764.1, annotated as GBV-B when this is actually GBV-C virus). These databases typically allow anyone to add reference genomes and there is no standard for annotating samples.

Further, many of the alignments are not accurate and raw results alone are not useful as-is clinically because an expert and/or clinician must analyze the results and manually interpret the returned alignment results to interpret the results of the sequencing reads. Typically, such interpretation required an expert genomicist or bioinformaticist or someone who is well-experienced with infectious diseases or laboratory medicine and who also understands the bioinformatics. Accordingly, such raw results are not helpful or usable in a clinical laboratory where such processes can provide great value to clinicians caring for patients.

To address this issue (non-specific identification), embodiments can filter and classify sequencing reads to obtain more accurate results and avoid false-positives, maximizing the specificity of detection. Embodiments can filter, classify, and display results to a clinician that can then quickly and easily obtain the results of the sequencing to identify a pathogen or other material that is being tested. Filtering can be helpful because the genome databases may be rife with misannotations and false positive matches.

Further, typical sequencing methods compare sequences to reference databases and then identifying hits based on what the sequence is aligned to. However, the actual interpretation of those hits is more subtle because an expert actually has to analyze the data to determine if some hits are not accurate and/or do not provide enough information to be classified as a particular pathogen. Embodiments can improve the classification of reads by applying a rapid taxonomic classification algorithm to alignment results to provide more accurate and clinically useful results. Additionally, embodiments can annotate the data within the genome reference database so that the system knows what reference genomes are a pathogen, a colonizer, what is considered likely contamination, etc. Embodiments can also provide a user-friendly visualization interface that can be used by clinicians to quickly and easily identify pathogens and other results.

Embodiments can provide a number of advantages including limiting the manual annotation and interpretation of sequencing analysis results to allow quick, efficient, and useful clinical and public health surveillance. For example, embodiments may be used where a patient is sick and a clinician may take one or more samples from the patient to determine what pathogens are present in the samples. The clinician may desire to know whether the patient has a viral infection, a bacterial infection, a fungal infection, a parasite, etc. Any type of biological sample that obtains DNA material could be used to identify potential causes of the disease. For instance, blood, cerebrospinal fluid, respiratory secretion, tissue, stool, etc., may be used to obtain a sequence read of the samples. Additionally, embodiments may be used in blood bank testing, food and water quality testing, environmental testing, animal testing, animal health, or any other area that may be assisted by quickly and efficiently determining potential sequence matches within a sample. However, a sample will likely include a mixture of multiple organisms. Some of these organisms may be from a person, viral, bacterial, and many samples are mostly mixtures of different organisms.

The mNGS assay has been clinically validated for diagnosis of encephalitis and meningitis in cerebrospinal fluid. The assay can incorporate: (1) the analytic wet bench process including handling of patient samples, nucleic acid extraction, mNGS library preparation, and sequence generation on an Illumina HiSeq instrument, and (2) bioinformatics analysis of sequence data and clinical interpretation by trained microbiologists/pathologists. Results of diagnostic mNGS testing from cerebrospinal fluid are reportable in the medical chart and can be used for clinical management. The clinical implementation of the mNGS assay has a direct, positive impact on patient outcomes by increasing the number and proportion of patients with an accurate, clinically actionable infectious disease diagnosis that allows for timely management and treatment. The assay is also able to detect rare, unexpected, or slow growing or uncultivable microorganisms, for which diagnosis is often delayed or missed. The mNGS assay may have particular utility in identifying culture-negative pathogens due to prior treatment with antibiotics. The mNGS assay may also be useful as a “rule-out” test for infectious diseases, which may impact management by increasing clinical confidence in working up and treating non-infectious causes of encephalitis/meningitis such as autoimmune disease with steroids and/or immunosuppressive medications

The mNGS data can also yield additional information besides whether or not a given microorganism or microorganism type is present or absence. The number and proportion of reads, typically expressed as a “reads per million (RPM)” metric normalized to a negative “no-template” control sample run in parallel, can provide some degree of quantitative or at least semi-quantitative information. In some cases, the pathogen genome coverage may be sufficient to facilitate (1) precise genotyping or strain identification, (2) analysis of single-nucleotide polymorphisms or mutations, (3) the generation of predicted antibiotic/antiviral resistance profiles.

I. Example Systems

Once a clinical sample has been obtained (e.g., from a patient or from an environmental samples, such as a water sample), the sample may be sequenced, and the sequencing results can be analyzed. Various systems and processes can be used. For example, nucleic acid extraction can be performed. A cDNA/DNA library preparation may include adding adapters, although other preparation for other types of sequencing may be performed, e.g., for nanopore sequencing or other single molecule sequencing. The library of templates can be fed into a sequencing device to provide sequencing information, e.g., base calls or raw signals that are used to determine base calls, thereby obtaining sequence reads. The analysis of the sequence reads can include host subtraction; adapter, quality, and low-complexity trimming; and alignment against reference databases such as NCBI (National Center for Biotechnology Information) GenBank.

FIG. 1 shows a block diagram of a system 100 for sample analysis, identification, and classification according to embodiments of the present invention. As shown, system 100 includes a sample analysis system 110 and a sequence identification computer 120. Sample analysis system 110 may be communicatively connected to the sequence identification computer 120 through a communications network 140, which may be any suitable communications network for communicating data. As another example, data may be communicated via a removable storage device, such as a USB drive.

The sample analysis system 110 may be configured to receive a sample 150 (e.g., after library preparation) from a clinician or other operator, sequence the sample to obtain a plurality of sequence reads of genetic material for the sample, and submit the plurality of sequence reads to sequence identification computer 120. A sequencing device 115 can correspond to any sequencing device, such as those produced by Illumina, Pacific Biosciences, or Oxford Nanopore. A processor 111 (e.g., a CPU) can control aspects of sequencing device 115, such as one or more cameras for taking images of the nucleic acids or electrical components, both receiving sequencing signals corresponding to nucleotides of the nucleic acids being sequenced.

A memory 112 (e.g., flash memory, hard drive, DRAM, cache, etc.) can store software for controlling processor 111, which can control sequencing device 115. In some embodiments, a sample collection module 113 can control robotic processes for obtaining a sample (e.g., via a syringe connected to a robotic arm), and perform any automated preparation processes. A sample sequencing module 114 can instruct processor 111 to perform the sequencing, using sequencing device 115. Sample analysis system 110 can perform analysis of the raw sequencing signals (e.g., fluorescent or electrical signals) to identify basecalls of sequence reads, or send such raw sequencing signals to sequence identification computer 120.

Sequence identification computer 120 can process, align, and identify genetic material that is present in the sequence reads to identify pathogens and/or other genetic material that is present in the sample. An alignment module 123 in memory 122 can instruct processor 121 to align the sequence reads to a plurality of reference genomes, e.g., as stored in reference genome database 130. One or more alignment techniques (e.g., a local aligner, such as SNAP) can be used to obtain alignment results, e.g., initial alignments results, as well as subsequent alignment results. A classification module 124 can classify the alignment results to include an accurate taxonomic classification level. A filtering module can filter the alignment results to remove false positives. Sequence identification computer 120 can return results for identifying pathogens and/or other genetic material present in the sample.

Reference genome database 130 can contains a plurality of reference genomic sequences that have been identified and classified as being associated with a particular biological organism. For example, reference sequences of different viruses, bacteria, fungi, human, animal, and/or any other reference DNA sequences of any other biological material may be stored in the classified reference genome database. Sequence identification computer 120 may apply one or more alignment techniques to align the received sequence reads from the sample to the classified reference genome sequences stored within reference genome database 130.

Although FIG. 1 shows the modules and describes the functionality of the modules as being completed by the two systems, in some embodiments, all or a portion of the functionality may be split differently and/or could be completed by a single device. For example, in some embodiments, sample analysis system 110 of FIG. 1 may be coupled to a reference genome database and/or the reference genome database may be stored as part of sample analysis system 110, and sample analysis system 110 may complete all of the functionality described herein.

System 100 can be used to identify and detect pathogens in clinical samples, e.g., for the purposes of surveillance, such as epidemiologic surveillance. For example, it may be desired to look at cases of sequences generated from patients with acute diarrheal illness or fever to identify and detect pathogens. Surveillance can also be performed to detect novel pathogens, thereby allowing pathogen discovery, e.g., when no match is found in reference genome database 130. Another application is use as a diagnostic tool, e.g., by identifying a known pathogen for which treatment is known.

Metagenomic sequencing can include alignment to multiple reference genomes, as may be stored in reference genome database 130. One example of such a database is GenBank, which has more than 20 gigabases in size. With such a large size, the alignment of aligning millions of sequence reads can provide many hits (matching alignments) to different genomes, which is why classification and/or filtering can be useful. There may actually be a pathogen in that sample that is identified by a hit, but there can be false hits, as may occur due to problems in the databases. Databases themselves are not well curated, and there are incorrect or errors in how they are constructed. For example, there may be sequences annotated as hepatitis C virus but may actually be a human sequence for instance and vice versa. The database itself can be cleaned up, as well, e.g., when misannotations are identified, as may be done using BLAST.

Additionally, a sample can include a large amount of background nucleic acids. Techniques can be used to filter out such background sequences, as may be done with a no template control (NTC) sample. NTC database 135 (also called a contaminant database) can store background sequences that have been identified in recent samples, e.g., samples within the last month. Levels (e.g., number of reads) that a background sequence has been quantified can be stored, as is described below.

For the taxonomic classification, individual sequences can be classified according to where they fit on the taxonomic level. For instance, if you have sequences that are aligned to certain regions (e.g., viral or bacterial genome) those individual reads may not be specific for that species. That is, it may not be specific to that viral genome of a viral species. In such a case, the read can be classified using a least common ancestor algorithm to the next higher taxonomic level. Once the sequence reads are classified, the system can then point out what species are specific to the sample, thereby allowing a good prediction that an actual viral species is in the sample.

In some embodiments, the classification can use the results of an initial alignment technique, e.g., SNAP, thereby allowing the classification to start while alignment is being performed for other sequence reads. The efficient alignment can relatively concurrent classification can allow a taxonomic classification in one to two hours. Accordingly, if the alignment results provide 100 hits, the classification can narrow to about 10 to 15. However, for clinical purposes, it is desirable to have only a few hits, e.g., two. The 10-15 hits may be real, but may not be clinically significant. For example, the hits may include microorganisms that are part of the laboratory background contaminants or part of the skin background as part of drawing blood. Techniques to filter out such contaminants are described in more detail below, as well as alignment and classification techniques.

II. Data Analysis

Once the system has obtained the raw sequence reads from the sequencing module, the raw sequence reads are analyzed to obtain an identification of matching reference sequences from one or more classified reference genome databases. Embodiments can provide data with respect to what sequences from different microbial pathogens are detected.

A. Example Pipeline

FIG. 2 shows an exemplary analysis process 200 for processing the sequence reads, identifying matching reference genome sequences, and determining matching viral and bacterial matches for the sample sequence reads. Process 200 may be performed with any suitable computer system that is connected to a sequencing machine, e.g., via a network or a removable storage device, such as a USB drive.

At step 210, the raw sequence reads arc received from the sequencing module. The sequencing process may return a large number (e.g., 1, 2, 3, 5, 10, 20, 50, or 100 million reads, or more) of DNA sequence reads from a sample. As examples, the raw sequence reads may be received in fasta, fastq, sam, or bam files. The sequence reads can be obtained in various ways, e.g., from a sequencing of nucleic acids or probe-based methods, such as hybridization arrays. The use of random sequencing that essentially sequences all nucleic acids can be preferable so that a significant number of pathogens (e.g., as many as reference genomes are available) can be detected.

At step 220, the raw sequence reads are preprocessed to trim the reads to remove low quality sequences, remove low complexity sequences (e.g., as would not be very informative), remove adaptors that can be retained at the ends, etc. Thus, the reads may be cleaned-up to ensure that the remaining reads are of high quality. If adaptors are used, the adaptors (adapters) may be ligated to the end of the nucleic acids and used as primers for sequencing, and thus may not relate to the actual nucleic acids. The low complexity reads (e.g., with many repeats) can be difficult to align, and thus removed. The quality of a read can be determined based on the quality of the basecalls, which may be determined based on the sequencing signals.

At step 230, an alignment module (e.g., alignment module 123) performs a sequence alignment analysis for the sample sequence reads to the reference genomes. In some embodiments, the alignment analysis may be performed in reference to a host genome database first and any sample sequence reads that align to the host may be removed from the sequence. For example, in embodiments that are identifying pathogens associated with a human sample, the sequence reads may be aligned with a reference human genome database and any matches may be removed (since no pathogens should be present in the human reference genome database). Comparison can be made to multiple reference human genomes.

Thus, embodiments can use computational subtraction to speed up the identification processes to remove any sequence reads that align with the host (and thus are not helpful for identifying biological material from other entities (i.e., pathogens)). In embodiments where a single large reference database is used, any results that are identified as matching with a human reference genome could be removed. However, by applying the full 100% of the sequence reads and not subtracting them before aligning with the full reference database, the process may be slower than if those sequence reads were removed before comparing to the large database of reference genomes.

The host database (e.g., human genome database) may include fewer genomes than the entire reference genome database. Accordingly, sequencing can be faster than comparing to all possible reference genomes. Thus, for relatively high quality sterile samples, applying such a host filtering method may remove most of the background reads quickly. For example, for a good sample, the removal of the host matches could take the analysis from potentially 100% of the reads down to maybe less than 10% of the reads. Thus, 90% of the sample sequence reads can be removed by alignment to a host (e.g., human) database.

Additionally, in some embodiments, a host database could include similarly related genomes that are not specifically from the host. For example, for the human example provided above, the system could also include primates that have similar genomic references to humans. This provides a more comprehensive host analysis and provides even faster and more effective subtraction of host matches.

Next, depending on the mode of the identification process, either a comprehensive mode or a fast mode can be performed to identify one or more matching reference genomes associated with the sample sequence reads. The fast analysis may only align the sample reads to a reference genome database of bacterial and viral reference genomes. Thus, the fast mode analysis may not identify all potential genomic matches in the sample sequence reads but will focus on the analysis on potential bacterial and viral matches due to the focus of the process on identifying pathogens within the samples. The comprehensive mode may align the sample reads to an entire nucleotide classified reference genome database that includes reference genomes from bacterial 240A, fungal 240B, parasitic 240C, viral 240D, and other 240E reference genomes.

At step 240 of the comprehensive mode, the alignment module performs a sequence alignment analysis for the sample sequence reads to the reference genome database, e.g., to reference genomes 240A-240E, Steps 241 and 242 of the fast mode may only perform alignment to the bacterial database 240A and viral database 240E, respectively. Any number of different sequencing algorithms may be used in embodiments of the present invention. For example, the “Scalable Nucleotide Alignment Program (SNAP)” algorithm includes a nucleotide aligner that takes raw sequence data and aligns it to nucleotide reference databases. SNAP is extremely fast and by using fast sequencing algorithms, the analysis processes of the present invention may return results extremely fast. While other analysis methods (i.e., “pipelines”) may analyze a sample sequence data within days, weeks, or months, embodiments of the present invention can analyze the sample sequence data in minutes to hours. Fast sequencing analysis is critical in a number of applications including, for example, infectious diseases analysis that can be paired with a next generation sequencing essay that can diagnose infectious diseases and get results to physicians regarding patient care within 8-12 hours, or as soon as possible

For example, in some embodiments, the system may align the sequence reads from the sample to all classified reference genomes in GenBank using a SNAP alignment algorithm. However, the GenBank is growing very fast (e.g., doubling every year) so it can be important to limit the number of reads being analyzed to improve the speed of the alignment analysis with the reference database (e.g., GenBank). For example, in some embodiments, a specific subset of the GenBank may be used to improve analysis speeds further. For example, if the clinician is not concerned about potential plant matches, such reference genomes may not be applied to the alignment analysis and/or a separate subset or different database may be used to avoid potential plant genetic matches. Accordingly, some embodiments may align to a tailored search for potential pathogens, bacteria, etc. and avoid aligning to reference sequences that are not a part of the tailored search.

At step 250, in the comprehensive mode, de novo contig assembly may be performed to obtain contigs associated with the results of the assignment from the classified reference genome database. Example assembly software include ABySS and Minimo.

At step 260, in the comprehensive mode, the reads and the contigs are used to align the translated nucleotides using another alignment algorithm (e.g., RAPSearch) and compared to a viral protein reference database. To generate translated nucleotides, the sequencing reals are translated in all 6 reading frames and the resulting amino acid sequences are then compared to a protein reference database. The output of RAPsearch is similar to that of BLAST but lists alignments that exceed a predesignated E-value significance threshold.

At step 270, the results of the alignment algorithms are obtained, and may include a summary table 270A that provides the matching viral, bacterial, fungal, parasitic, and other reference genomes that have aligned to the sample sequence reads. Thus, in response to the alignment, the system can provide a table that shows alignments to bacterial and viral reference genomes. Summary table 270A may list all the hits to the different species in the sample and the different genera in the different families of organisms found in the sample, e.g., according to taxonomic classification 270C. For example, the system may divide the results into viruses, bacteria, non-cordate eukaryotes (basically fungi and parasites that do not have a backbone), human, as well as an “other” category. Cordates are all higher level eukaryotes, eukaryotic organisms that have a backbone. Eukaryotic organisms that do not have a backbone are more often microorganisms. These categories are used in order for the system to identify fungi and parasites as well as other microorganisms that do not have backbone (e.g., invertebrates like worms). The system may find matches that are not necessarily microorganisms, for instance, the system may be capable of diagnosing a tape worm infection.

Thus, at the end of the alignment process and in an initial table, the system may provide matches for different bacteria and may include the number hits for each type of bacteria, virus, etc. Furthermore, in some embodiments, de novo assembly may be applied to actually recreate the genome and coverage maps 270B of the genome may be provide to indicate how well the reads cover the genome of one the genomes in the list.

Accordingly, the results of the alignment process of FIG. 2 can provide all possible matches for one or more DNA sequence reads from a sample. For example, summary table 270A may include fifty bacteria and/or twenty viruses. Next the question becomes which one of these or are all of these causing the potential infection? Some of these results may be contaminates, some might be colonizers, etc., so the results must be interpreted to identify a potential candidate as the pathogen that is causing an infection. Traditionally, if a system detects three different bacteria, a clinician has to interpret the matches to determine which one of those three bacteria or which combination of bacteria are the most important in terms of causing the disease. Accordingly, it is important to be able to limit as many potential false-positives and incorrect matches as possible to assist clinicians and users of the system to identify the most likely one or more pathogens responsible for an illness.

Accordingly, embodiments of the present invention provide (1) filtering, (2) taxonomic classification, and (3) best match algorithms for identifying and providing the most useful data to a clinician. Additionally, embodiments provide additional features including RNA sequence removal and on-the-fly annotation of database entries to further assist interpretation of alignment results.

B. Flowchart of Example Data Analysis

FIG. 3 shows a streamlined process 300 for identify, matching, and interpreting sample sequence read alignments with classified reference genome databases according to embodiments of the present invention. Aspects of process 300 may be performed in a similar manner as process 200. Process 300 may be performed using system 100.

At step 310, a system may receive a sample from a patient and/or other biological item and/or entity. The sample may be provided by a clinician or other operator of the sample analysis system. The sample may form nucleic acids that have already been prepared into a library for sequencing.

At step 320, the system may obtain sequence reads of nucleic acids from the sample. Any suitable method for obtaining DNA sequence reads may be used, e.g., as described herein.

At step 330, the system may preprocess the sequence reads. The preprocessing can include trimming sequence reads and removing some reads, e.g., low quality or low complexity reads.

At step 340, the system may apply an alignment algorithm to the preprocessed sequence reads and obtain alignment results for the sequences reads. Each read can be aligned to a classified reference genome database, e.g., including millions of reference genome sequences. The alignment may return millions of “hits” or alignments to the classified reference genome database.

Some of the matches may not be perfect and may include sequence reads that match a portion of a reference genome sequence. A quality value may be returned with the alignment results that indicates a measurement and/or magnitude of similarity between the sequence read and the reference genome sequence. Additionally, individual reads may align with more than one reference genome. Accordingly, alignment results may include an identifier of the reference genome that was matched, classification information associated with the annotated information provided when a reference genome was uploaded to the classified reference genome database, and a similarity measurement indicating the quality and/or closeness of the alignment between the sequence read and the reference genome.

At step 350, the system may provide a taxonomic re-classification of the alignment results for reads that have aligned to multiple reference genome sequences with different taxonomic classifications. For example, the system may compare levels of classification between multiple aligned reference genomes for a particular sequence read and may remove one or more classification levels that do not match between the reference genome sequences aligned with a particular sequence read. Additional details and steps associated with the rapid taxonomic classification process is provided in reference to FIG. 4 below.

At step 360, the system may filter the alignment results to remove false positives and mis-annotated or mis-classified reference genome sequences that aligned to the sequence reads of the sample. For example, the system may select a best match sequence read for an identified reference genome and may apply a second alignment technique to the best matching sequence read to identify if the same reference genome aligned with the second independent alignment. If the alignment of the sequence read does not match the same reference genome, the read and the reference genome may be removed from the alignment results as the reference genome is likely a false positive hit. Additional details and steps associated with the filtering steps are provided in reference to FIG. 5 below.

In some embodiments, the filtering may include removal of a reference genome as a viable pathogen when the taxonomic classification is listed in the results of a no template control (NTC) for the current experiment or a database of previous NTC results. Further details are provided in a section below.

At step 370, the system can use the updated alignment results including the filtered and taxonomically re-classified sequence reads to identify one or more best matches for the sequence reads from the sample. For example, the best match may be selected based on the reference genome sequence that aligned with the most sequence reads.

At step 380, the system provides the results of the best matching analysis to the analysis system being operated by the clinician. The results may be provided via a network connection. The results may include the most probable one or more pathogens that may be responsible for a patient's illness.

At step 390, the analysis system provides and/or displays the results to the clinician. The results can be provided in any variety of ways, e.g., visually or by audio. The results can indicate a treatment to be provided, thereby providing a therapeutic intervention.

C. Example Best Match (GI Picker)

An objective of embodiments can be to pick the most informative one or more genome identifiers (GI) out of the GIs matched by all the sample reads of a single taxonomic assignment, a particular species. As part of the analysis, a score can be determined for each of the matched reference sequences, specifically each GI. In some implementations, the scores are determined for only GIs at a given taxonomic level, e.g., species. In other implementations, the total scores are determined for GIs across various (e.g., all relevant) taxonomic levels, with a top score being picked across all the GIs. The total scores can be determined from respective scores of various properties.

In some embodiments, each of the matched reference sequences can be scored as to any or all of length, coverage, identity, and percent identity. Thus, four scores can be determined and used to determine a total score. Examples of how the scores can be determined are provided below.

Length is the length of the reference sequence. Thus, a length score is larger for longer reference sequences. Therefore, a species of microorganism that has a larger genome would have a longer length score. In some embodiments, matches can be ranked by whether or not complete versus partial genomes are available in the database, where use of only a partial genome would affect the length score.

For each reference sequence, a coverage score can be determined. In one embodiment, the coverage score can be determined as a cumulative score representing a sum of the counts of aligned reads at each genomic position (e.g., base position) along the reference sequence. In another embodiment, the coverage score can be determined using a sum of genomic positions that have at least one read aligned to that position, e.g., regardless if there is a match that position. The coverage score may be a percent of the genome that is “covered” by at least one aligned read, and thus the sum of genomic positions that have at least one read aligned to that position can be divided by the total number of positions in the database for that reference genome.

For reference sequence, the identity score can be determined as a cumulative score representing a sum of fractional scores at each genomic position. A fractional score can be calculated as the number of aligned reads with a nucleotide matching the reference sequence divided by the number of aligned reads whether or not matching the reference sequences (i.e., the total number of aligned reads covering the position). Thus, a genomic position will have a larger fractional score when there are fewer mismatches of reads at that position.

Percent identity can be calculated by dividing the identity score by the coverage score.

Once the individual scores are determined, an overall score can be calculated for each reference sequence (taxonomic identifier). In one embodiment, the overall score can be determined by adding ‘length’ and ‘identity’ and multiplying the sum by ‘percent identity’. The reference sequences can be ranked by the overall score, e.g., in descending order. And, a top score(s) can be selected, e.g., by choosing the first reference sequence in the list. In some embodiments, more than one reference sequence can be selected, e.g., all having a total score above a score threshold. In other embodiments, the top N or X % of taxonomic identifiers can be identified.

Accordingly, for each sample, a score can be determined for each relevant GI (e.g., a set of taxonomy identifiers) at one or more taxonomic levels (e.g., for each species, genus, family, etc.). The scores may be determined only at one taxonomic level, or at a given level and all GIs at lower levels. As an example for each GI, the reads aligned to that GI can be grouped. For each GI, the coverage and identity score can be determined in the following manner. For each alignment of a read, the following calculation can be performed for each alignment position: map the read to the genomic position of the GI, the position coverage score is increased, and the position identity score is increased if there is a match. Across all genomic positions for the GI, the total coverage score can be increased by one for each position that has a position coverage score greater than zero. The total identity score can be determined as the sum of the fractional scores, computed as the position identity score divided by the position coverage score. The percent identity score can be determined as the total identity score divided by the total coverage score.

In some embodiments, the following pseudo code can be used:

-   -   for each sample         -   for each species             -   group by GI all reads and alignment data             -   for each GI                 -   score covered and identity                 -   for each alignment                 -   for each alignment position                 -   map to GI position                 -   covs[position]+=1 if match/mismatch                 -   idents[position]+=1 if match                 -   sum across all GI positions                 -   covered+=1 if covs[position]>0                 -   identity+=idents[position]/covs[position]                 -   score percent identity=identity/covered             -   sort descending by (length+identity)*percent identity             -   choose first in sorted list

D. Mapping Reads to Selected GI

In various embodiments, an output (e.g., a file) can contain the picked reference sequence (GI) and a consensus sequence of the reads mapped to it. The term “coverage” can be used in two contexts: (1) one of the sub-metrics that make up the overall score for picking a GI; and (2) an overall metric, inferred from a coverage map file, for how comprehensively a set of mapped reads cover the picked GI.

For picking a GI, a particular species can have multiple genomic identifiers for different reference sequences in a database. The database of reference sequences used can vary or be constructed from some combination of separate databases. On database is the NCBI ‘nt’ database. It is highly redundant, and often mis-annotated. A single species can have thousands of entries, or GIs. A GI can be picked using embodiments of the best match algorithm above.

After picking a GI for each species, scores can be determined for each of the selected GIs (e.g., one for each species). The scores can be determined using all reads assigned to the species level, and can include reads assigned to higher levels. The scores can be determined as described above, and scores exceeding a threshold, e.g., an absolute number, a value at which the top N scores are above (where N is an integer), or a value at which the top X % of scores are above.

Assigned reads at the subspecies, species, genus, or family level can also be mapped to the selected GI for visualization. Mapping of reads at higher taxonomic levels (genus or family) will increase the coverage of the GI (the percentage of nucleotides that have at least one read mapped to it), but can also incur a risk of erroneous mapping from different species or genera in the sample. Real-time visualization of the coverage maps and pairwise identity plots can facilitate expert interpretation for clinical results reporting.

III. Rapid Taxonomic Classification

Some embodiments may analyze clinical or environmental metagenomic data (e.g., metadata) to classify the origin of each of the millions of next-generation sequencing reads (e.g. bacterial, viral, human, etc.) that are returned by the alignment results. However, individuals reads may be informative only at a given taxonomic level (species, genus, family, etc.). For instance, a read in the conserved matrix gene of influenza viruses may be identical (or nearly identical) in influenza A and influenza B, and thus this read cannot be used to distinguish between influenza A and B. This is problematic if one is interested in species or strain level identification, which is critically important in infectious disease diagnosis (examples: Bacillus anthracia/anthrax versus Bacillus cereus; enterovirus versus rhinovirus; Ebola Zaire versus Ebola Reston).

As such, some embodiments may provide a taxonomic classification algorithm that removes some of the hits from consideration through the metagenomic data provided through the classification of the reference genomes. The taxonomic classification algorithm may apply a least common ancestor (LCA) approach where taxonomic identifiers for each of the associated reference genomes are analyzed to identify the lowest level of shared ancestry between two or more reference genomes. For instance, if a sequence has hits to both a human and a virus, and the match is conserved between a human and a virus, then the algorithm would move the result up to the proper taxonomic level, which would in this case be above the kingdom level. In such embodiments, the system would assume that it mapped to the human and would remove the result since the taxonomic classification moved up to the kingdom level.

As another example, if the system returned hits for a single virus but different versions of the virus (e.g., influenza A and influenza B), and the hits are indistinguishable between influenza A and influenza B, then those sequences would be assigned to the influenza genus, and not to the species level (i.e., Influenza A and B). Thus, the system would move the classification for this reference genome up to a higher level (e.g., family enterovirus instead of mapping to two different enterovirus species—influenza A and influenza B).

Such classification allows a clinician to look at the results at different levels and provide potential options for treatment. Using the influenza A or influenza B example above, the result may indicate that the system matched either influenza A or influenza B instead of influenza A. For example, if the system is analyzing an influenza A sample (i.e., a sample taken from a patient infected with influenza A), the sequence results will likely return that 90% of the sample sequencing reads are going to align influenza A and 10% are going to align to influenza B because of the indistinguishable portions of the genetic material between influenza A and B. Reads may be from regions that are highly conserved between these two species of influenza. Thus, a clinician may not know, looking at the raw data, whether this patient is infected with both influenza A and B, is infected with influenza A alone and it just happens to be that B was misaligned to influenza B, or vice-versa. However, by applying the taxonomic classifier and by properly classifying those reads that were initially assigned to influenza B, the system moves up the classification to the genus level, and the remaining reads assigned to influenza A would allow a clinician to make a determination that the patient has an influenza A infection instead of the erroneous call, which would be this a dual infection with influenza A and B.

Thus, after shared reference genomes have been classified at the appropriate level, the remaining reads in the alignment results should be specific at each taxonomic level and thus, the remaining hits that are species-specific are known to be particular to that species. By removing the alignment results that are shared between multiple reference genomes sharing a genus, family, etc., the results will be more accurate and usable by the clinician. The reads that remain after shared reference genomes of the same genus have been reassigned to a genus as opposed to a particular species associated with the aligned reference genome. Identifying species-specific reads can be extremely informative to clinicians where the clinicians know that the read is specific to the particular species, sub-species, strains, etc. It can also be useful to know the higher taxonomic level information including the number of reads at a particular genus, family, etc. Thus, embodiments provide more informative and accurate information for a clinician.

A. Classification Method

FIG. 4 shows an exemplary process 400 of taxonomic classification for find the appropriate taxonomic classification level for reads and aligned reference genomes that may share portions of genomic sequences between various species, genus, family, subspecies, strain, lineage, etc. Thus, embodiments allow for appropriate identification of possible reads aligned to a reference genome to allow the system to provide the best possible identification of potential pathogens causing illness and/or present in a sample from a host organism. Process 400 may be performed by system 100.

At step 410, the system receives the alignment results from applying an alignment algorithm to the sequence reads from the sample. In some embodiments, the alignment results may have been previously filtered using techniques described herein. The rapid taxonomic classification process may be performed concurrently with the alignment process and thus, before the filtering process has been completed. Thus, the system can receive a plurality of sequence reads obtained from a sequencing of DNA molecules from the sample of biological material where the sample includes DNA molecules from a plurality of organisms and aligns the sequence reads using an alignment technique to align the plurality of sequence reads to a plurality of classified reference genomes in a database.

The system may obtain initial alignment results that include, for each of at least a portion of the sequence reads, a matching reference genome to which the sequence read aligns. The initial alignment results may further include classification information for each of the matching reference genomes where the classification information includes a taxonomic identifier including multiple classification levels for each reference genome. For example, a reference genome (e.g., enterovirus D) may include a species (e.g., enterovirus D), germs (e.g., enterovirus), and family (e.g., picornaviridae) classification level within a taxonomic identifier. Many other classification levels may be also be assigned and included in the taxonomic identifier and corresponding classification information with the initial alignment results for each aligned reference genome.

At step 420, the system identifies a set of the sequence reads that aligned to two or more of the classified reference genomes. For example, in FIG. 7, the human rhinovirus and the human enterovirus K1105 both have a single hit. However, as can be seen in results table 900 of FIG. 9, after classifying the initial results, both of these reference genomes have been re-classified. Accordingly, these reference genomes matched with a read that aligned to multiple reference genomes. Step 420 can ensure that sequence reads that align to only one reference genome are not subjected to unnecessary further analysis, as those sequence reads can be classified based on the single matching reference genome.

Steps 430-490 can be performed for each of the sequence reads identified in step 420.

At step 430, the system identifies two or more matching reference genomes of the classified reference genomes to which the sequence read aligns with at least the minimum alignment threshold. For example, using the example provided above in reference to FIG. 7 and assuming that a sequence read matched to both of these reference genomes, the system may select a read that aligned to both of these species-specific reference genomes. The minimum alignment threshold can ensure that a matching reference genome is sufficiently similar to the sequence read at a given location. The minimum alignment threshold may be used for all sequence reads, or can vary for different sequence reads, e.g., based on the specific alignment quality for each of the matching reference genomes or as a result of taking the highest N matching reference genomes (e.g., as determined by an alignment score, such as a number of mismatches).

At step 440, the system selects a set of reference genomes for a sequence read that matches to multiple reference genomes. For example, using the example provided above, the system may select the human rhinovirus and the human enterovirus K1105 reference genomes that may have been associated with the same read. The selected set of matching reference genomes may correspond to those identified in step 430. In other embodiments, other factors may be used to determine the selection of the set of reference genomes to use. Such other factors can include the magnitude of the differences between the “best” and “second-best” alignments, the genomic sequence diversity corresponding to the detected microorganism (e.g. viral genomes are most diverse than bacterial genomes), and the completeness of the available reference databases, including the availability of “near-neighbor reference genomes (i.e. genomes that are closely matched in the reference database).

At step 450, the system assigns a taxonomy identifier from the classification information to each of the two or more of the classified reference genomes. The taxonomy identifier may include at least two levels of classification (e.g., species, genus, family, etc.) and the levels may have a hierarchy such that there is a lower level (e.g., species) and at least one higher level (e.g., genus, family, etc.). For example, using the example of FIG. 7, the system may assign the taxonomic identifier including a species of human rhinovirus species, a genus of enterovirus, and a family of picornaviridae to the first reference genome and a species of human enterovirus K1105_171204, a genus of enterovirus, and a family of picornaviridae to the second reference genome.

At step 460, the system compares the assigned taxonomy identifier of each of the two or more classified reference genomes at each of the at least two levels of classification. For example, using the example of FIG. 7, the system may compare the species, genus, and family of the human rhinovirus sp and the human enterovirus K1105_171204 reference genomes. The system may begin at the lowest level (e.g., species or sub-species level) and compare whether they are the same.

At step 470, the system removes each level of the at least two levels from the assigned taxonomy identifier that do not match between the two or more classified reference genomes. For example, using the using the example of FIG. 7, the system may determine that the sub-species and species between the human rhinovirus sp and the human enterovirus K1105_171204 reference genomes are not same and may remove those levels from the results. Thus, if differences are present between the levels of classification assigned to each of the multiple reference genomes, the lower levels of classification between the multiple reference genomes may be removed until all the levels match.

At step 480, the system may assign to the sequence read the lowest level of the at least two levels of the assigned taxonomy identifiers that is shared between the two or more classified reference genomes. For example, using the example of FIG. 7, the system may determine that the lowest shared level that is shared is the genus level of enterovirus for the human rhinovirus sp and the human enterovirus K1105_171204 reference genomes. Accordingly, the system may assign the lowest level of the taxonomic identifier for the two reference genomes to be the genus level of enterovirus for the two reference genomes. In such a case, the identifier for the genus level would be assigned to the sequence read. If none of the at least two levels of the taxonomy identifier match, an unassigned state may be associated with the sequence read.

As shown in FIG. 9, the species results for both the human rhinovirus sp and the human enterovirus K1105_171204 reference genomes have been removed from the results table, and instead the reads have been combined into a genus match for the enterovirus that has a read of 10,461 sequence read hits. Accordingly, the results indicate that 10,461 sequence reads aligned with reference genomes that were shared between various species within the enterovirus genus, but that were not specific to a particular species. Accordingly, the results are simplified for a clinician.

At step 490, the system determines whether all of the reads matching multiple reference genomes sequences have had taxonomic classification processing applied. If there are additional reads matching multiple reference genomes to be analyzed, the system may repeat the process for the next sequence read in the initial alignment results (and may return to step 430). This process may continue until all of the sequence reads matching to two or more reference genome sequences have been analyzed to ensure the reference genomes have the appropriate level of classification.

At step 495, the system provides an identification of one or more taxonomy identifiers corresponding to one or more candidate pathogens based on numbers of corresponding sequence reads assigned to each of the plurality of taxonomy identifiers. The identification can be the taxonomy identifiers with or without additional information. For example, the number of reads assigned to each of the identifiers may or may not be included. The fact that the identifier(s) are being provided can indicate that the identifier(s) are candidate pathogens, i.e., a sufficiently high likelihood of existing in the sample based on the number of reads assigned. Various criteria can be used to determine whether a pathogen is identified as a candidate pathogen, e.g., using best match algorithm for picking one or more identifiers (as described herein) and criteria for determining whether the best matched taxonomy identifier provides a sufficient match to be identified as a candidate pathogen.

In some embodiments, sequence reads having an assigned taxonomy identifier can be used to identify the one or more taxonomy identifiers corresponding to one or more candidate pathogens. For each of a plurality of taxonomy identifiers, a total score can be determined based on at least coverage of a reference genome corresponding to the taxonomy identifier. The total scores can be ranked, and one or more taxonomy identifiers exceeding a threshold can be identified.

In one embodiment, updated alignment results can be provided to a clinician where the updated alignment results include a number of reads for each of the plurality of taxonomy identifiers. As shown in FIGS. 7, 9, and 10, the effect of taxonomic classification is to properly classify reads to their correct taxonomic level. Two reads (751 and 752 in FIG. 7), one aligning to rhinovirus A 752, and the other aligning to human enterovirus K1105 751, are correctly bumped up to the genus level (genus Enterovirus) in FIG. 9. By doing the taxonomic classification, the interpretation of the data is that there is only one species of enterovirus identified with a large number of sequence read alignments (enterovirus D) in the sample (yellow shading), and not two different species (rhinovirus A and enterovirus D, for example). Table 1000 of FIG. 10 shows a results table 1000 for classification and filtering.

In some embodiments, the classification of a read can begin once alignment results for that read are obtained. Thus, the alignment results for all reads are not required before classification can begin. Accordingly, the two processes may be performed in parallel (e.g., in separate threads) with the output of the alignment thread being used by the classification thread, with the classification thread operating slightly delayed as needing to wait for the alignment results for a first read.

Once the candidate pathogen(s) are identified, additional available clinical and laboratory data can be helpful in determining whether or not the detected organism is pathogenic (i.e. causing disease) in the host organisms (e.g. human patient). The detection of the presence of a potential pathogen in a clinical sample does not necessarily mean that it is causing disease; the potential pathogen could be a colonizer, for instance, or a bystander and have nothing to do with the host organism's illness. If the candidate pathogen is deemed to be pathogenic by clinical and other criteria, the detection can be used to guide clinical interventions, which can include (1) drug therapy (e.g. prescribing or administration of a targeted antimicrobial agent), (2) drug discontinuation (e.g. discontinuing a drug that was administered empirically in the absence of a definitive diagnosis), (3) vaccination, if a vaccine is available and efficacious after infection (e.g. rabies), and (4) medical procedures (e.g. valve replacement in cases of fungal endocarditis, for which antifungal therapy alone is ineffective). The failure to detect a candidate pathogen may also be clinically useful to “rule out” infection as the cause of illness, which can guide clinicians to work up and treat for non-infectious causes (e.g. administering intravenous immunoglobulin and steroids for autoimmune disease).

Accordingly, when the sample of biological material is from a host organism, a clinical intervention can be performed for the host organism based on the identification of the one or more taxonomy identifiers corresponding to the one or more candidate pathogens. The clinical intervention can include the examples above. The intervention can include actually performing/administering any of the above procedures or prescribing them.

B. Minimum Alignment Threshold

The alignment results can provide a distance (e.g., an edit distance) between the sequence read and the matching region of a reference genome. The distance can provide a measure of how many nucleotide differences are there between the actual sequence read and the reference. An edit distance can be the minimum number of editing operations (e.g., insertion, deletion, and substitution) to transform the read into the reference genome, or vice versa. Using the edit distance can help to perform the classification quickly.

Once a distance is obtained, the reference genomes can be ranked by the distance. The taxonomic classification can place the read at different species, genus, or family level based on edit distances to different matching reference genomes. Thus, the taxonomic locations for each of the reference sequences can be determined.

The desired edit distance for determining matched alignments will vary depending on the sequence diversity of the organism. Thus, the minimum alignment threshold can vary depending on which reference genomes are involved. For instance, viruses have more divergent genomes (meaning there is a lot of diversity). For viruses, a relative large threshold can be used, e.g., up to 12 mismatches in the fixed read length. For bacteria and fungi, a smaller distance may be tolerated because bacterial genomes that fit in the same species are highly identical; we want to be able to not only identify the species but we want to be specific.

For instance, if there are two different subtypes of enteroviruses, 70 and 71, there can actually be up to 30% difference by sequence between those two viruses. On the other hand, two bacterial species, such as staph aureus and staph epidermidis, can be 99% identical across the genome. Edit distances can be used as a criteria for classifying these reads, and different taxonomic groups you would require different edit differences. A suitable matching threshold distance can be determined empirically meaning from clinical data or data generated from positive and negative controls where it is known what is in the sample, e.g., what pathogens are in the sample

IV. Filtering Reference Genome Sequence Hits

One embodiment of the present invention is directed to filtering out sequences that appear to be false positives and/or are incorrectly annotated and/or classified. For example, some reads may initially appear to align to a particular classified reference genome using a first alignment algorithm (e.g., a global alignment algorithm like, for example, SNAP or Needleman & Wunsch) but a more comprehensive alignment (e.g., using a local alignment algorithm like, for example, “Basic Local Alignment Search Tool” (BLASTn) or Smith & Waterman) shows that even a best read can actually align to a different reference sequence using different alignment algorithms.

Furthermore, some reference genomes with a GenBank or other large database of classified reference genomes may be mis-annotated taxonomically (i.e., are classified incorrectly). Additionally, some reference genomes may include portions of sequences that can be assigned to multiple taxonomies. If a sequence in the dataset aligns to the “erroneous” portion of that reference genome sequence, it will be mis-assigned by an alignment process. For example, an HIV viral sequence with flanking human integration sites may be annotated as HIV and may cause erroneous matches and identification of potential pathogens. As such, embodiments are directed at filtering these “erroneous” reference genomes using a filtering algorithm.

Accordingly, embodiments of the present invention can be applied after a sequencing alignment analysis to help filter, categorize, and interpret the results of the sequencing analysis to identify a smaller list of potential pathogens to a more manageable and easily interpretable result. Thus, embodiments may be directed at filtering results once a system has the raw reads, raw contigs, and assembled contigs are determined. Embodiments can take the reads and annotate them with what the hits were from the alignment analysis. For example, the system may take the results of the analysis and determine how accurate and/or good of a hit a result is. The system can take that annotated raw output and further filter it to remove mis-annotated results and false positive aligned reference genomes. Generated result tables showing the alignments and the coverage maps can be determined and displayed after filtering.

A. Filtering Method

FIG. 5 shows an exemplary filtering process 500 for removing false-positive and/or other erroneous reference genome sequences from alignment results. Before the process shown in FIG. 5, the system may receive a plurality of sequence reads obtained from a sequencing of DNA molecules from a sample of biological material that includes multiple different organisms. In some embodiments, the alignment results may have been previously taxonomically classified, e.g., using process 400.

At step 510, alignment results are obtained in response to applying an initial alignment technique (e.g., using a global alignment algorithm) to a plurality of sequence reads from a sample to align the sequence reads to reference genomes in a global database. The alignment results include a matching reference genome for each of at least a portion of the sequence reads in the sample that the sequence read aligns. For example, the alignment technique may include a global alignment algorithm that takes each sequence read and searches a database of classified reference genomes for aligned portions of the sequence reads to portions of the classified reference genomes. The results may include multiple aligned reference genomes for each of the reads, and each reference genome sequence may have many different reads that align to a portion of the reference genome.

The alignment results may include a plurality of sequence reads, sequence read identifiers, a plurality of aligned classified reference genomes, reference genome identifiers from the database of reference genomes, an alignment and/or similarity measurement for each of the alignments, a taxonomy identifier and/or other classification information associated with each of the reference genomes, and any other suitable information that may be used in alignment and identification processes. For example, the source and/or database identifier where the aligned reference genome was stored, the source where the reference genome was provided from (e.g., hospital, clinic, company, provider, etc.), and the data and/or time of the upload may also be provided.

FIG. 6 shows an exemplary summary 620 of sequence alignment results for three different samples. Column 610 shows read type. The first sample 630 and second sample 640 are known to be entrovirus D68, and the third sample 650 is a control that aligns with hepatitis C. The three samples shown in FIG. 6 include a result table that shows the number of initial sequence reads (e.g., for sample 1 there are 2,159,196 sequence reads), the number of sequence reads after preprocessing (e.g., for sample 1 the preprocessing results in 1,933,046 sequence reads), the number of sequence reads that align to human reference sequence genomes (e.g., for sample 1, 1,725,425 reads align to human reference sequence genomes), the number of reads that matched to a reference genome (e.g., for sample 1, 133,932 reads aligned with a reference genome from the database of reference genomes), and the number of bacteria and viruses matched (e.g., 133,827 and 21, respectively for sample 1). Further, the sequence statistics shows the percentage of human matches (e.g., 89.3% for sample 1) and the percentage of reads that remained after preprocessing (e.g., 89.5% for sample 1).

Additionally, FIG. 7 shows an exemplary result table 700 for initial sequence reads. The result table includes the multiple reference genomes that matched with at least one sequence read, a taxonomic identifier for each reference genome including a species (710), genus (720), and family (730) classification for each reference genome. The species, genus, and family classifications are hierarchical with the species level being the lowest level, the genus level being above the species level, and the family level being above both the species and the genus level. The taxonomic identifier and the corresponding classification information may be provided by the classification information (e.g., an annotation entry) of an aligned reference genome from the database that was used during alignment.

Further, the results table may include tag (740) for each reference genome that may identify the database that the reference genome was aligned with and the type of reference genome (e.g., human, plant, bacteria, etc.). Finally, the result table may include the number of reads that aligned to each reference genome for each sample. For example, for the first sample (column 750) shown in FIG. 7, there are 133,004 reads aligned to Enterovirus D; 247 reads aligned to Enterovirus sp., 574 sequence reads aligned to Human Enterovirus, 1 read aligned to human enterovirus K1105_171204, and 1 read aligned to human rhinovirus sp. Accordingly, initial alignment results may include any alignment matches to any of the reference genomes that are present in the database. The second sample (column 760) has a similar pattern as for the first sample 750. The negative control (column 770) shows a match for hepatitis C, as is expected.

At step 520, the system identifies all reference genomes that matched to a read from the sample. For example, for the alignment results of sample #1 shown in FIG. 7, five different reference genomes were aligned with sequence reads from the sample.

At step 530, the system identifies an optimally-aligning sequence read that aligns to the matching reference genome with an optimal alignment score that exceeds or is equal to alignments scores of other sequence reads that align to the matching reference genome. The optimally-aligning sequence read and optimal alignment score can be identified based on the initial alignment results. The optimally-aligned sequence read may be determined by analyzing the reads that aligned with the reference genome and comparing their alignment value for the closest match. Any process may be used to identify the best match, for example, the read with the most significant match may be defined in SNAP by minimum edit distance and RAPSearch by minimum e-value.

At step 540, the system applies a second alignment technique for the optimally aligning sequence read to the plurality of classified reference genomes to obtain a plurality of new alignment scores. The second alignment technique may include a different alignment technique, a different reference database, and/or a combination thereof. For example, a local alignment algorithm may be applied for the second alignment technique that may take longer to process for each sequence read but because there are fewer sequence reads to process (e.g., only the identified optimally-aligning sequence reads from step 530), the processing time may not overly delay the identification process.

For example, the first alignment technique may include a global alignment algorithm (e.g., SNAP/RAPSearch) that provides faster analysis, but may be sensitive to poor quality sequence regions. The second alignment technique may include applying a local aligner (e.g., BLASTn, Bowtie2, etc.) that may be much slower (1,000×-10000× slower), but that is less sensitive to poor reads and may provide a different alignment of the sequence read that can be compared to the results of the first alignment technique. Thus, the use of the two different alignment techniques allows the system to correct rare “false positives” that occur from the first alignment technique (e.g., SNAP alignment). By only selecting one “test hit” read per reference genome for the second alignment technique (e.g., BLASTn confirmatory alignment), the system can greatly speed up the turnaround time of the process.

At step 550, the system compares the re-alignment results obtained using the second alignment technique for the optimally-aligned read with the initial alignment results obtained using the first alignment technique. For example, different reference genomes may be aligned with the sequence read using the second alignment technique that were not aligned and/or different alignment scores may be returned with the re-alignment than the original alignment.

At step 560, the system determines whether the alignment results using the first alignment technique match the alignment results using the second alignment technique. The comparison of the first alignment results and the second alignment results to determine a match may be accomplished through any suitable manner. For example, the system may determine whether any of the new alignment scores exceeds the optimal alignment score for the optimally aligning sequence read. In other embodiments, the system may determine whether the same reference genome was aligned with the sequence read. Further, other embodiments may determine whether the aligned reference genome sequences share the same taxonomic identifier (e.g., the two reference genome sequences share the same family, genus, and species classifications).

At step 570, if the alignment results for the sequence read are determined not to match, the system may remove the reference genome sequence from the results and all corresponding sequence reads that aligned with the reference genome sequence. In some embodiments, if a sequence read aligned with multiple different reference genomes, only the association of that read with the currently analyzed reference genome sequence may be removed and the read itself (along with the alignment result to the other reference genome sequence) may not be removed. Thus, for the exemplary comparison method of comparing the alignment scores discussed above, when one new alignment score exceeds the optimal alignment score for the optimally aligning sequence read, the system may remove the matching reference genome from the set of matching reference genomes and the corresponding sequence reads associated with the reference genome. In some embodiments, if the two reference genomes share a same taxonomic identifier, the matching reference is not removed.

For example, FIGS. 7-8 shows a table of example results where filtering has been applied to remove mis-annotated and false positive alignments. FIG. 7 shows the raw initial alignment results and FIG. 8 shows the results 800 after filtering has been applied. As shown in FIG. 8, the effect of filtering is to remove two mis-assigned reads that occur as a result of misannotations in the database (761 and 771). These reads are annotated as “hepatitis B” and “simian virus 40” but actually are human integration sites that have been mis-annotated as viral. Using a local aligner (e.g., BLASTn), the system determined that these two reads are actually human reads (returning human reference genomes), and are therefore removed from the output when the alignment results do not match with the initial alignment results of “hepatitis B” and “simian virus 40”.

At step 580, if the alignment results for the sequence read are determined to match, the system may consider the reference genome a real match that is not mis-annotated and may maintain the reference genome and aligned sequence reads in the alignment results.

At step 590, the system determines whether all of the reference genomes have had the filtering process applied. If there are additional reference genomes to be analyzed, the system may repeat the process for the next aligned reference genome in the initial alignment results (and may return to step 530). This process may continue until all of the aligned reference genomes have been analyzed for false-positive matches and are removed and/or confirmed as matches.

At step 595, once the system has applied the filtering process to all of the aligned reference genomes, the process may provide the updated alignment results (e.g., not including the selected reference genomes from step 570) to the sample alignment system for display to an operator/clinician. In some embodiments, just the set of matching reference genomes is output. In some embodiments, additional processing may be applied to the updated alignment results before returning the updated alignment results to the sample analysis computer, as described herein.

Accordingly, embodiments may apply a filtering algorithm to the results that identifies all the false positive hits, and systematically remove false positive hits corresponding to that genome. The system can identify the particular genomes that provided a sufficient quality alignment. The system can then further analyze those reads to remove further alignments that do not appear to be high quality. The filtering algorithm can apply an even more stringent criteria to remove hits that are likely not real and instead are due to a misannotation of a reference genome database and/or a misclassification.

For example, a particular read might have hit two, five, or more hits. If a genome had a single hit then it is likely an accurate hit, but if a read hit ten different reference genomes, then it is not likely that the hit is of a high quality. Thus, the filtering could remove all those hits and a new summary table and/or coverage map for the genome. Further, if the coverage map for a particular genome has now drastically reduced due to removing the false positive hits, then the system may remove the genome as a potential pathogen or cause of the illness. Accordingly, embodiments of the present invention may filter what was fifty results (e.g., fifty different bacteria) down to the five bacteria that are the most high quality and likely real results. This allows for much easier and faster clinical analysis and/or interpretation.

In some embodiments, methods 400 and 500 can be combined in a single pipeline, e.g., in a research pipeline. In other embodiments, method 500 may not be performed, e.g., a clinical pipeline might not include method 500. For example, the taxonomic classification may be sufficiently accurate that filtering may not be needed. And, such a pipeline can operate faster without the filtering step, particularly one that uses BLASTn. The turnaround time for a clinical pipeline may be desired within 4-6 hours, whereas more time is acceptable for a research pipeline. Thus, filtering, RAPSearch translated nucleotide alignment, and contig assembly may be retained in a research pipeline, but dropped for a clinical pipeline.

B. Ribosomal RNA Sequence Removal

Additionally, in some embodiments, ribosomal RNA sequences may be removed from the alignment results. Sequences derived from ribosomal regions are difficult to speciate/classify accurately due to conservation of the ribosomal sequences. Accordingly, embodiments may remove ribosomal sequences from the alignment output by aligning sequence reads that had been assigned to bacteria against bacterial ribosomal 16S/23S reference genome sequences. Additionally, the aligned sequence reads assigned to eukaryotes (fungi and parasites) may be further aligned to 18S/28S GIs. Any sequence reads that align to these bacterial or eukaryotic ribosomal reference genomes may be removed from the alignment results.

C. “On-the-Fly” Annotation of Database of Classified Reference Genomes: Automated Clean-Up of the Reference Database

Further, due to the rapidly expanding number of reference genomes stored in the classified reference genomic sequence database, such databases may double every 3 years or faster. However, there are very few controls over who or the quality of those providing reference genomic sequences. Accordingly, embodiments can includes automated scripts to (1) download all sequences from the GenBank (currently ˜72 GB of data), (2) remove all entries that are poorly annotated by screening out key words in the description of a reference genome sequence entry (e.g., “uncultured”, “unclassified”, “environmental”), (3) remove all vector sequences, and (4) chop up the classified reference genome sequence database (e.g., GenBank) into “chunks” and convert them into reference databases to be used by local and global alignment algorithms (e.g., SNAP and RAPSearch). A database may be chopped up so as to fit into memory.

V. Removal of Background

As mentioned above, a negative control can be used to identify contaminants or other microorganisms that are not clinically relevant. Such remove of background microorganisms is particularly relevant to clinical applications, where diagnosis and treatment of a patient are a goal. Embodiments can provide one or more criteria for discriminating between a pathogen and a contaminant.

In some embodiments, one or more control samples can be analyzed in parallel (simultaneously) with one or more patient samples (e.g., 5-8 patient samples). A positive control can consist of several different organisms spiked into negative matrix, e.g., spinal cord fluid. A negative control sample (also a called a no-template control) has no spiked organisms. The negative control can be just the buffer that is used for a PCR reaction, which is part of the library preparation for sequencing.

A. Types of Controls

The positive control can be used to confirm that the microorganisms in the spiked contamination are identified. The positive control can ensure that the procedure is robust. For instance, if the sample is like a bloody fluid, the heme in blood may actually inhibit PCR. There can be a negative result, but it might be a false negative. Thus, if the positive control is negative, then the sample can be identified as a false negative. The external positive control can ensure that there are a certain number of reads for each of the spiked organisms (e.g., 7 different organisms), in order for the run to pass quality control.

With the no template control (e.g., a buffer for a PCR reaction or other library preparation step), any reads classified (i.e., hits) as aligning to a microorganism can be identified as corresponding to background. As examples, such false hits can correspond to: reagent contamination, contamination introduced in a laboratory, or contamination introduced from other samples (cross-contamination). In this manner, the number of hits can be reduced, and proper treatment for an actual pathogen is more likely to be performed.

An internal spiked control can include organisms spiked into a patient sample. For example, a DNA phage and/or an RNA phage can be spiked in specific concentrations in a clinical sample. Embodiments can check to make sure that the phages can be detected in the DNA library, e.g., a sufficient number of sequences from the DNA phage. Similarly, a check can make sure that the phages are detected in the RNA library, e.g., a sufficient number of sequences from the RNA phage. This control can be an internal control on every sample, in addition to external control.

B. Normalization

After analyzing the data from the sequencing of the samples (e.g., aligning, classifying, filtering), a list of potential candidate pathogens can be obtained. The list can be a table showing a number of reads seen for each of a plurality of species or other taxonomic classification. The number of reads can be normalized. For instance, the number of reads can be determined per million reads. If there are 10 hits and 10 million reads, the RPM (reads per million) will be one. The total reads can be raw reads that were generated per clinical sample. But, there can be further normalization.

An RPM can also be determined in the negative control. An RPM ratio can be determined from the reads per million (i.e., of a particular taxonomic classification) in the sample divided by your reads per million in the negative control. For example, assume there are 10 normalized reads for a virus HIV) in the sample, and two normalized reads in the no-template control. In that case, the RPM ratio would be five. If there are no reads in the negative control, then there would be a division by zero error. In such cases, a value of one can be used as the RPM from the NTC sample. In other words, the RPM ratio can be the reads per million in the sample.

The RPM ratio allows for easier discrimination of what is really in the sample and what is background. A threshold (e.g., 10) can be used for the discrimination. Thus, pathogens with an RPM ratio of greater than 10 can be identified as clinically significant. Using this criteria, very good sensitivity and specificity can be obtained. The use of the RPM ratio reduces the number of false positive. In particular, very good test performance is obtained relative to using the Bayesian probability, which is based on just the number of reads in a sample. Such probabilistic method can be used when determining an error model, and do not involve normalization from reads in a negative control.

Tables 1-4 below shows revised accuracy data with and without discrepancy testing for results-based and sample-based testing. Discrepancy testing refers to the use of an orthogonal (i.e. different kind of) test to resolve discrepant results between the mNGS assay and conventional clinical microbiology laboratory testing. For instance, if a sample is found to be positive for Mycobacterium tuberculosis by mNGS but is negative by culture (because some Mycobacterium tuberculosis strains grow slowly or not at all in culture depending on pathogen titer and sample type), a Mycobacterium tuberculosis PCR test can be used as an orthogonal test for discrepancy testing, and the results of the orthogonal test are taken as correct for the purposes of determining accuracy.

For sample-based testing, for each sample, performance is evaluated only as it pertains to detection of the original organism reported by the reference lab. For samples completely negative by reference lab testing, the ability to detect all 5 organism types is evaluated.

For result-based testing, for each sample, performance is evaluated as it pertains to detection of all 5 organism types (bacteria, fungi, DNA virus, RNA virus and parasites). Only results for acceptable samples (passed quality control metrics of >5,000,000 reads per library and >10 RPM of either an internally spiked DNA control, T1 bacteriophage, for DNA libraries, or an internally spiked RNA control, M2 bacteriophage, for RNA libraries; sufficient sample volume) are provided in the tables below.

Table 1 shows sample-based accuracy (most clinically significant result per sample) without discrepancy testing.

TABLE 1 Sample-based accuracy (most clinically significant result per sample) ClinMicro Reference Lab Total Pos Neg Untested Total mNGS Pos 59 7 7 73 Neg 25 90 115 Total 84 97 7 Sensitivity TP/(TP + FN) 70.2 Specificity TN/(TN + FP) 92.8 Accuracy (TP + TN)/ 82.3 (TN + FN + FP + FN)

Table 2 shows results-based accuracy (includes multiple positive i negative results per sample) without discrepancy testing.

TABLE 2 Results-based accuracy (includes multiple positive/negative results per sample) ClinMicro Reference Lab Total Pos Neg Untested Total mNGS Pos 60 7 31 98 Neg 27 404 431 Total 87 411 31 Sensitivity TP/(TP + FN) 69.0 Specificity TN/(TN + FP) 98.3 Accuracy (TP + TN)/ 93.2 (TN + FN + FP + FN)

Table 3 shows revised sample-based accuracy after discrepancy testing.

TABLE 3 Revised accuracy after discrepancy testing, sample-based ClinMicro Reference Lab Total Pos Neg Untested Total mNGS Pos 68 7 22 97 Neg 11 327 338 Total 79 334 22 Sensitivity TPR TP/(TP + FN) 86.1 Specificity TNR TN/(TN + FP) 97.9

Table 4 shows revised results-based accuracy after discrepancy testing. Only results for acceptable samples are provided.

TABLE 4 Revised accuracy after discrepancy testing, results-based ClinMicro Reference Lab Total Pos Neg Untested Total mNGS Pos 47 6 18 71 Neg 8 303 311 Total 55 309 18 Sensitivity TPR TP/(TP + FN) 85.5 Specificity TNR TN/(TN + FP) 98.1

After discrepancy testing and including results passing IC with sufficient sample volume, overall sensitivity increases to 86.1% (sample-based) and 85.5% (results-based), while maintaining high specificity. Given the difficulty making an etiologic diagnosis in many cases of encephalitis/meningitis, a test sensitivity of 80-90% with specificity >95% is very useful in patient management decisions. He above data shows that the overall test accuracy is acceptable for clinical use

The tables can be displayed in a visualization program. The tables can be provided in a user interface, in a normalized or non-normalized manner. A normalized table is easier to identify pathogens than to look at a list of hundreds of hits. For example, many hits can be removed after normalization to obtain the RPM ratio, and use of a threshold.

C. No-Template Control Database

In addition to the use of the no-template control for a given experimental run, embodiments can use a database of no-template controls (NTC database). For every run, the database can be reviewed before a call is made about a pathogen being in the sample, e.g., to make sure that organism is not in the NTC database. The NTC database includes historical data about what pathogens were identified in NTC samples, e.g., within a specific amount of time, such as 30 days.

For example, a herpes virus might be identified in an NTC sample a month ago (e.g., due to contamination resulting from the same lab doing herpes virus testing). Herpes virus 1 is a common cause of meningitis, and thus can be a dangerous pathogen, if it actually existed. The NTC database can include herpes virus 1, until it is not seen in an NTC sample for at least 30 days. In some embodiments, if the taxonomic classification is found on the skin like poly papilloma viruses, they can be excluded outright. But, they can still be included as part of the NTC database.

Very often, background includes skin flora. For instance, papilloma viruses may be seen in clinical samples in the no-template control. If we do see such a virus in a no-template control, it is added to the database. Embodiments may not actually report organisms that were found in the no-template control unless they are present at a much higher level. For example, assume the RPM was 10 in the NTC for one run, then the taxonomic classification gets added to the database with an RPM of 10. For the next run, if there was a herpes virus detected at 10 RPM (i.e., an RPM ratio of 1), it would not be reported when a threshold of greater than 10 for the RPM ration was used. It would have to be at an RPM of 100 to provide an RPM ration of 10 times the highest level that was present in the NTC database.

As it is not desirable to keep these pathogens in the database because that would affect the performance of the test, as the NTC database decreases the sensitivity, since the bar is set higher to make a positive call for the pathogens in the NTC database. Thus, the NTC database is an evolving database. For example, if a pathogen has not been present for longer than three months (or other specified time period), it then gets removed from the database. Steps can be performed to remove contaminants from NTC samples, e.g., replacing reagents. The amount (e.g., RPM) in the NTC database for a particular contaminant can be updated from the initial value to a new value if the new value is greater than the initial value. Accordingly, in one embodiment, the amount in the NTC database for a contaminant can be the maximum that is observed.

D. Flowchart for Use of NTC

FIG. 11 is a flowchart of a method 1100 for identifying pathogens in a test sample of biological material according to embodiments of the present invention. Method 1100 may be performed by a computer system. The test sample includes DNA molecules from a plurality of organisms, e.g., from a patient and from microorganisms.

At block 1110, a plurality of sequence reads are obtain from a sequencing of DNA molecules from the test sample of biological material. The sequence reads can be received at the computer system in a variety of ways, e.g., via a network or a removable storage device. The computer system can also receive information about sequence reads from a negative control sample. For example, the computer system can receive an amount of sequence reads for a plurality of difference reference genomes, e.g., as classified into various taxonomic levels.

At block 1120, an alignment technique is used to align the plurality of sequence reads to a plurality of classified reference genomes in a database. The reference genomes may be considered classified in that the genomes are classified to at least one classification level in a taxonomy. Alignment results can include, for each of at least a portion of the sequence reads, at least one matching reference genome to which the sequence read aligns. The alignment results can include classification information for each of the matching reference genomes.

Blocks 1130-1160 may be performed for each of a group of one or more matching reference genomes.

At block 1130, a first amount of sequence reads from the test sample that align to the matching genome is determined. As an example, the first amount may be a total number of sequence reads aligned to the particular matching reference genome. The first amount may be normalized based on the number of sequence reads obtained from the test sample, or the total number of sequence reads aligned to at least one reference genome

At block 1140, a second amount of sequence reads from a negative control sample that align to the matching genome is determined. In one embodiment, the second amount of sequence reads can be determined from a negative control database. In another embodiment, the second amount of sequence reads can be determined from a parallel sequencing of a negative control sample. Amounts from a parallel negative control sample and the negative control sample can be combined to provide a single amount.

At block 1150, a ratio of the first amount and the second amount is determined. The ratio can take various forms, such as the first mount divided by the second amount, or the second amount divided by the first amount. Further, a numerator or denominator can include a sum of the two amounts.

At block 1160, the ratio is compared to a threshold to determine whether the ratio exceeds the threshold. A set of one or more matching reference genomes can have the ratio exceed the threshold. The threshold can be determined based on empirical data (e.g., values of samples with a known composition), as will be appreciated by one skilled in the art. The selection of the threshold can be based on a desired balance between sensitivity and specificity.

At block 1170, an output is provided that identifies the set of one or more snatching reference genomes as potential pathogens in the test sample. As examples, the output can be a list of the set of matching reference genomes, e.g., presented in classification levels. Such a list could include more matching reference genomes, where the set is indicated, as may be done with a marker or the ratio. For instance, an RPM ratio can be provided.

VI. Further Data Analysis

Once the false-positives and inaccurate reads have been removed, the system may apply post-filtering analysis and provide the results to the sample analysis system for display the clinician. In some embodiments, the filtered and classified results may be used to generate new coverage maps. The filtered and classified results may be analyzed to determine a best match and whether a sufficient match obtained. For example, a best match score (or any one of the individual score) can be required to be above a specified threshold. In some embodiments, if the coverage is below a certain level, than the system will dismiss the organism as being a possibility. In other embodiments, all possible identification results could be provided to allow a clinician to make an identification.

Further, in some embodiments, a best match algorithm may be applied to select, for instance, the best match among a variety of closely related matches. For instance, the results may have hits to multiple different strains of influenza. The best match algorithm may choose the most closely matched sequence in the reference database and select that classified element as the best match for the sample sequence.

Examples of a best match algorithm are described herein, e.g., in section II.C. In another example, for each species, the reference genome sequence with the greatest number of reads assigned to it can be used as a reference for mapping of all of the species-specific reads. A “consensus sequence” may be generated from the mapping. The global pairwise identity using the Needleman-Wunsch algorithm may then calculated between the “consensus sequence” and each reference genome sequence. The reference genome sequence with the greatest global pairwise identity may be selected as the top reference genome sequence. The reference genome sequences under consideration can also be prioritized as follows: (1) complete genomes, (2) complete sequences, or (3) partial sequences/individual genes.”

A positive pathogen identification may be reported where for each top reference genome a positive identification metric may be evaluated and provided to the clinician in a return message to the sample analysis computer. For example, the positive identification metric may include coverage of at least 3 genes/open reading frames (ORFs), or where there are fewer than 3 ORFs in the virus, then coverage that exceeds greater than 3 times the read length (e.g., for 100 bp sequencing, coverage of greater than 300 bp).

Additionally, systems may implement a tagging mechanism for the reference genomes that are matched to the sample sequence reads in order to make pathogen identification easier. For instance, the system may add the “host” to viral reference genome sequences (GIs), where the “host” is bacteria, for instance, the system knows that the virus is a phage (which the system may want to mask out sequences to phages). Similarly, the system can apply arbitrary labels such as “pathogens”, “colonizers”, and “laboratory contaminants” to specific bacterial GIs that can be helpful in downstream visualization of results for clinical interpretation. These results may be provided with the final results of the alignment, filtering, and classification process.

VII. Example for Neurobrucellosis

A diagnosis of brucellosis can be difficult because routine culture and serological methods exhibit variable sensitivity and specificity. At present, the standard laboratory diagnosis of brucellosis is based on isolation of the bacteria from clinical specimens and/or serological detection of Brucella antibodies.

Neurobrucellosis is a known complication of systemic Brucella infection. It remains a difficult diagnosis to make and can mimic other fastidious infections such as TB [1]. Manifestations of neurobrucellosis are widely variable and include meningoencephalitis, cerebrovascular disease, peripheral and cranial neuropathies, or myelitis. Neuroimaging studies of neurobrucellosis can also vary greatly among individuals [6]. A previous study from Turkey reported a high prevalence of neurobrucellosis at 37% among 128 patients with brucellosis. Although culture is the gold standard for diagnosis, Brucella species are relatively fastidious and slow-growing; cultures fail to recover the organism 30-90% of the time, and also present a risk of laboratory-acquired infection [7]. Serology is more sensitive for detection, but can lead to false-positive results and may not distinguish between active and prior infection.

Molecular methods based on detection of nucleic acid such as PCR [8] and now potentially metagenomic NGS [13] can offer increased sensitivity and specificity over conventional diagnostic testing. We present the use of a metagenomic next-generation sequencing assay to diagnose a case of neurobrucellosis from cerebrospinal fluid (CSF), resulting in the institution of appropriate antibiotic treatment and a favorable clinical outcome.

A. Case Report

Initially, microbiological testing was positive for Epstein-Barr virus (EBV) and human herpesvirus 7 (HHV-7) from CSF. Testing for other pathogens, including Brucella by IgM antibody, was negative. The patient was treated with parenteral acyclovir, followed by oral ganciclovir to complete a 14-day course. Two weeks after hospital discharge, she developed back pain and worsening headache. Upon hospitalization, the patient's vital signs were normal without fever. Physical examination was remarkable for multifocal myoclonus throughout her body. There was hyperreflexia and decreased proprioceptive sensation of bilateral lower limbs. A complete blood count was notable for a white blood cell count of 4.82×10⁹/L (44% neutrophils, 44% lymphocytes, 10% monocytes).

A repeat lumbar puncture was performed, with additional microbiological testing unrevealing. Cytologic examination of the CSF revealed no malignant cells. Despite a negative tuberculin skin test and QuantiFERON-TB test, she was started on 4-drug therapy with isoniazid, rifampin, pyrazinamide, and ethambutol given a high concern for TB disease based on the patient's suggestive CSF profile and risk factors for TB.

On day 8 of hospitalization, the result of CSF for Mycobacterium tuberculosis (MTB) by PCR was negative. Given the lack of response to empiric TB therapy, there was a concern for drug-resistant TB. Ethambutol was thus changed to ethionamide, and levofloxacin was added to the regimen. She improved substantially after changing her antibiotic regimen, and was discharged home on 5 anti-TB medications. At follow up 1 week and 1 month after discharge, her headache had resolved, but she continued to have fatigue, mild back pain, and intermittent episodes of shaking of her extremities. When the result of the mycobacterial culture was finalized as negative at 6 weeks, she returned to Mexico to continue her anti-TB therapy with INH and rifampin alone.

B. Results

RNA and DNA mNGS libraries from 600 μL of the patient's CSF sample were constructed. From a total of 23,638,587 raw reads in the DNA library, there were 277 (0.0012%) reads corresponding to the Brucella genus in the DNA library, corresponding to an RPM ratio of 15.6, with all species-specific reads aligning to Brucella melitensis (see Tables 1200 of FIG. 12, 1300 of FIG. 13, and 1 and 2 of the Appendix). Notably, no reads aligning to Brucella among 9,161,626 reads in the corresponding RNA library from CSF were detected, and Brucella reads were absent in both negative “no template control” (NTC) and positive control samples. Other bacterial reads present in DNA library were also present in the NTC sample processed in parallel on the same run, and thus did not meet the established threshold for reporting using an RPM ratio of ≥10 and were attributed to laboratory reagent contamination.

Similarly, no viruses, fungi, or parasites met criteria for reporting. In contrast, all 7 microorganisms in the positive control were detected at levels above the established reporting threshold. Importantly, no Brucella isolates or positive clinical samples from suspected or confirmed cases had been present in the clinical laboratory prior to the mNGS testing, nor had Brucella sequences ever been detected in the NTC. The presence of Brucella in the patient's cerebrospinal fluid was confirmed by Brucella-specific PCR testing and Sanger sequencing of the amplicon. The PCR reaction for 177 bp region of Brucella IS711 gene was run on a 2% agarose gel. Other samples included: a no-template-control; 7 organism positive control mixture of CMV, HIV, Streptococcus agalactiae, Klebsiella pneumoniae, Cryptococcus neoformans; water; and 1 kb Plus DNA ladder (Invitrogen).

After detection of Brucella reads in her CSF sample, the patient was contacted and instructed to seek further medical evaluation. She returned to the hospital in Los Angeles for a second admission. Although she had completed INH and rifampin therapy one week prior, she reported persistent back pain, nausea, and fatigue. Repeat MRI of the brain and spine was remarkable, and CSF was normal. A confirmatory serum Brucella agglutinin titer was positive at 1:80. Because of persistent symptoms, mNGS and PCR testing showing Brucella, and positive confirmatory serology, she was diagnosed with chronic neurobrucellosis and started on targeted therapy with doxycycline and rifampin. Two weeks after starting therapy, she reported that her symptoms had fully resolved. Notably, Brucella spp. was not isolated from either CSF or blood culture, nor was Brucella DNA detected in CSF submitted for 16S rRNA gene sequencing.

This highlights the clinical impact of mNGS for diagnosis of infections such as neurobrucellosis that are challenging to accurately diagnose and treat. Although initially covered for TB meningitis using antibiotics that partially treat neurobrucellosis (rifampin+/−levofloxacin), the patient continued to be symptomatic and was not placed on targeted therapy for Brucella until comprehensive metagenomic sequencing revealed the presence of bacterial DNA in her CSF and neurobrucellosis was confirmed by Brucella agglutinin testing.

FIG. 12 shows table 1200 of sequencing statistics according to embodiments of the present invention. Results are shown from both the patient's sample, no template control (NTC), and the positive batch control (PC) sample.

FIG. 13 shows a table 1300 of the number of reads from DNA libraries and RNA libraries aligning to viral sequences according to embodiments of the present invention. Table 1300 shows the number of reads from DNA libraries and RNA libraries aligning to viral sequences. Results are shown from both the patient's sample, NTC, and the PC sample.

In table 1300, the reported viruses are cytomegalovirus (CMV) in the “DNA PC” column and human immunodeficiency virus 1 (HIV-1) in the “RNA PC” column. These are the viruses spiked into the PC sample. No viruses are reported from the patient CSF because human papillomaviruses, part of skin flora and presumed to be a contaminant, are not reported by the mNGS assay.

Supplementary Table 1 of the Appendix shows the number of reads from the DNA library aligning to bacterial sequences. Results are shown from both the patient's sample, NTC, and the PC sample. In Supplementary Table 1, “*” indicates re-classified to a higher taxonomic rank because reads aligned equally well to multiple different organisms that shared the same species, genus or family. The abbreviations are as follows: NTC, no template control; PC, positive control; and CSF, cerebrospinal fluid.

Supplementary Table 2 of the Appendix shows the number of reads from the DNA library aligning to fungal and parasite sequences. Results are shown from both the patient's sample, NTC, and the PC sample. In Supplementary Table 2, “*” indicates re-classified to a higher taxonomic rank because reads aligned equally well to multiple different organisms that shared the same species, genus or family.

Supplementary Table 3 of the Appendix shows bacterial taxa identified by embodiments using an RPM ratio metric for reads shown in Supplementary Table 1. In Supplementary Table 3, “*” indicates re-classified to a higher taxonomic rank because reads aligned equally well to multiple different organisms that shared the same species, genus or family; “#” indicates positive result according to an RPM (reads per million) ratio≥10, where the RPM ratio=RPM(sample)/RPM(NTC); and “&” indicates reads correspond to the Streptococcus agalactiae PC. The abbreviations are as follows: PC, positive control; CSF, cerebrospinal fluid; mNGS, metagenomic next-generation sequencing; and RPM, reads per million.

In Supplementary Table 3, cells that show a positive result by mNGS assay are highlighted in yellow. Under the column designated “DNA PC”, the two bacterial species reported at the pre-established threshold criterion of >10 RPM are Streptococcus agalactiae and Klebsiella pneumoniae, which were spiked into the PC sample. Other positive rows are not reported because they represent either (1) classification at a higher taxonomic rank (genus) or (2) another bacterial species in the same genus (i.e. Streptococcus suis) that does not meet the pre-established requirement of having at least 1/10 of the number of RPM corresponding to the predominant species (Streptococcus agalactiae) to call a co-infection from at least 2 different species. Under the column designated “DNA Patient CSF”, the only bacterial taxa with >10 RPM is Brucella genus, which is the reported organism reported in the patient's CSF.

Supplementary Table 4 of the Appendix shows fungal and parasitic taxa identified by embodiments using an RPM ratio metric for reads shown in Supplementary Table 1. In Supplementary Table 4, “*” indicates re-classified to a higher taxonomic rank because reads aligned equally well to multiple different organisms that shared the same species, genus or family; and “#” indicates positive result according to an RPM ratio≥10, where the RPM ratio=RPM(sample)/RPM(NTC). The abbreviations are: PC, positive control; CSF, cerebrospinal fluid; mNGS, metagenomic next-generation sequencing; RPM, reads per million.

In Supplementary Table 4, cells that show a positive result by mNGS assay are highlighted in yellow. Under the column “DNA PC”, the two fungal/parasitic species reported at the pre-established threshold criterion of >10 RPM are Cryptococcus neoformans, Aspergillus niger, and Toxoplasma gondii, which were spiked into the PC sample. Other positive rows are not reported because they represent classification at a higher taxonomic rank (genus or family).

Supplementary Tables 3 and 4, which are used to interpret and report mNGS results, show the effect of taxonomic classification and normalization using an “RPM ratio” metric in comparison to the NTC in simplifying the clinical interpretation, as compared to Supplementary Tables 1 and 2, respectively.

Of particular relevance for this case, two chronic granulomatous infections, tuberculosis and brucellosis, have not only overlapping clinical or radiographic features but also histologic characteristics in common [10]. As such, misdiagnosis of tuberculosis in patients with brucellosis has been reported in the literature [10]. This matter is further complicated by the fact that neither a negative CSF mycobacterial culture nor tuberculosis PCR-based assay excludes the diagnosis of TB meningitis if the clinical suspicion is high. Additionally, false-positive Brucella seroreactivity, ELISA and agglutination titer, in patients with active TB have also been reported. In this patient's case, the negative Brucella IgM but positive IgG was incorrectly attributed to false-positive Brucella seroreactivity in the setting of TB and not to active Brucella infection. It is possible that the patient may not have mounted a detectable IgM antibody response, or that Brucella IgM levels had waned by the time of hospital admission. Confirmatory agglutinin testing may have been helpful in making the diagnosis of Brucella earlier.

This patient had temporary clinical improvement after initiation of anti-TB medications. Rifampin and levofloxacin are two anti-TB agents that are also active against Brucella spp. [1]. However, rifampin, while active against Brucella, should always be used in combination with other agents (e.g. doxycycline, trimethoprim-sulfamethoxazole, or quinolones) since monotherapy, as was inadvertently administered to this patient initially, has been associated with high relapse rates [1]. In hindsight, the patient's prolonged and indolent course was more likely to be associated with Brucella, since more rapid clinical deterioration would have been expected for patients with active TB who are inadequately treated [11]. But, Brucella at that time was not high on the differential. It is only after discharge and because of the patient's persistent symptoms that an alternative diagnosis was considered. Taken altogether, knowledge of these pitfalls is essential for clinicians to reduce diagnostic errors.

Metagenomic next-generation sequencing (mNGS) is an emerging approach in diagnostic microbiology with the ability to detect all microorganisms—viruses, bacteria, fungi, and parasites—in a single assay [2, 4, 9, 2, 3]. Here mNGS was used to provide an accurate diagnosis of neurobrucellosis and to guide the institution of targeted therapy, leading to complete resolution of the patient's illness.

C. Implementations Details

1. Metagenomic Library Construction

DNA and RNA metagenomic libraries were constructed from the patient's CSF sample as previously described [2, 3]. After bead-beating using Lysis matrix B (MP Biomedicals, Santa Ana Calif.) at 6 m/s for 30 seconds, total nucleic acid was extracted using the Qiagen EZ1 Viral kit (Qiagen, Valencia Calif.). Half of the nucleic acid from CSF was treated with Turbo DNase (Ambion, Waltham Mass.), followed by reverse transcription of the RNA to cDNA using random hexamers and NGS library preparation using the Nextera XT DNA Library Prep Kit (Illumina, San Diego Calif.). The remaining half was treated with the NEBNext Microbiome DNA Enrichment Kit to enrich for microbial DNA (New England Biosciences, Ipswich Mass.)), followed by Nextera XT library preparation. Dual-indexed, barcoded NGS libraries were quantitated on the BioAnalyzer (Agilent, Santa Clara Calif.) and run on the Illumina HiSeq (1×160 bp run).

The patient's CSF sample was processed and sequenced as part of a recently developed standardized operating procedure (SOP) for clinical mNGS testing from patient samples in the University of California, San Francisco (UCSF) Clinical Microbiology Laboratory, a CLIA (Clinical Laboratory Improvement Amendments)-licensed laboratory (Naccache, et al., manuscript in preparation). For each sequencing run, the SOP includes running two external controls: (1) a negative “no-template” control (NTC) sample consisting of elution buffer, (2) and a positive control sample consisting of a quantified mixture of 7 representative pathogens (CMV, HIV, Streptococcus agalactiae, Klebsiella pneumoniae, Cryptococcus neoformans, Aspergillus niger, and Toxoplasma gondii). Each is spiked into negative CSF matrix at a concentration 1-2 log above the estimated limits of detection for that microorganism by probit analysis (Naccache, et al., manuscript in preparation).

2. Bioinformatics Analysis

Metagenomic NGS data was analyzed for pathogens using a modified version of the SURPI (“sequence-based ultra-rapid pathogen identification”) computational pipeline, which identifies pathogen sequences on the basis of nucleotide alignments to National Center for Biotechnology information (NCBI) nt reference database (March 2015 build) [4]. A clinical version of the SURPI pipeline, named SURPI+, which employs taxonomic classification for more accurate read assignments and establish normalized metrics and thresholds for clinical results reporting, was used for automated interpretation. Briefly, for identification of viral sequences, both mNGS data from bath DNA and RNA sample libraries were used, whereas only mNGS data from DNA libraries alone were used for identification of sequences from bacteria, fungi, and parasites. An edit distance cutoff of 16, indicating the number of single nucleotide insertions, deletions, or mismatches allowed between the read and the reference sequence, was used for virus detection, whereas a more stringent edit distance cutoff of 1 was used far bacterial, fungal, and parasitic detection [2, 3]. Following alignment, a rapid taxonomic classification algorithm based on the lowest common ancestor algorithm was used to assign viral, bacterial, and non-chordate eukaryotic (fungal or parasitic) NGS reads to the species, genus, or family level, as previously described [2, 3].

3. Results Reporting

As part of clinical validation of mNGS testing for pathogen detection, we have established threshold cutoffs for automated reporting of positive results (Naccache, et al., manuscript in preparation). Briefly, for reporting of bacteria, fungi, and parasites, the cutoff is defined as an RPM (reads per million) ratio of ≥10, where the RPM ratio is defined as the RPM_(sample)/RPM_(NTC) for any given taxon (species, genus, or family). If the taxon is not present in the NTC, then the RPM_(NTC) is 1. For reporting of viruses, the criteria can include coverage of ≥2 non-contiguous/non-overlapping gene regions. Viruses with non-vertebrate hosts, that are found in the NTC, or that constitute normal body flora (e.g. anelloviruses) are not reported. A “gene region” can be defined as a non-overlapping and non-contiguous region of length<=read length (e.g., 141 bp).

4. Brucella PCR Confirmation

Confirmation of Brucella detection was performed by PCR using the following primer set targeting the IS711 gene [5]: BrucellaGenus_F-2_GCCTTGGATCTGAGCCGTT (SEQ ID NO:1); BrucellaGenus_IS711_R GGCCTACCGCTGCGAAT (SEQ ID NO:2). The reaction was carried out using the Qiagen One-Step RT-PCR Kit in a 25 μL total reaction volume by addition of 4 μL Q solution, 4 μL. 5× buffer, 1 μL dNTP, 1 μL enzyme, 1 μL of each primer at 12 pmol, and 2 μL of template extracted DNA, and water. Cycling conditions were as follows: 40 cycles of 94° C., 30 s/53° C., 30 s/72° C., 30. The recovered sequence was confirmed to be Bruce by Sanger sequencing of the PCR amplicon.

VIII. Exemplary Computer System

FIG. 14 shows a block diagram of an example computer system 10 usable with system and methods according to embodiments of the present invention. For example, computer system 10 can be used in implementing sample analysis system 110 and/or sequence identification computer 120.

Any of the computer systems mentioned herein may utilize any suitable number of subsystems. Examples of such subsystems are shown in FIG. 14 in computer apparatus 10. In some embodiments, a computer system includes a single computer apparatus, where the subsystems can be the components of the computer apparatus. In other embodiments, a computer system can include multiple computer apparatuses, each being a subsystem, with internal components. A computer system can include desktop and laptop computers, tablets, mobile phones and other mobile devices.

The subsystems shown in FIG. 14 are interconnected via a system bus 75. Additional subsystems such as a printer 74, keyboard 78, storage device(s) 79, monitor 76, which is coupled to display adapter 82, and others are shown. Peripherals and input/output (I/O) devices, which couple to I/O controller 71, can be connected to the computer system by any number of means known in the art such as input/output (I/O) port 77 (e.g., USB, FireWire®). For example, I/O port 77 or external interface 81 (e.g. Ethernet, Wi-Fi, etc.)can be used to connect computer system 10 to a wide area network such as the Internet, a mouse input device, or a scanner. The interconnection via system bus 75 allows the central processor 73 to communicate with each subsystem and to control the execution of a plurality of instructions from system memory 72 or the storage device(s) 79 (e.g., a fixed disk, such as a hard drive, or optical disk), as well as the exchange of information between subsystems. The system memory 72 and/or the storage device(s) 79 may embody a computer readable medium. Another subsystem is a data collection device 85, such as a camera, microphone, accelerometer, and the like. Any of the data mentioned herein can be output from one component to another component and can be output to the user.

A computer system can include a plurality of the same components or subsystems, e.g., connected together by external interface 81, by an internal interface, or via removable storage devices that can be connected and removed from one component to another component. In some embodiments, computer systems, subsystem, or apparatuses can communicate over a network. In such instances, one computer can be considered a client and another computer a server, where each can be part of a same computer system. A client and a server can each include multiple systems, subsystems, or components.

Aspects of embodiments can be implemented in the form of control logic using hardware (e.g. an application specific integrated circuit or field programmable gate array) and/or using computer software with a generally programmable processor in a modular or integrated manner. As used herein, a processor includes a single-core processor, multi-core processor on a same integrated chip, or multiple processing units on a single circuit board or networked. Based on the disclosure and teachings provided herein, a person of ordinary skill in the art will know and appreciate other ways and/or methods to implement embodiments of the present invention using hardware and a combination of hardware and software.

Any of the software components or functions described in this application may be implemented as software code to be executed by a processor using any suitable computer language such as, for example, Java, C, C++. C#, Objective-C, Go, Swift, or scripting language such as Perl or Python using, for example, conventional or object-oriented techniques. The software code may be stored as a series of instructions or commands on a computer readable medium for storage and/or transmission. A suitable non-transitory computer readable medium can include random access memory (RAM), a read only memory (ROM), a magnetic medium such as a hard-drive or a floppy disk, or an optical medium such as a compact disk (CD) or DVD (digital versatile disk), flash memory, and the like. The computer readable medium may be any combination of such storage or transmission devices.

Such programs may also be encoded and transmitted using carrier signals adapted for transmission via wired, optical, and/or wireless networks conforming to a variety of protocols, including the Internet. As such, a computer readable medium may be created using a data signal encoded with such programs. Computer readable media encoded with the program code may be packaged with a compatible device or provided separately from other devices (e.g., via Internet download). Any such computer readable medium may reside on or within a single computer product (e.g. a hard drive, a CD, or an entire computer system), and may be present on or within different computer products within a system or network. A computer system may include a monitor, printer, or other suitable display for providing any of the results mentioned herein to a user.

Any of the methods described herein may be totally or partially performed with a computer system including one or more processors, which can be configured to perform the steps. Thus, embodiments can be directed to computer systems configured to perform the steps of any of the methods described herein, potentially with different components performing a respective steps or a respective group of steps. Although presented as numbered steps, steps of methods herein can be performed at a same time or in a different order. Additionally, portions of these steps may be used with portions of other steps from other methods. Also, all or portions of a step may be optional. Additionally, any of the steps of any of the methods can be performed with modules, units, circuits, or other means for performing these steps.

The specific details of particular embodiments may be combined in any suitable manner without departing from the spirit and scope of embodiments of the invention. However, other embodiments of the invention may be directed to specific embodiments relating to each individual aspect, or specific combinations of these individual aspects.

The above description of example embodiments of the invention has been presented for the purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise form described, and many modifications and variations are possible in light of the teaching above.

A recitation of “a”, “an” or “the” is intended to mean “one or more” unless specifically indicated to the contrary. The use of “or” is intended to mean an “inclusive or,” and not an “exclusive or” unless specifically indicated to the contrary. Reference to a “first” component does not necessarily require that a second component be provided. Moreover reference to a “first” or a “second” component does not limit the referenced component to a particular location unless expressly stated.

All patents, patent applications, publications, and descriptions mentioned herein are herein incorporated by reference in their entirety for all purposes. None is admitted to be prior art. Such references include:

-   1. Pappas G, Akritidis N, Bosilkovski M, Tsianos E. Brucellosis. The     New England journal of medicine 2005; 352(22): 2325-36. -   2. Greninger A L, Messacar K, Dunnebacke T, et al. Clinical     metagenomic identification of Balamuthia mandrillaris encephalitis     and assembly of the draft genome: the continuing case for reference     genome sequencing. Genome medicine 2015; 7(1): 113. -   3. Greninger A L, Naccache S N, Messacar K, et al. A novel outbreak     enterovirus D68 strain associated with acute flaccid myelitis cases     in the USA (2012-14): a retrospective cohort study. The Lancet     Infectious diseases 2015; 15(6): 671-82. -   4. Naccache S N, Federman S, Veeraraghavan N, et al. A     cloud-compatible bioinformatics pipeline for ultrarapid pathogen     identification from next-generation sequencing of clinical samples.     Genome research 2014; 24(7): 1180-92. -   5. Hinic V, Brodard I, Thomann A, et al. Novel identification and     differentiation of Brucella melitensis, B. abortus, B. suis, B.     ovis, B. canis, and B. neotomae suitable for both conventional and     real-time PCR systems. Journal of microbiological methods 2008;     75(2): 375 -   6. Al-Sous M W, Bolilega S, Al-Kawi M Z, Alwatban J, McLean D R.     Neurobrucellosis: clinical and neuroimaging correlation. AJNR     American journal of neuroradiology 2004; 25(3): 395-401. -   7. Moyer N P H L, Murray P R, Baron E J, Pfaller M A, Tenover F C,     Youlken R H. Brucella. Manual of Clinical Microbiology. ASM Press,     1995:549-55. -   8. Yu W L, Nielsen K. Review of detection of Brucella spp. by     polymerase chain reaction. Croatian medical journal 2010; 51(4):     306-13. -   9. Chiu C Y, Miller S. Next-Generation Sequencing. In: Persing D H,     Tenover F C, Hayden R T, Ieven G, Miller M B, Nolte F S. Molecular     Microbiology, Diagnostic Principles and Practice, 3rd Edition.     Washington, D.C.: ASM Press, 2016:68-79 -   10. Dasari S, Naha K, Prabhu M. Brucellosis and tuberculosis:     clinical overlap and pitfalls. Asian Pacific journal of tropical     medicine 2013; 6(10): 823-5. -   11. Verdon R, Chevret S, Laissy J P, Wolff M. Tuberculous meningitis     in adults: review of 48 cases. Clinical infectious diseases: an     official publication of the Infectious Diseases Society of America     1996; 22(6): 982-8. -   12. Naccache S N, Peggs K S, Mattes F M, et al. Diagnosis of     neuroinvasive astrovirus infection in an immunocompromised adult     with encephalitis by unbiased next-generation sequencing. Clinical     infectious diseases: an official publication of the Infectious     Diseases Society of America 2015; 60(6): 919-23. -   13. Wilson M R, Naccache S N, Samayoa E, et al. Actionable diagnosis     of neuroleptospirosis by next-generation sequencing. The New England     journal of medicine 2014; 370(25): 2408-17.

APPENDIX

SUPPLEMENTARY TABLE 1 # of Reads DNA DNA DNA Patient species genus family NTC PC CSF * Brucella Brucellaceae 0 0 251 Brucella melitensis Brucella Brucellaceae 0 0 26 Streptococcus agalactiae Streptococcus Streptococcaceae 9 50,894 0 Klebsiella pneumoniae Klebsiella Enterobacteriaceae 426 39,013 4 * * Enterobacteriaceae 246,369 21,342 2,350 Escherichia coli Escherichia Enterobacteriaceae 170,966 11,959 1,601 Propionibacterium acnes Propionibacterium Propionibacteriaceae 126,533 5,329 1,158 * Streptococcus Streptococcaceae 236 3,556 2 * Klebsiella Enterobacteriaceae 49 3,542 0 * * * 41,826 3,194 457 Streptococcus suis Streptococcus Streptococcaceae 20 561 0 Streptococcus salivarius Streptococcus Streptococcaceae 188 333 13 Streptococcus pyogenes Streptococcus Streptococcaceae 0 299 0 Serratia marcescens Serratia Enterobacteriaceae 59 233 0 Lactococcus lactis Lactococcus Streptococcaceae 338 223 2 * Escherichia Enterobacteriaceae 3,865 210 35 Pseudomonas sp. TKP Pseudomonas Pseudomonadaceae 941 205 69 * Pseudomonas Pseudomonadaceae 499 185 11 Streptococcus Streptococcus Streptococcaceae 0 177 0 macedonicus Enterobacter cloacae Enterobacter Enterobacteriaceae 235 160 1 Klebsiella oxytoca Klebsiella Enterobacteriaceae 102 129 4 Thermoanaerobacterium Thermoanaerobacterium Thermoanaerobacterales 268 107 2 thermosaccharolyticum Family III. Incertae Sedis Streptococcus sp. VT 162 Streptococcus Streptococcaceae 90 100 0 Enterococcus faecium Enterococcus Enterococcaceae 526 96 85 Pseudomonas protegens Pseudomonas Pseudomonadaceae 17 87 3 Micrococcus luteus Micrococcus Micrococcaceae 6,424 78 3 Streptococcus infantarius Streptococcus Streptococcaceae 94 75 0 Staphylococcus Staphylococcus Staphylococcaceae 8,168 74 13 epidermidis Enterobacter asburiae Enterobacter Enterobacteriaceae 250 65 0 Cupriavidus metallidurans Cupriavidus Burkholderiaceae 934 61 0 Pseudomonas putida Pseudomonas Pseudomonadaceae 300 57 1 Enterococcus Enterococcus Enterococcaceae 371 54 79 casseliflavus Streptococcus mitis Streptococcus Streptococcaceae 55 53 0 Streptococcus equi Streptococcus Streptococcaceae 0 51 0 Salmonella enterica Salmonella Enterobacteriaceae 1 48 0 Streptococcus oralis Streptococcus Streptococcaceae 51 47 0 Klebsiella variicola Klebsiella Enterobacteriaceae 0 42 0 Burkholderia lata Burkholderia Burkholderiaceae 85 40 3 Pseudomonas stutzeri Pseudomonas Pseudomonadaceae 1,094 40 32 Streptococcus Streptococcus Streptococcaceae 63 39 3 pneumoniae Streptococcus Streptococcus Streptococcaceae 3 36 0 dysgalactiae * Burkholderia Burkholderiaceae 171 35 18 Pseudomonas fluorescens Pseudomonas Pseudomonadaceae 50 35 4 Acinetobacter guillouiae Acinetobacter Moraxellaceae 207 35 4 Veillonella parvula Veillonella Veillonellaceae 151 33 7 Xanthomonas campestris Xanthomonas Xanthomonadaceae 164 32 24 Exiguobacterium sp. Exiguobacterium 271 32 0 AT1b Pseudomonas Pseudomonas Pseudomonadaceae 932 32 0 pseudoalcaligenes Streptococcus Streptococcus Streptococcaceae 0 27 0 pasteurianus Rothia dentocariosa Rothia Micrococcaceae 270 25 2 * * Streptococcaceae 0 25 0 Delftia acidovorans Delftia Comamonadaceae 28 25 0 * Propionibacterium Propionibacteriaceae 773 23 8 Acinetobacter baumannii Acinetobacter Moraxellaceae 853 23 5 Streptococcus Streptococcus Streptococcaceae 393 23 2 parasanguinis Pseudomonas sp. Pseudomonas Pseudomonadaceae 15 22 1 WCS374 Enterobacter aerogenes Enterobacter Enterobacteriaceae 0 22 0 Stenotrophomonas Stenotrophomonas Xanthomonadaceae 78 22 0 maltophilia Alicyclobacillus Alicyclobacillus Alicyclobacillaceae 79 21 2 acidocaldarius Haemophilus influenzae Haemophilus Pasteurellaceae 145 19 13 Rothia mucilaginosa Rothia Micrococcaceae 379 19 3 Staphylococcus xylosus Staphylococcus Staphylococcaceae 68 18 3 Acidovorax sp. JS42 Acidovorax Comamonadaceae 92 18 3 * Rahnella Enterobacteriaceae 0 18 0 Streptococcus Streptococcus Streptococcaceae 30 18 0 pseudopneumoniae Streptococcus Streptococcus Streptococcaceae 114 17 4 thermophilus Pseudomonas aeruginosa Pseudomonas Pseudomonadaceae 114 17 1 Corynebacterium Corynebacterium Corynebacteriaceae 0 16 0 kroppenstedtii Staphylococcus Staphylococcus Staphylococcaceae 1,320 16 0 haemolyticus Serratia sp. SCBI Serratia Enterobacteriaceae 0 15 0 * Enterococcus Enterococcaceae 123 14 5 Staphylococcus warneri Staphylococcus Staphylococcaceae 326 14 3 Burkholderia cepacia Burkholderia Burkholderiaceae 63 13 9 * Corynebacterium Corynebacteriaceae 167 13 7 Comamonas testosteroni Comamonas Comamonadaceae 175 13 1 Bifidobacterium Bifidobacterium Bifidobacteriaceae 0 13 0 thermophilum * Lactobacillus Lactobacillaceae 50 13 0 Thermoanaerobacterium Thermoanaerobacterium Thermoanaerobacterales 56 13 0 xylanolyticum Family III. Incertae Sedis Ralstonia pickettii Ralstonia Burkholderiaceae 147 13 0 Meiothermus ruber Meiothermus Thermaceae 540 13 0 Acidovorax ebreus Acidovorax Comamonadaceae 10 12 0 Micrococcus sp. V7 Micrococcus Micrococcaceae 36 12 0 Leuconostoc Leuconostoc Leuconostocaceae 335 11 1 mesenteroides Bacillus halodurans Bacillus Bacillaceae 9 11 0 Corynebacterium Corynebacterium Corynebacteriaceae 20 11 0 variabile * Acidovorax Comamonadaceae 48 11 0 Bifidobacterium bifidum Bifidobacterium Bifidobacteriaceae 619 11 0 Rhizobium sp. IRBG74 Rhizobium Rhizobiaceae 182 10 4 * Acinetobacter Moraxellaceae 569 10 1 Acinetobacter Acinetobacter Moraxellaceae 0 10 0 calcoaceticus Chroococcidiopsis Chroococcidiopsis 28 10 0 thermalis Streptococcus sanguinis Streptococcus Streptococcaceae 47 10 0 Acidovorax sp. KKS102 Acidovorax Comamonadaceae 74 10 0 Raoultella ornithinolytica Raoultella Enterobacteriaceae 94 10 0 Ochrobactrum anthropi Ochrobactrum Brucellaceae 270 10 0 Lactobacillus johnsonii Lactobacillus Lactobacillaceae 9 9 2 Methylobacterium populi Methylobacterium Methylobacteriaceae 30 9 3 Rhodococcus equi Rhodococcus Nocardiaceae 62 9 4 Lactobacillus helveticus Lactobacillus Lactobacillaceae 4 9 0 * Serratia Enterobacteriaceae 30 9 0 Burkholderia cenocepacia Burkholderia Burkholderiaceae 15 8 1 * Staphylococcus Staphylococcaceae 1,001 8 6 Enterobacter sp. R4-368 Enterobacter Enterobacteriaceae 0 8 0 Propionibacterium Propionibacterium Propionibacteriaceae 31 8 0 propionicum Streptococcus gordonii Streptococcus Streptococcaceae 35 8 0 Corynebacterium Corynebacterium Corynebacteriaceae 48 8 0 singulare Burkholderia ambifaria Burkholderia Burkholderiaceae 0 7 1 * Micrococcus Micrococcaceae 142 7 2 Fervidobacterium Fervidobacterium Thermotogaceae 83 7 1 nodosum Aeromonas media Aeromonas Aeromonadaceae 433 7 3 Cronobacter sakazakii Cronobacter Enterobacteriaceae 0 7 0 Myroides profundi Myroides Flavobacteriaceae 0 7 0 Methylobacterium oryzae Methylobacterium Methylobacteriaceae 0 7 0 * Xanthomonas Xanthomonadaceae 13 7 0 Thermoanaerobacterium Thermoanaerobacterium Thermoanaerobacterales 73 7 0 saccharolyticum Family III. Incertae Sedis Pseudomonas mendocina Pseudomonas Pseudomonadaceae 84 7 0 Corynebacterium Corynebacterium Corynebacteriaceae 315 7 0 ureicelerivorans Lactobacillus crispatus Lactobacillus Lactobacillaceae 74 6 3 Alicycliphilus Alicycliphilus Comamonadaceae 132 6 2 denitrificans Gardnerella vaginalis Gardnerella Bifidobacteriaceae 1,140 6 2 * Gemella 0 6 0 * Ralstonia Burkholderiaceae 0 6 0 Eggerthella lenta Eggerthella Coriobacteriaceae 3 6 0 * * Rhizobiaceae 7 6 0 Prevotella denticola Prevotella Prevotellaceae 16 6 0 Prevotella intermedia Prevotella Prevotellaceae 18 6 0 Psychrobacter sp. PRwf-1 Psychrobacter Moraxellaceae 31 6 0 Azospira oryzae Azospira Rhodocyclaceae 38 6 0 Acinetobacter Acinetobacter Moraxellaceae 47 6 0 haemolyticus * Delftia Comamonadaceae 54 6 0 Burkholderia contaminans Burkholderia Burkholderiaceae 7 5 4 Arthrobacter arilaitensis Arthrobacter Micrococcaceae 20 5 7 Dermacoccus Dermacoccus Dermacoccaceae 31 5 6 nishinomiyaensis Pantoea ananatis Pantoea Enterobacteriaceae 40 5 5 Staphylococcus Staphylococcus Staphylococcaceae 148 5 8 saprophyticus Staphylococcus pasteuri Staphylococcus Staphylococcaceae 69 5 1 Rahnella aquatilis Rahnella Enterobacteriaceae 0 5 0 Rahnella sp. Y9602 Rahnella Enterobacteriaceae 0 5 0 Campylobacter concisus Campylobacter Campylobacteraceae 0 5 0 Geobacillus sp. WCH70 Geobacillus Bacillaceae 0 5 0 * Frankia Frankiaceae 0 5 0 Lactobacillus casei Lactobacillus Lactobacillaceae 0 5 0 Thiomonas intermedia Thiomonas 0 5 0 Streptococcus gallolyticus Streptococcus Streptococcaceae 5 5 0 Thioalkalivibrio Thioalkalivibrio Ectothiorhodospiraceae 10 5 0 sulfidiphilus * Bradyrhizobium Bradyrhizobiaceae 13 5 0 Bifidobacterium longum Bifidobacterium Bifidobacteriaceae 27 5 0 Corynebacterium falsenii Corynebacterium Corynebacteriaceae 32 5 0 Delftia sp. Cs1-4 Delftia Comamonadaceae 40 5 0 Acinetobacter sp. M131 Acinetobacter Moraxellaceae 79 5 0 Prevotella Prevotella Prevotellaceae 151 5 0 melaninogenica * * Comamonadaceae 217 5 0 Leuconostoc carnosum Leuconostoc Leuconostocaceae 13 4 1 Pectobacterium Pectobacterium Enterobacteriaceae 74 4 3 carotovorum * Myroides Flavobacteriaceae 0 4 0 * Erwinia Enterobacteriaceae 0 4 0 Gordonia sp. KTR9 Gordonia Gordoniaceae 0 4 0 Paenibacillus sp. FSL R7- Paenibacillus Paenibacillaceae 0 4 0 0273 Paracoccus sp. N81106 Paracoccus Rhodobacteraceae 0 4 0 Sphingobium fuliginis Sphingobium Sphingomonadaceae 0 4 0 * * 1 4 0 * Geobacillus Bacillaceae 10 4 0 Pseudomonas sp. Pseudomonas Pseudomonadaceae 11 4 0 VLB120 Pelagibacterium Pelagibacterium Hyphomicrobiaceae 22 4 0 halotolerans Streptococcus intermedius Streptococcus Streptococcaceae 34 4 0 Propionibacterium Propionibacterium Propionibacteriaceae 45 4 0 freudenreichii * Enterobacter Enterobacteriaceae 72 4 0 Nakamurella multipartita Nakamurella Nakamurellaceae 32 3 7 Haemophilus parasuis Haemophilus Pasteurellaceae 0 3 0 Fusobacterium nucleatum Fusobacterium Fusobacteriaceae 0 3 0 Citrobacter freundii Citrobacter Enterobacteriaceae 1 3 0 Ruminococcus sp. SR1/5 Ruminococcus Ruminococcaceae 9 3 0 Pseudoxanthomonas Pseudoxanthomonas Xanthomonadaceae 10 3 0 spadix Lactococcus garvieae Lactococcus Streptococcaceae 15 3 0 Neisseria elongata Neisseria Neisseriaceae 18 3 0 Acidovorax citrulli Acidovorax Comamonadaceae 19 3 0 Novosphingobium Novosphingobium Sphingomonadaceae 20 3 0 pentaromativorans Citrobacter koseri Citrobacter Enterobacteriaceae 26 3 0 Methylobacterium Methylobacterium Methylobacteriaceae 32 3 0 aquaticum Pseudomonas Pseudomonas Pseudomonadaceae 34 3 0 denitrificans Rhodococcus erythropolis Rhodococcus Nocardiaceae 36 3 0 Lactobacillus reuteri Lactobacillus Lactobacillaceae 48 3 0 Bacteroides fragilis Bacteroides Bacteroidaceae 70 3 0 Lactobacillus plantarum Lactobacillus Lactobacillaceae 137 3 0 * Bacillus Bacillaceae 196 3 0 Pseudomonas Pseudomonas Pseudomonadaceae 0 2 1 rhizosphaerae Achromobacter Achromobacter Alcaligenaceae 41 2 1 xylosoxidans Lactobacillus amylovorus Lactobacillus Lactobacillaceae 0 2 0 Propionibacterium Propionibacterium Propionibacteriaceae 0 2 0 acidipropionici Leuconostoc gelidum Leuconostoc Leuconostocaceae 0 2 0 Weissella thailandensis Weissella Leuconostocaceae 0 2 0 Pandoraea sp. RB-44 Pandoraea Burkholderiaceae 0 2 0 Escherichia vulneris Escherichia Enterobacteriaceae 0 2 0 Yersinia intermedia Yersinia Enterobacteriaceae 0 2 0 Flavobacteriaceae Flavobacteriaceae 0 2 0 bacterium 3519-10 * Rhodococcus Nocardiaceae 0 2 0 Streptococcus sp. (N1) Streptococcus Streptococcaceae 0 2 0 * Methylobacterium Methylobacteriaceae 8 2 0 Sphingomonas sp. MM-1 Sphingomonas Sphingomonadaceae 9 2 0 Rhizobium etli Rhizobium Rhizobiaceae 18 2 0 * Agrobacterium Rhizobiaceae 41 2 0 Thermus scotoductus Thermus Thermaceae 53 2 0 Methylobacterium Methylobacterium Methylobacteriaceae 60 2 0 extorquens Streptococcus sp. I-P16 Streptococcus Streptococcaceae 63 2 0 Pantoea sp. PSNIH1 Pantoea Enterobacteriaceae 76 2 0 Methylobacterium Methylobacterium Methylobacteriaceae 96 2 0 radiotolerans [ Ruminococcus ] torques Blautia Lachnospiraceae 0 1 7 Lactobacillus sakei Lactobacillus Lactobacillaceae 0 1 2 * * Rhodocyclaceae 0 1 1 Bacillus licheniformis Bacillus Bacillaceae 34 1 5 Corynebacterium Corynebacterium Corynebacteriaceae 0 1 0 accolens Corynebacterium sp. Corynebacterium Corynebacteriaceae 0 1 0 ATCC 6931 Serratia symbiotica Serratia Enterobacteriaceae 0 1 0 Lactobacillus delbrueckii Lactobacillus Lactobacillaceae 0 1 0 Bacillus sp. YP1 Bacillus Bacillaceae 0 1 0 Klebsiella sp. PG122E Klebsiella Enterobacteriaceae 0 1 0 * Rathayibacter Microbacteriaceae 0 1 0 Pseudoalteromonas sp. Pseudoalteromonas Pseudoalteromonadaceae 0 1 0 P30 Staphylococcus sp. Staphylococcus Staphylococcaceae 0 1 0 CDC25 Corynebacterium Corynebacterium Corynebacteriaceae 1 1 0 resistens Shigella dysenteriae Shigella Enterobacteriaceae 9 1 0 * * Xanthomonadaceae 11 1 0 Agrobacterium fabrum Agrobacterium Rhizobiaceae 19 1 0 Gordonia Gordonia Gordoniaceae 22 1 0 polyisoprenivorans Pseudomonas balearica Pseudomonas Pseudomonadaceae 25 1 0 * * Pseudomonadaceae 25 1 0 Ruminococcus bromii Ruminococcus Ruminococcaceae 29 1 0 Brachybacterium faecium Brachybacterium Dermabacteraceae 31 1 0 Acinetobacter johnsonii Acinetobacter Moraxellaceae 37 1 0 Micrococcus sp. A1 Micrococcus Micrococcaceae 95 1 0 Filifactor alocis Filifactor Peptostreptococcaceae 132 1 0 Pantoea vagans Pantoea Enterobacteriaceae 162 1 0 Haemophilus Haemophilus Pasteurellaceae 241 1 0 parainfluenzae Pantoea rwandensis Pantoea Enterobacteriaceae 0 0 5 Corynebacterium Corynebacterium Corynebacteriaceae 0 0 4 vitaeruminis Pseudomonas poae Pseudomonas Pseudomonadaceae 8 0 3 * * Brucellaceae 0 0 2 Lactobacillus fermentum Lactobacillus Lactobacillaceae 0 0 2 Anabaena variabilis Anabaena Nostocaceae 0 0 2 Sphingobacterium sp. Sphingobacterium Sphingobacteriaceae 0 0 2 ML3W sugarcane isolate 74-1 0 0 2 * * Geodermatophilaceae 1 0 2 Megasphaera elsdenii Megasphaera Veillonellaceae 0 0 1 Pseudoxanthomonas Pseudoxanthomonas Xanthomonadaceae 7 0 1 suwonensis Corynebacterium Corynebacterium Corynebacteriaceae 157 0 17 glutamicum Sphingomonas taxi Sphingomonas Sphingomonadaceae 19 0 2 Pseudomonas graminis Pseudomonas Pseudomonadaceae 12 0 1 Bradyrhizobium sp. Bradyrhizobium Bradyrhizobiaceae 21 0 1 BTAi1 Enterococcus hirae Enterococcus Enterococcaceae 25 0 1 Corynebacterium sp. L2- Corynebacterium Corynebacteriaceae 34 0 1 79-05 Arthrobacter Arthrobacter Micrococcaceae 36 0 1 phenanthrenivorans Corynebacterium maris Corynebacterium Corynebacteriaceae 87 0 1 Gordonia bronchialis Gordonia Gordoniaceae 90 0 1 Kytococcus sedentarius Kytococcus Dermacoccaceae 130 0 1 Kosakonia cowanii Kosakonia Enterobacteriaceae 1 0 0 Xenorhabdus bovienii Xenorhabdus Enterobacteriaceae 1 0 0 Paracoccus haeundaensis Paracoccus Rhodobacteraceae 1 0 0 Methylobacterium sp. 238 Methylobacterium Methylobacteriaceae 1 0 0 Acinetobacter sp. BW3 Acinetobacter Moraxellaceae 1 0 0 Aeromonas sobria Aeromonas Aeromonadaceae 1 0 0 Bacillus lehensis Bacillus Bacillaceae 1 0 0 Ralstonia solanacearum Ralstonia Burkholderiaceae 1 0 0 Citrobacter sp. FPO3 Citrobacter Enterobacteriaceae 1 0 0 Citrobacter sp. I91-3 Citrobacter Enterobacteriaceae 1 0 0 Erwinia amylovora Erwinia Enterobacteriaceae 1 0 0 Klebsiella milletis Klebsiella Enterobacteriaceae 1 0 0 Salmonella bongori Salmonella Enterobacteriaceae 1 0 0 Serratia grimesii Serratia Enterobacteriaceae 1 0 0 Yersinia pestis Yersinia Enterobacteriaceae 1 0 0 * Enterobacteriaceae 1 0 0 Lactobacillus brevis Lactobacillus Lactobacillaceae 1 0 0 Kocuria sp. starX Kocuria Micrococcaceae 1 0 0 Acinetobacter sp. 26 Acinetobacter Moraxellaceae 1 0 0 Peptoclostridium difficile Peptoclostridium Peptostreptococcaceae 1 0 0 Sorangium cellulosum Sorangium Polyangiaceae 1 0 0 Pseudomonas sp. NSi14 Pseudomonas Pseudomonadaceae 1 0 0 Synergistetes oral clone 1 0 0 03 5 D05 bacterium EBAD26 1 0 0 bacterium NLAE-zl-G351 1 0 0 rumen bacterium 1 0 0 enrichment culture clone Y74 unidentified marine 2 0 0 bacterioplankton Escherichia albertii Escherichia Enterobacteriaceae 2 0 0 Mobiluncus curtisii Mobiluncus Actinomycetaceae 2 0 0 * Caulobacter Caulobacteraceae 2 0 0 Methylotenera versatilis Methylotenera Methylophilaceae 2 0 0 Propionibacterium sp. Propionibacterium Propionibacteriaceae 3 0 0 KPL1849 bacterium EBAD25 3 0 0 Bacillus sp. Pc3 Bacillus Bacillaceae 3 0 0 Acinetobacter sp. EVA14 Acinetobacter Moraxellaceae 3 0 0 Agrobacterium sp. Agrobacterium Rhizobiaceae 3 0 0 Alistipes shahii Alistipes Rikenellaceae 3 0 0 * Thermus Thermaceae 3 0 0 * Methylibium 3 0 0 butyrate-producing 3 0 0 bacterium SSC/2 Escherichia fergusonii Escherichia Enterobacteriaceae 4 0 0 Enterobacter sp. Ni15 Enterobacter Enterobacteriaceae 4 0 0 Capnocytophaga ochracea Capnocytophaga Flavobacteriaceae 4 0 0 Thauera sp. 6NLG Thauera Rhodocyclaceae 4 0 0 Desulfovibrio alaskensis Desulfovibrio Desulfovibrionaceae 4 0 0 Variovorax sp. Alb14 Variovorax Comamonadaceae 4 0 0 * Shigella Enterobacteriaceae 4 0 0 * * Micromonosporaceae 4 0 0 * Micromonospora Micromonosporaceae 4 0 0 Thermobifida fusca Thermobifida Nocardiopsaceae 4 0 0 Turneriella parva Turneriella Leptospiraceae 5 0 0 [Clostridium] sticklandii Peptoclostridium Peptostreptococcaceae 5 0 0 Acinetobacter sp. Ooi24 Acinetobacter Moraxellaceae 5 0 0 Ochrobactrum sp. SJY1 Ochrobactrum Brucellaceae 5 0 0 Carnobacterium sp. Carnobacterium Carnobacteriaceae 5 0 0 WN1359 Iamia majanohamensis Iamia Iamiaceae 5 0 0 Saccharomonospora Saccharomonospora Pseudonocardiaceae 5 0 0 viridis Rhizobium sp. Rhizobium Rhizobiaceae 5 0 0 Staphylococcus sp. CDC3 Staphylococcus Staphylococcaceae 5 0 0 Shigella sonnei Shigella Enterobacteriaceae 6 0 0 Pseudomonas syringae Pseudomonas Pseudomonadaceae 6 0 0 Burkholderia Burkholderia Burkholderiaceae 6 0 0 vietnamiensis Shigella boydii Shigella Enterobacteriaceae 6 0 0 Bacillus Bacillus Bacillaceae 6 0 0 weihenstephanensis Erythrobacter litoralis Erythrobacter Erythrobacteraceae 6 0 0 Pseudoalteromonas Pseudoalteromonas Pseudoalteromonadaceae 6 0 0 haloplanktis Pseudomonas sp. FGI182 Pseudomonas Pseudomonadaceae 6 0 0 * Rhizobium Rhizobiaceae 6 0 0 * Rickettsia Rickettsiaceae 6 0 0 Sphingobium yanoikuyae Sphingobium Sphingomonadaceae 6 0 0 Stenotrophomonas Stenotrophomonas Xanthomonadaceae 6 0 0 rhizophila * Leuconostoc Leuconostocaceae 7 0 0 Aquincola tertiaricarbonis Aquincola 7 0 0 Nocardiopsis dassonvillei Nocardiopsis Nocardiopsaceae 7 0 0 Carnobacterium Carnobacterium Carnobacteriaceae 7 0 0 maltaromaticum * Haemophilus Pasteurellaceae 7 0 0 Bordetella parapertussis Bordetella Alcaligenaceae 7 0 0 * Dietzia Dietziaceae 7 0 0 Shewanella sp. W3-18-1 Shewanella Shewanellaceae 7 0 0 Sphingomonas sp. NP5 Sphingomonas Sphingomonadaceae 7 0 0 Staphylococcus Staphylococcus Staphylococcaceae 7 0 0 gallinarum Micavibrio Micavibrio 7 0 0 aeruginosavorus Paracoccus denitrificans Paracoccus Rhodobacteraceae 8 0 0 [ Cellvibrio ] gilvus Cellulomonas Cellulomonadaceae 8 0 0 Corynebacterium Corynebacterium Corynebacteriaceae 8 0 0 jeikeium * * Staphylococcaceae 8 0 0 Meiothermus silvanus Meiothermus Thermaceae 8 0 0 Asticcacaulis excentricus Asticcacaulis Caulobacteraceae 8 0 0 * Atopobium Coriobacteriaceae 8 0 0 Streptococcus constellatus Streptococcus Streptococcaceae 8 0 0 Microcystis aeruginosa Microcystis 8 0 0 agricultural soil bacterium 8 0 0 SC-I-13 [ Ruminococcus ] obeum Blautia Lachnospiraceae 9 0 0 Thermus thermophilus Thermus Thermaceae 9 0 0 Shigella flexneri Shigella Enterobacteriaceae 9 0 0 * Mycobacterium Mycobacteriaceae 9 0 0 Pseudomonas savastanoi Pseudomonas Pseudomonadaceae 9 0 0 Staphylococcus capitis Staphylococcus Staphylococcaceae 9 0 0 * Cupriavidus Burkholderiaceae 9 0 0 Dyadobacter fermentans Dyadobacter Cytophagaceae 9 0 0 Dietzia sp. CQ4 Dietzia Dietziaceae 9 0 0 * * Methylophilaceae 9 0 0 * Neisseria Neisseriaceae 9 0 0 Corynebacterium Corynebacterium Corynebacteriaceae 10 0 0 aurimucosum * 10 0 0 Pseudomonas fulva Pseudomonas Pseudomonadaceae 10 0 0 Chromohalobacter Chromohalobacter Halomonadaceae 10 0 0 salexigens Brevundimonas diminuta Brevundimonas Caulobacteraceae 10 0 0 Streptococcus lutetiensis Streptococcus Streptococcaceae 10 0 0 Bordetella petrii Bordetella Alcaligenaceae 10 0 0 Erythrobacter sp. JP13.1 Erythrobacter Erythrobacteraceae 10 0 0 Methylobacillus Methylobacillus Methylophilaceae 10 0 0 glycogenes Candidatus Rhodoluna Candidatus Rhodoluna Microbacteriaceae 10 0 0 lacicola Arthrobacter sp. JBH1 Arthrobacter Micrococcaceae 10 0 0 Aggregatibacter Aggregatibacter Pasteurellaceae 10 0 0 aphrophilus Thauera sp. B4 Thauera Rhodocyclaceae 10 0 0 Lysobacter dokdonensis Lysobacter Xanthomonadaceae 10 0 0 Clostridiales genomosp. 10 0 0 BVAB3 Streptococcus sp. I-G2 Streptococcus Streptococcaceae 11 0 0 Pseudomonas mandelii Pseudomonas Pseudomonadaceae 11 0 0 Bradyrhizobium sp. Bradyrhizobium Bradyrhizobiaceae 11 0 0 S23321 Phenylobacterium Phenylobacterium Caulobacteraceae 11 0 0 zucineum Pseudomonas mosselii Pseudomonas Pseudomonadaceae 11 0 0 Staphylococcus Staphylococcus Staphylococcaceae 11 0 0 lugdunensis Proteus mirabilis Proteus Enterobacteriaceae 11 0 0 * * Neisseriaceae 11 0 0 Arthrobacter sp. J3-40 Arthrobacter Micrococcaceae 11 0 0 * Pantoea Enterobacteriaceae 12 0 0 Corynebacterium Corynebacterium Corynebacteriaceae 12 0 0 efficiens * Halomonas Halomonadaceae 12 0 0 Trueperella pyogenes Trueperella Actinomycetaceae 12 0 0 Streptomyces coelicolor Streptomyces Streptomycetaceae 12 0 0 Kocuria rhizophila Kocuria Micrococcaceae 13 0 0 Bacillus cereus Bacillus Bacillaceae 13 0 0 Tannerella forsythia Tannerella Porphyromonadaceae 13 0 0 * Alkalibacterium Carnobacteriaceae 13 0 0 Atopobium parvulum Atopobium Coriobacteriaceae 13 0 0 Serinicoccus profundi Serinicoccus Intrasporangiaceae 13 0 0 * Leptotrichia Leptotrichiaceae 13 0 0 * * Planococcaceae 13 0 0 Planomicrobium Planomicrobium Planococcaceae 13 0 0 okeanokoites Jannaschia sp. CCS1 Jannaschia Rhodobacteraceae 13 0 0 Paracoccus aestuarii Paracoccus Rhodobacteraceae 13 0 0 Rhodobacter blasticus Rhodobacter Rhodobacteraceae 13 0 0 Agrobacterium Agrobacterium Rhizobiaceae 14 0 0 tumefaciens Shewanella sp. ANA-3 Shewanella Shewanellaceae 14 0 0 Pseudomonas cichorii Pseudomonas Pseudomonadaceae 14 0 0 Halomonas sp. A3H3 Halomonas Halomonadaceae 14 0 0 Serratia liquefaciens Serratia Enterobacteriaceae 14 0 0 Sphingopyxis alaskensis Sphingopyxis Sphingomonadaceae 14 0 0 * Brevundimonas Caulobacteraceae 14 0 0 Deinococcus deserti Deinococcus Deinococcaceae 14 0 0 Desulfovibrio vulgaris Desulfovibrio Desulfovibrionaceae 14 0 0 Propionibacterium sp. Propionibacterium Propionibacteriaceae 14 0 0 NTS31307302 Paracoccus alcaliphilus Paracoccus Rhodobacteraceae 14 0 0 Vibrio parahaemolyticus Vibrio Vibrionaceae 14 0 0 Candidatus 14 0 0 Saccharibacteria oral taxon TM7x Geitlerinema sp. PCC Geitlerinema 15 0 0 7407 * Actinomyces Actinomycetaceae 15 0 0 Brevundimonas Brevundimonas Caulobacteraceae 15 0 0 vesicularis Acinetobacter sp. YS0810 Acinetobacter Moraxellaceae 15 0 0 * Prevotella Prevotellaceae 15 0 0 Methyloceanibacter Methyloceanibacter 15 0 0 caenitepidi Leuconostoc citreum Leuconostoc Leuconostocaceae 16 0 0 * Bacteroides Bacteroidaceae 16 0 0 Pseudomonas alcaligenes Pseudomonas Pseudomonadaceae 16 0 0 Methylibium Methylibium 16 0 0 petroleiphilum Moraxella catarrhalis Moraxella Moraxellaceae 16 0 0 Sphingopyxis sp. Kp5.2 Sphingopyxis Sphingomonadaceae 16 0 0 Pandoraea apista Pandoraea Burkholderiaceae 16 0 0 * Cellulomonas Cellulomonadaceae 16 0 0 Olsenella uli Olsenella Coriobacteriaceae 16 0 0 Acinetobacter oleivorans Acinetobacter Moraxellaceae 16 0 0 Sphingomonas sp. 133 Sphingomonas Sphingomonadaceae 16 0 0 Meiothermus taiwanensis Meiothermus Thermaceae 16 0 0 Deinococcus geothermalis Deinococcus Deinococcaceae 17 0 0 Lactobacillus Lactobacillus Lactobacillaceae 17 0 0 sanfranciscensis Acinetobacter sp. Acinetobacter Moraxellaceae 17 0 0 LUH5605 Staphylococcus simulans Staphylococcus Staphylococcaceae 17 0 0 Arsenophonus nasoniae Arsenophonus Enterobacteriaceae 17 0 0 Buchnera aphidicola Buchnera Enterobacteriaceae 17 0 0 Weissella koreensis Weissella Leuconostocaceae 17 0 0 Psychrobacter sp. G Psychrobacter Moraxellaceae 17 0 0 Amycolicicoccus Amycolicicoccus Mycobacteriaceae 17 0 0 subflavus Staphylococcus hyicus Staphylococcus Staphylococcaceae 17 0 0 Morganella morganii Morganella Enterobacteriaceae 18 0 0 Kyrpidia tusciae Kyrpidia Alicyclobacillaceae 18 0 0 Ramlibacter tataouinensis Ramlibacter Comamonadaceae 18 0 0 Weeksella virosa Weeksella Flavobacteriaceae 18 0 0 Acinetobacter junii Acinetobacter Moraxellaceae 18 0 0 Acinetobacter sp. 26B2 Acinetobacter Moraxellaceae 18 0 0 Mycobacterium abscessus Mycobacterium Mycobacteriaceae 18 0 0 Neisseria gonorrhoeae Neisseria Neisseriaceae 18 0 0 Sphingomonas wittichii Sphingomonas Sphingomonadaceae 19 0 0 Bacteroides dorei Bacteroides Bacteroidaceae 19 0 0 Corynebacterium Corynebacterium Corynebacteriaceae 19 0 0 halotolerans Pelomonas aquatica Pelomonas Comamonadaceae 19 0 0 Janibacter sp. TYM3221 Janibacter Intrasporangiaceae 19 0 0 * Arthrobacter Micrococcaceae 19 0 0 Mycobacterium gordonae Mycobacterium Mycobacteriaceae 19 0 0 Pimelobacter simplex Pimelobacter Nocardioidaceae 19 0 0 Pseudomonas sp. Pseudomonas Pseudomonadaceae 19 0 0 OM2164 Streptomyces albus Streptomyces Streptomycetaceae 20 0 0 Halomonas halocynthiae Halomonas Halomonadaceae 20 0 0 Nitrosomonas sp. AL212 Nitrosomonas Nitrosomonadaceae 20 0 0 Sphingobacterium Sphingobacterium Sphingobacteriaceae 20 0 0 mizutaii Vibrio cholerae Vibrio Vibrionaceae 20 0 0 Rickettsia felis Rickettsia Rickettsiaceae 21 0 0 * * Moraxellaceae 21 0 0 Bacteroidetes bacterium 21 0 0 WD317 Corynebacterium casei Corynebacterium Corynebacteriaceae 22 0 0 Corynebacterium Corynebacterium Corynebacteriaceae 22 0 0 marinum * * Flavobacteriaceae 22 0 0 Caulobacter segnis Caulobacter Caulobacteraceae 22 0 0 Lactobacillus gasseri Lactobacillus Lactobacillaceae 22 0 0 * Meiothermus Thermaceae 22 0 0 Rhizobium sp. NT-26 Rhizobium Rhizobiaceae 23 0 0 Bacillus coagulans Bacillus Bacillaceae 23 0 0 * Sphingomonas Sphingomonadaceae 23 0 0 * Brevibacterium Brevibacteriaceae 23 0 0 Nitrosomonas europaea Nitrosomonas Nitrosomonadaceae 23 0 0 Pseudomonas Pseudomonas Pseudomonadaceae 24 0 0 alkylphenolia Terrabacter sp. DBF63 Terrabacter Intrasporangiaceae 24 0 0 beta proteobacterium CB 24 0 0 Moraxella ovis Moraxella Moraxellaceae 24 0 0 Shewanella baltica Shewanella Shewanellaceae 24 0 0 Mycobacterium gilvum Mycobacterium Mycobacteriaceae 25 0 0 * Exiguobacterium 25 0 0 * Ochrobactrum Brucellaceae 25 0 0 Geodermatophilus Geodermatophilus Geodermatophilaceae 25 0 0 obscurus * Devosia Hyphomicrobiaceae 25 0 0 Moraxella osloensis Moraxella Moraxellaceae 25 0 0 Exiguobacterium sp. 11- Exiguobacterium 25 0 0 28 Nocardioides sp. JS614 Nocardioides Nocardioidaceae 26 0 0 Nocardioides sp. USM2 Nocardioides Nocardioidaceae 26 0 0 Burkholderia gladioli Burkholderia Burkholderiaceae 27 0 0 Renibacterium Renibacterium Micrococcaceae 27 0 0 salmoninarum Pseudomonas syringae Pseudomonas Pseudomonadaceae 27 0 0 group genomosp. 3 Bifidobacterium Bifidobacterium Bifidobacteriaceae 27 0 0 pseudolongum toluene-degrading 27 0 0 bacterium UCR 021t Corynebacterium imitans Corynebacterium Corynebacteriaceae 28 0 0 Corynebacterium callunae Corynebacterium Corynebacteriaceae 28 0 0 Bosea sp. WAO Bosea Bradyrhizobiaceae 28 0 0 Xanthobacter Xanthobacter Xanthobacteraceae 29 0 0 autotrophicus Corynebacterium Corynebacterium Corynebacteriaceae 29 0 0 diphtheriae Bacteroides Bacteroides Bacteroidaceae 29 0 0 thetaiotaomicron Caulobacter vibrioides Caulobacter Caulobacteraceae 29 0 0 * * Pasteurellaceae 30 0 0 Finegoldia magna Finegoldia Peptoniphilaceae 30 0 0 Anaerococcus prevotii Anaerococcus Peptoniphilaceae 30 0 0 Azorhizobium Azorhizobium Xanthobacteraceae 30 0 0 caulinodans Erysipelothrix Erysipelothrix Erysipelotrichaceae 31 0 0 rhusiopathiae Porphyromonas Porphyromonas Porphyromonadaceae 31 0 0 asaccharolytica * Sphingopyxis Sphingomonadaceae 31 0 0 Eubacterium rectale Eubacterium Eubacteriaceae 32 0 0 Acinetobacter venetianus Acinetobacter Moraxellaceae 33 0 0 Variovorax paradoxus Variovorax Comamonadaceae 33 0 0 Acinetobacter sp. ED45- Acinetobacter Moraxellaceae 33 0 0 25 Bradyrhizobium Bradyrhizobium Bradyrhizobiaceae 34 0 0 diazoefficiens Bradyrhizobium Bradyrhizobium Bradyrhizobiaceae 35 0 0 japonicum Megamonas hypermegale Megamonas Veillonellaceae 35 0 0 * Methylobacillus Methylophilaceae 35 0 0 Nocardiopsis alba Nocardiopsis Nocardiopsaceae 35 0 0 Modestobacter marinus Modestobacter Geodermatophilaceae 36 0 0 Corynebacterium Corynebacterium Corynebacteriaceae 37 0 0 doosanense Blastococcus saxobsidens Blastococcus Geodermatophilaceae 38 0 0 Anoxybacillus Anoxybacillus Bacillaceae 38 0 0 flavithermus Aeromonas caviae Aeromonas Aeromonadaceae 40 0 0 Bacillus subtilis Bacillus Bacillaceae 40 0 0 Elizabethkingia anophelis Elizabethkingia Flavobacteriaceae 40 0 0 Staphylococcus hominis Staphylococcus Staphylococcaceae 44 0 0 Ruminococcus bicirculans Ruminococcus Ruminococcaceae 45 0 0 Paracoccus marcusii Paracoccus Rhodobacteraceae 45 0 0 * Psychrobacter Moraxellaceae 46 0 0 Sphingomonas Sphingomonas Sphingomonadaceae 47 0 0 sanxanigenens * * Leuconostocaceae 47 0 0 * * Burkholderiaceae 49 0 0 Bacteroides vulgatus Bacteroides Bacteroidaceae 50 0 0 Rhodopseudomonas Rhodopseudomonas Bradyrhizobiaceae 51 0 0 palustris Pantoea agglomerans Pantoea Enterobacteriaceae 52 0 0 * * Sphingomonadaceae 53 0 0 Mycobacterium kansasii Mycobacterium Mycobacteriaceae 53 0 0 * Streptomyces Streptomycetaceae 55 0 0 Enterococcus faecalis Enterococcus Enterococcaceae 55 0 0 Acinetobacter sp. NFM2 Acinetobacter Moraxellaceae 57 0 0 Shewanella putrefaciens Shewanella Shewanellaceae 62 0 0 Bifidobacterium Bifidobacterium Bifidobacteriaceae 63 0 0 adolescentis * Bifidobacterium Bifidobacteriaceae 66 0 0 Porphyromonas gingivalis Porphyromonas Porphyromonadaceae 68 0 0 Neisseria meningitidis Neisseria Neisseriaceae 69 0 0 Rhodococcus Rhodococcus Nocardiaceae 70 0 0 pyridinivorans Aeromonas salmonicida Aeromonas Aeromonadaceae 73 0 0 Planococcus sp. PAMC Planococcus Planococcaceae 75 0 0 21323 Pseudomonas simiae Pseudomonas Pseudomonadaceae 81 0 0 Faecalibacterium Faecalibacterium Ruminococcaceae 95 0 0 prausnitzii Acinetobacter lwoffii Acinetobacter Moraxellaceae 142 0 0 Exiguobacterium sp. Exiguobacterium 146 0 0 N139 Streptococcus anginosus Streptococcus Streptococcaceae 156 0 0 Thauera sp. MZ1T Thauera Rhodocyclaceae 159 0 0 * Shewanella Shewanellaceae 226 0 0 * Aeromonas Aeromonadaceae 332 0 0 Staphylococcus aureus Staphylococcus Staphylococcaceae 407 0 0 Aeromonas hydrophila Aeromonas Aeromonadaceae 537 0 0 Aeromonas veronii Aeromonas Aeromonadaceae 563 0 0 Bifidobacterium breve Bifidobacterium Bifidobacteriaceae 892 0 0 Bacillus megaterium Bacillus Bacillaceae 2,131 0 0

SUPPLEMENTARY TABLE 2 Readcount DNA DNA DNA Patient species genus family NTC PC CSF Cryptococcus neoformans Filobasidiella Tremellaceae 0 104,135 0 Aspergillus niger Aspergillus Aspergillaceae 8 21,725 0 Toxoplasma gondii Toxoplasma Sarcocystidae 0 9,228 0 * Filobasidiella Tremellaceae 0 519 0 * * Sarcocystidae 0 266 0 Neospora caninum Neospora Sarcocystidae 0 50 0 Hammondia triffittae Hammondia Sarcocystidae 0 24 0 * * Aspergillaceae 0 20 0 * Hammondia Sarcocystidae 0 20 0 * Aspergillus Aspergillaceae 50 173 0 Hammondia hammondi Hammondia Sarcocystidae 7 20 0 Aspergillus fumigatus Aspergillus Aspergillaceae 0 15 0 Aspergillus awamori Aspergillus Aspergillaceae 0 7 0 Bartheletia paradoxa Bartheletia Bartheletiaceae 0 7 0 Cryptococcus gattii Filobasidiella Tremellaceae 0 7 0 Aspergillus kawachii Aspergillus Aspergillaceae 0 6 0 Spirometra erinaceieuropaei Spirometra Diphyllobothriidae 3 5 3 Mucor racemosus Mucor Mucoraceae 0 5 0 Aspergillus oryzae Aspergillus Aspergillaceae 0 4 0 Setosphaeria turcica Setosphaeria Pleosporaceae 0 4 0 Caenorhabditis remanei Caenorhabditis Rhabditidae 0 3 0 Aspergillus tubingensis Aspergillus Aspergillaceae 0 3 0 Anisakis simplex Anisakis Anisakidae 7 2 0 Penicillium solitum Penicillium Aspergillaceae 0 1 0 Plectosphaerella sp. 93 OA-2013 Plectosphaerella Plectosphaerellaceae 0 1 0 * * Tremellaceae 0 1 0 Malassezia globosa Malassezia Malasseziaceae 231 35 2 Gongylonema pulchrum Gongylonema Gongylonematidae 27 4 0 Wallemia sebi Wallemia 68 10 1 * * * 2,821 394 17 Candida parapsilosis Candida Debaryomycetaceae 323 45 6 Elaeophora elaphi Elaeophora Onchocercidae 77 10 1 Lichtheimia hongkongensis Lichtheimia Lichtheimiaceae 70 6 0 Brugia timori Brugia Onchocercidae 117 3 0 Parastrongyloides trichosuri Parastrongyloides Strongyloididae 293 5 1 Sordaria macrospora Sordaria Sordariaceae 7 0 1 Bipolaris sorokiniana Bipolaris Pleosporaceae 34 0 3 Malassezia restricta Malassezia Malasseziaceae 1 0 0 Cladosporium oxysporum Cladosporium Cladosporiaceae 1 0 0 Strongylocentrotus purpuratus Strongylocentrotus Strongylocentrotidae 5 0 0 Caenorhabditis elegans Caenorhabditis Rhabditidae 5 0 0 Penicillium griseoroseum Penicillium Aspergillaceae 7 0 0 Fusarium graminearum Fusarium Nectriaceae 7 0 0 Leptosphaeria biglobosa Leptosphaeria Leptosphaeriaceae 8 0 0 Chaetomium globosum Chaetomium Chaetomiaceae 8 0 0 Myceliophthora thermophila Myceliophthora Chaetomiaceae 8 0 0 Alternaria tenuissima Alternaria Pleosporaceae 8 0 0 * Saccharomyces Saccharomycetaceae 8 0 0 Aphanomyces euteiches Aphanomyces Saprolegniaceae 8 0 0 Ophiostoma piliferum Ophiostoma Ophiostomataceae 9 0 0 Rhodotorula taiwanensis Rhodotorula 9 0 0 Trichosporon domesticum Trichosporon 10 0 0 Phoma herbarum Phoma Didymellaceae 12 0 0 Elysiphe alphitoides Erysiphe Erysiphaceae 13 0 0 * Leptosphaeria Leptosphaeriaceae 15 0 0 Zymoseptoria tritici Zymoseptoria Mycosphaerellaceae 15 0 0 Debaryomyces hansenii Debaryomyces Debaryomycetaceae 15 0 0 Alternaria alternata Alternaria Pleosporaceae 16 0 0 * Phoma Didymellaceae 17 0 0 * Umbilicaria Umbilicariaceae 18 0 0 * Penicillium Aspergillaceae 20 0 0 Wuchereria bancrofti Wuchereria Onchocercidae 20 0 0 Cladosporium cladosporioides Cladosporium Cladosporiaceae 20 0 0 Saccharomyces bayanus Saccharomyces Saccharomycetaceae 20 0 0 Fusarium phaseoli Fusarium Nectriaceae 21 0 0 Dasybranchus sp. DH1 Dasybranchus Capitellidae 23 0 0 Neofusicoccum parvum Neofusicoccum Botryosphaeriaceae 30 0 0 * Cladosporium Cladosporiaceae 30 0 0 Penicillium rubens Penicillium Aspergillaceae 33 0 0 Ustilago maydis Ustilago Ustilaginaceae 38 0 0 Mucor circinelloides Mucor Mucoraceae 68 0 0 Albugo laibachii Albugo Albuginaceae 115 0 0

SUPPLEMENTARY TABLE 3 Normalized Reads Per million DNA Patient species genus family DNA PC CSF * Brucella Brucellaceae 0.0 15.7^(#) Brucella melitensis Brucella Brucellaceae 0.0 1.6 Streptococcus agalactiae Streptococcus Streptococcaceae 3,264.1^(#) 0.0 Klebsiella pneumoniae Klebsiella Enterobacteriaceae 52.9^(#) 0.0 * * Enterobacteriaceae 0.1 0.0 Escherichia coli Escherichia Enterobacteriaceae 0.0 0.0 Propionibacterium acnes Propionibacterium Propionibacteriaceae 0.0 0.0 * Streptococcus Streptococcaceae 8.7 0.0 * Klebsiella Enterobacteriaceae 41.6^(#) 0.0 * * * 0.0 0.0 Streptococcus suis Streptococcus Streptococcaceae 16.2^(#&) 0.0 Streptococcus salivarius Streptococcus Streptococcaceae 1.0 0.0 Streptococcus pyogenes Streptococcus Streptococcaceae 30.7^(#&) 0.0 Serratia marcescens Serratia Enterobacteriaceae 2.3 0.0 Lactococcus lactis Lactococcus Streptococcaceae 0.4 0.0 * Escherichia Enterobacteriaceae 0.0 0.0 Pseudomonas sp. TKP Pseudomonas Pseudomonadaceae 0.1 0.0 * Pseudomonas Pseudomonadaceae 0.2 0.0 Streptococcus macedonicus Streptococcus Streptococcaceae 18.2^(#&) 0.0 Enterobacter cloacae Enterobacter Enterobacteriaceae 0.4 0.0 Klebsiella oxytoca Klebsiella Enterobacteriaceae 0.7 0.0 Thermoanaerobacterium Thermoanaerobacterium Thermoanaerobacterales 0.2 0.0 thermosaccharolyticum Family III. Incertae Sedis Streptococcus sp. VT 162 Streptococcus Streptococcaceae 0.6 0.0 Enterococcus faecium Enterococcus Enterococcaceae 0.1 0.1 Pseudomonas protegens Pseudomonas Pseudomonadaceae 3.0 0.1 Micrococcus luteus Micrococcus Micrococcaceae 0.0 0.0 Streptococcus infantarius Streptococcus Streptococcaceae 0.5 0.0 Staphylococcus epidermidis Staphylococcus Staphylococcaceae 0.0 0.0 Enterobacter asburiae Enterobacter Enterobacteriaceae 0.2 0.0 Cupriavidus metallidurans Cupriavidus Burkholderiaceae 0.0 0.0 Pseudomonas putida Pseudomonas Pseudomonadaceae 0.1 0.0 Enterococcus casseliflavus Enterococcus Enterococcaceae 0.1 0.1 Streptococcus mitis Streptococcus Streptococcaceae 0.6 0.0 Streptococcus equi Streptococcus Streptococcaceae 5.2 0.0 Salmonella enterica Salmonella Enterobacteriaceae 4.9 0.0 Streptococcus oralis Streptococcus Streptococcaceae 0.5 0.0 Klebsiella variicola Klebsiella Enterobacteriaceae 4.3 0.0 Burkholderia lata Burkholderia Burkholderiaceae 0.3 0.0 Pseudomonas stutzeri Pseudomonas Pseudomonadaceae 0.0 0.0 Streptococcus pneumoniae Streptococcus Streptococcaceae 0.4 0.0 Streptococcus dysgalactiae Streptococcus Streptococcaceae 3.7 0.0 * Burkholderia Burkholderiaceae 0.1 0.0 Pseudomonas fluorescens Pseudomonas Pseudomonadaceae 0.4 0.0 Acinetobacter guillouiae Acinetobacter Moraxellaceae 0.1 0.0 Veillonella parvula Veillonella Veillonellaceae 0.1 0.0 Xanthomonas campestris Xanthomonas Xanthomonadaceae 0.1 0.1 Exiguobacterium sp. AT1b Exiguobacterium 0.1 0.0 Pseudomonas Pseudomonas Pseudomonadaceae 0.0 0.0 pseudoalcaligenes Streptococcus pasteurianus Streptococcus Streptococcaceae 2.8 0.0 Rothia dentocariosa Rothia Micrococcaceae 0.1 0.0 * * Streptococcaceae 2.6 0.0 Delftia acidovorans Delftia Comamonadaceae 0.5 0.0 * Propionibacterium Propionibacteriaceae 0.0 0.0 Acinetobacter baumannii Acinetobacter Moraxellaceae 0.0 0.0 Streptococcus parasanguinis Streptococcus Streptococcaceae 0.0 0.0 Pseudomonas sp. WCS374 Pseudomonas Pseudomonadaceae 0.8 0.0 Enterobacter aerogenes Enterobacter Enterobacteriaceae 2.3 0.0 Stenotrophomonas maltophilia Stenotrophomonas Xanthomonadaceae 0.2 0.0 Alicyclobacillus Alicyclobacillus Alicyclobacillaceae 0.2 0.0 acidocaldarius Haemophilus influenzae Haemophilus Pasteurellaceae 0.1 0.0 Rothia mucilaginosa Rothia Micrococcaceae 0.0 0.0 Staphylococcus xylosus Staphylococcus Staphylococcaceae 0.2 0.0 Acidovorax sp. JS42 Acidovorax Comamonadaceae 0.1 0.0 * Rahnella Enterobacteriaceae 1.8 0.0 Streptococcus Streptococcus Streptococcaceae 0.3 0.0 pseudopneumoniae Streptococcus thermophilus Streptococcus Streptococcaceae 0.1 0.0 Pseudomonas aeruginosa Pseudomonas Pseudomonadaceae 0.1 0.0 Corynebacterium Corynebacterium Corynebacteriaceae 1.6 0.0 kroppenstedtii Staphylococcus haemolyticus Staphylococcus Staphylococcaceae 0.0 0.0 Serratia sp. SCBI Serratia Enterobacteriaceae 1.5 0.0 * Enterococcus Enterococcaceae 0.1 0.0 Staphylococcus wameri Staphylococcus Staphylococcaceae 0.0 0.0 Burkholderia cepacia Burkholderia Burkholderiaceae 0.1 0.1 * Corynebacterium Corynebacteriaceae 0.0 0.0 Comamonas testosteroni Comamonas Comamonadaceae 0.0 0.0 Bifidobacterium thermophilum Bifidobacterium Bifidobacteriaceae 1.3 0.0 * Lactobacillus Lactobacillaceae 0.2 0.0 Thermoanaerobacterium Thermoanaerobacterium Thermoanaerobacterales 0.1 0.0 xylanolyticum Family III. Incertae Sedis Ralstonia pickettii Ralstonia Burkholderiaceae 0.1 0.0 Meiothermus ruber Meiothermus Thermaceae 0.0 0.0 Acidovorax ebreus Acidovorax Comamonadaceae 0.7 0.0 Micrococcus sp. V7 Micrococcus Micrococcaceae 0.2 0.0 Leuconostoc mesenteroides Leuconostoc Leuconostocaceae 0.0 0.0 Bacillus halodurans Bacillus Bacillaceae 0.7 0.0 Corynebacterium variabile Corynebacterium Corynebacteriaceae 0.3 0.0 * Acidovorax Comamonadaceae 0.1 0.0 Bifidobacterium bifidum Bifidobacterium Bifidobacteriaceae 0.0 0.0 Rhizobium sp. IRBG74 Rhizobium Rhizobiaceae 0.0 0.0 * Acinetobacter Moraxellaceae 0.0 0.0 Acinetobacter calcoaceticus Acinetobacter Moraxellaceae 1.0 0.0 Chroococcidiopsis thermalis Chroococcidiopsis 0.2 0.0 Streptococcus sanguinis Streptococcus Streptococcaceae 0.1 0.0 Acidovorax sp. KKS102 Acidovorax Comamonadaceae 0.1 0.0 Raoultella ornithinolytica Raoultella Enterobacteriaceae 0.1 0.0 Ochrobactrum anthropi Ochrobactrum Brucellaceae 0.0 0.0 Lactobacillus johnsonii Lactobacillus Lactobacillaceae 0.6 0.1 Methylobacterium populi Methylobacterium Methylobacteriaceae 0.2 0.0 Rhodococcus equi Rhodococcus Nocardiaceae 0.1 0.0 Lactobacillus helveticus Lactobacillus Lactobacillaceae 0.9 0.0 * Serratia Enterobacteriaceae 0.2 0.0 Burkholderia cenocepacia Burkholderia Burkholderiaceae 0.3 0.0 * Staphylococcus Staphylococcaceae 0.0 0.0 Enterobacter sp. R4-368 Enterobacter Enterobacteriaceae 0.8 0.0 Propionibacterium Propionibacterium Propionibacteriaceae 0.1 0.0 propionicum Streptococcus gordonii Streptococcus Streptococcaceae 0.1 0.0 Corynebacterium singulare Corynebacterium Corynebacteriaceae 0.1 0.0 Burkholderia ambifaria Burkholderia Burkholderiaceae 0.7 0.1 * Micrococcus Micrococcaceae 0.0 0.0 Fervidobacterium nodosum Fervidobacterium Thermotogaceae 0.0 0.0 Aeromonas media Aeromonas Aeromonadaceae 0.0 0.0 Cronobacter sakazakii Cronobacter Enterobacteriaceae 0.7 0.0 Myroides profundi Myroides Flavobacteriaceae 0.7 0.0 Methylobacterium oryzae Methylobacterium Methylobacteriaceae 0.7 0.0 * Xanthomonas Xanthomonadaceae 0.3 0.0 Thermoanaerobacterium Thermoanaerobacterium Thermoanaerobacterales 0.1 0.0 saccharolyticum Family III. Incertae Sedis Pseudomonas mendocina Pseudomonas Pseudomonadaceae 0.0 0.0 Corynebacterium Corynebacterium Corynebacteriaceae 0.0 0.0 ureicelerivorans Lactobacillus crispatus Lactobacillus Lactobacillaceae 0.0 0.0 Alicycliphilus denitrificans Alicycliphilus Comamonadaceae 0.0 0.0 Gardnerella vaginalis Gardnerella Bifidobacteriaceae 0.0 0.0 * Gemella 0.6 0.0 * Ralstonia Burkholderiaceae 0.6 0.0 Eggerthella lenta Eggerthella Coriobacteriaceae 0.6 0.0 * * Rhizobiaceae 0.5 0.0 Prevotella denticola Prevotella Prevotellaceae 0.2 0.0 Prevotella intermedia Prevotella Prevotellaceae 0.2 0.0 Psychrobacter sp. PRwf-1 Psychrobacter Moraxellaceae 0.1 0.0 Azospira oryzae Azospira Rhodocyclaceae 0.1 0.0 Acinetobacter haemolyticus Acinetobacter Moraxellaceae 0.1 0.0 * Delftia Comamonadaceae 0.1 0.0 Burkholderia contaminans Burkholderia Burkholderiaceae 0.4 0.2 Arthrobacter arilaitensis Arthrobacter Micrococcaceae 0.1 0.1 Dermacoccus Dermacoccus Dermacoccaceae 0.1 0.1 nishinomiyaensis Pantoea ananatis Pantoea Enterobacteriaceae 0.1 0.0 Staphylococcus saprophyticus Staphylococcus Staphylococcaceae 0.0 0.0 Staphylococcus pasteuri Staphylococcus Staphylococcaceae 0.0 0.0 Rahnella aquatilis Rahnella Enterobacteriaceae 0.5 0.0 Rahnella sp. Y9602 Rahnella Enterobacteriaceae 0.5 0.0 Campylobacter concisus Campylobacter Campylobacteraceae 0.5 0.0 Geobacillus sp. WCH70 Geobacillus Bacillaceae 0.5 0.0 * Frankia Frankiaceae 0.5 0.0 Lactobacillus casei Lactobacillus Lactobacillaceae 0.5 0.0 Thiomonas intermedia Thiomonas 0.5 0.0 Streptococcus gallolyticus Streptococcus Streptococcaceae 0.5 0.0 Thioalkalivibrio sulfidiphilus Thioalkalivibrio Ectothiorhodospiraceae 0.3 0.0 * Bradyrhizobium Bradyrhizobiaceae 0.2 0.0 Bifidobacterium longum Bifidobacterium Bifidobacteriaceae 0.1 0.0 Corynebacterium falsenii Corynebacterium Corynebacteriaceae 0.1 0.0 Delftia sp. Cs1-4 Delftia Comamonadaceae 0.1 0.0 Acinetobacter sp. M131 Acinetobacter Moraxellaceae 0.0 0.0 Prevotella melaninogenica Prevotella Prevotellaceae 0.0 0.0 * * Comamonadaceae 0.0 0.0 Leuconostoc carnosum Leuconostoc Leuconostocaceae 0.2 0.0 Pectobacterium carotovorum Pectobacterium Enterobacteriaceae 0.0 0.0 * Myroides Flavobacteriaceae 0.4 0.0 * Erwinia Enterobacteriaceae 0.4 0.0 Gordonia sp. KTR9 Gordonia Gordoniaceae 0.4 0.0 Paenibacillus sp. FSL R7-0273 Paenibacillus Paenibacillaceae 0.4 0.0 Paracoccus sp. N81106 Paracoccus Rhodobacteraceae 0.4 0.0 Sphingobium fuliginis Sphingobium Sphingomonadaceae 0.4 0.0 * * 0.4 0.0 * Geobacillus Bacillaceae 0.2 0.0 Pseudomonas sp. VLB120 Pseudomonas Pseudomonadaceae 0.2 0.0 Pelagibacterium halotolerans Pelagibacterium Hyphomicrobiaceae 0.1 0.0 Streptococcus intermedius Streptococcus Streptococcaceae 0.1 0.0 Propionibacterium Propionibacterium Propionibacteriaceae 0.1 0.0 freudenreichii * Enterobacter Enterobacteriaceae 0.0 0.0 Nakamurella multipartita Nakamurella Nakamurellaceae 0.1 0.1 Haemophilus parasuis Haemophilus Pasteurellaceae 0.3 0.0 Fusobacterium nucleatum Fusobacterium Fusobacteriaceae 0.3 0.0 Citrobacter freundii Citrobacter Enterobacteriaceae 0.3 0.0 Ruminococcus sp. SR1/5 Ruminococcus Ruminococcaceae 0.2 0.0 Pseudoxanthomonas spadix Pseudoxanthomonas Xanthomonadaceae 0.2 0.0 Lactococcus garvieae Lactococcus Streptococcaceae 0.1 0.0 Neisseria elongata Neisseria Neisseriaceae 0.1 0.0 Acidovorax citrulli Acidovorax Comamonadaceae 0.1 0.0 Novosphingobium Novosphingobium Sphingomonadaceae 0.1 0.0 pentaromativorans Citrobacter koseri Citrobacter Enterobacteriaceae 0.1 0.0 Methylobacterium aquaticum Methylobacterium Methylobacteriaceae 0.1 0.0 Pseudomonas denitrificans Pseudomonas Pseudomonadaceae 0.1 0.0 Rhodococcus erythropolis Rhodococcus Nocardiaceae 0.0 0.0 Lactobacillus reuteri Lactobacillus Lactobacillaceae 0.0 0.0 Bacteroides fragilis Bacteroides Bacteroidaceae 0.0 0.0 Lactobacillus plantarum Lactobacillus Lactobacillaceae 0.0 0.0 * Bacillus Bacillaceae 0.0 0.0 Pseudomonas rhizosphaerae Pseudomonas Pseudomonadaceae 0.2 0.1 Achromobacter xylosoxidans Achromobacter Alcaligenaceae 0.0 0.0 Lactobacillus amylovorus Lactobacillus Lactobacillaceae 0.2 0.0 Propionibacterium Propionibacterium Propionibacteriaceae 0.2 0.0 acidipropionici Leuconostoc gelidum Leuconostoc Leuconostocaceae 0.2 0.0 Weissella thailandensis Weissella Leuconostocaceae 0.2 0.0 Pandoraea sp. RB-44 Pandoraea Burkholderiaceae 0.2 0.0 Escherichia vulneris Escherichia Enterobacteriaceae 0.2 0.0 Yersinia intermedia Yersinia Enterobacteriaceae 0.2 0.0 Flavobacteriaceae bacterium Flavobacteriaceae 0.2 0.0 3519-10 * Rhodococcus Nocardiaceae 0.2 0.0 Streptococcus sp. (N1) Streptococcus Streptococcaceae 0.2 0.0 * Methylobacterium Methylobacteriaceae 0.1 0.0 Sphingomonas sp. MM-1 Sphingomonas Sphingomonadaceae 0.1 0.0 Rhizobium etli Rhizobium Rhizobiaceae 0.1 0.0 * Agrobacterium Rhizobiaceae 0.0 0.0 Thermus scotoductus Thermus Thermaceae 0.0 0.0 Methylobacterium extorquens Methylobacterium Methylobacteriaceae 0.0 0.0 Streptococcus sp. I-P16 Streptococcus Streptococcaceae 0.0 0.0 Pantoea sp. PSNIH1 Pantoea Enterobacteriaceae 0.0 0.0 Methylobacterium Methylobacterium Methylobacteriaceae 0.0 0.0 radiotolerans [Ruminococcus] torques Blautia Lachnospiraceae 0.1 0.4 Lactobacillus sakei Lactobacillus Lactobacillaceae 0.1 0.1 * * Rhodocyclaceae 0.1 0.1 Bacillus licheniformis Bacillus Bacillaceae 0.0 0.1 Corynebacterium accolens Corynebacterium Corynebacteriaceae 0.1 0.0 Corynebacterium sp. ATCC Corynebacterium Corynebacteriaceae 0.1 0.0 6931 Serratia symbiotica Serratia Enterobacteriaceae 0.1 0.0 Lactobacillus delbrueckii Lactobacillus Lactobacillaceae 0.1 0.0 Bacillus sp. YP1 Bacillus Bacillaceae 0.1 0.0 Klebsiella sp. PG122E Klebsiella Enterobacteriaceae 0.1 0.0 * Rathayibacter Microbacteriaceae 0.1 0.0 Pseudoalteromonas sp. P30 Pseudoalteromonas Pseudoalteromonadaceae 0.1 0.0 Staphylococcus sp. CDC25 Staphylococcus Staphylococcaceae 0.1 0.0 Corynebacterium resistens Corynebacterium Corynebacteriaceae 0.1 0.0 Shigella dysenteriae Shigella Enterobacteriaceae 0.1 0.0 * * Xanthomonadaceae 0.1 0.0 Agrobacterium fabrum Agrobacterium Rhizobiaceae 0.0 0.0 Gordonia polyisoprenivorans Gordonia Gordoniaceae 0.0 0.0 Pseudomonas balearica Pseudomonas Pseudomonadaceae 0.0 0.0 * * Pseudomonadaceae 0.0 0.0 Ruminococcus bromii Ruminococcus Ruminococcaceae 0.0 0.0 Brachybacterium faecium Brachybacterium Dermabacteraceae 0.0 0.0 Acinetobacter johnsonii Acinetobacter Moraxellaceae 0.0 0.0 Micrococcus sp. A1 Micrococcus Micrococcaceae 0.0 0.0 Filifactor alocis Filifactor Peptostreptococcaceae 0.0 0.0 Pantoea vagans Pantoea Enterobacteriaceae 0.0 0.0 Haemophilus parainfluenzae Haemophilus Pasteurellaceae 0.0 0.0 Pantoea rwandensis Pantoea Enterobacteriaceae 0.0 0.3 Corynebacterium Corynebacterium Corynebacteriaceae 0.0 0.2 vitaeruminis Pseudomonas poae Pseudomonas Pseudomonadaceae 0.0 0.1 * * Brucellaceae 0.0 0.1 Lactobacillus fermentum Lactobacillus Lactobacillaceae 0.0 0.1 Anabaena variabilis Anabaena Nostocaceae 0.0 0.1 Sphingobacterium sp. ML3W Sphingobacterium Sphingobacteriaceae 0.0 0.1 sugarcane isolate 74-1 0.0 0.1 * * Geodermatophilaceae 0.0 0.1 Megasphaera elsdenii Megasphaera Veillonellaceae 0.0 0.1 Pseudoxanthomonas Pseudoxanthomonas Xanthomonadaceae 0.0 0.1 suwonensis Corynebacterium Corynebacterium Corynebacteriaceae 0.0 0.0 glutamicum Sphingomonas taxi Sphingomonas Sphingomonadaceae 0.0 0.0 Pseudomonas graminis Pseudomonas Pseudomonadaceae 0.0 0.0 Bradyrhizobium sp. BTAi1 Bradyrhizobium Bradyrhizobiaceae 0.0 0.0 Enterococcus hirae Enterococcus Enterococcaceae 0.0 0.0 Corynebacterium sp. L2-79-05 Corynebacterium Corynebacteriaceae 0.0 0.0 Arthrobacter Arthrobacter Micrococcaceae 0.0 0.0 phenanthrenivorans Corynebacterium maris Corynebacterium Corynebacteriaceae 0.0 0.0 Gordonia bronchialis Gordonia Gordoniaceae 0.0 0.0 Kytococcus sedentarius Kytococcus Dermacoccaceae 0.0 0.0 Kosakonia cowanii Kosakonia Enterobacteriaceae 0.0 0.0 Xenorhabdus bovienii Xenorhabdus Enterobacteriaceae 0.0 0.0 Paracoccus haeundaensis Paracoccus Rhodobacteraceae 0.0 0.0 Methylobacterium sp. 238 Methylobacterium Methylobacteriaceae 0.0 0.0 Acinetobacter sp. BW3 Acinetobacter Moraxellaceae 0.0 0.0 Aeromonas sobria Aeromonas Aeromonadaceae 0.0 0.0 Bacillus lehensis Bacillus Bacillaceae 0.0 0.0 Ralstonia solanacearum Ralstonia Burkholderiaceae 0.0 0.0 Citrobacter sp. FPO3 Citrobacter Enterobacteriaceae 0.0 0.0 Citrobacter sp. I91-3 Citrobacter Enterobacteriaceae 0.0 0.0 Erwinia amylovora Erwinia Enterobacteriaceae 0.0 0.0 Klebsiella milletis Klebsiella Enterobacteriaceae 0.0 0.0 Salmonella bongori Salmonella Enterobacteriaceae 0.0 0.0 Serratia grimesii Serratia Enterobacteriaceae 0.0 0.0 Yersinia pestis Yersinia Enterobacteriaceae 0.0 0.0 * Enterobacteriaceae 0.0 0.0 Lactobacillus brevis Lactobacillus Lactobacillaceae 0.0 0.0 Kocuria sp. starX Kocuria Micrococcaceae 0.0 0.0 Acinetobacter sp. 26 Acinetobacter Moraxellaceae 0.0 0.0 Peptoclostridium difficile Peptoclostridium Peptostreptococcaceae 0.0 0.0 Sorangium cellulosum Sorangium Polyangiaceae 0.0 0.0 Pseudomonas sp. NSi14 Pseudomonas Pseudomonadaceae 0.0 0.0 Synergistetes oral clone 03 5 0.0 0.0 D05 bacterium EBAD26 0.0 0.0 bacterium NLAE-zl-G351 0.0 0.0 rumen bacterium enrichment 0.0 0.0 culture clone Y74 unidentified marine 0.0 0.0 bacterioplankton Escherichia albertii Escherichia Enterobacteriaceae 0.0 0.0 Mobiluncus curtisii Mobiluncus Actinomycetaceae 0.0 0.0 * Caulobacter Caulobacteraceae 0.0 0.0 Methylotenera versatilis Methylotenera Methylophilaceae 0.0 0.0 Propionibacterium sp. Propionibacterium Propionibacteriaceae 0.0 0.0 KPL1849 bacterium EBAD25 0.0 0.0 Bacillus sp. Pc3 Bacillus Bacillaceae 0.0 0.0 Acinetobacter sp. EVA14 Acinetobacter Moraxellaceae 0.0 0.0 Agrobacterium sp. Agrobacterium Rhizobiaceae 0.0 0.0 Alistipes shahii Alistipes Rikenellaceae 0.0 0.0 * Thermus Thermaceae 0.0 0.0 * Methylibium 0.0 0.0 butyrate-producing 0.0 0.0 bacterium SSC/2 Escherichia fergusonii Escherichia Enterobacteriaceae 0.0 0.0 Enterobacter sp. Ni15 Enterobacter Enterobacteriaceae 0.0 0.0 Capnocytophaga ochracea Capnocytophaga Flavobacteriaceae 0.0 0.0 Thauera sp. 6NLG Thauera Rhodocyclaceae 0.0 0.0 Desulfovibrio alaskensis Desulfovibrio Desulfovibrionaceae 0.0 0.0 Variovorax sp. Alb14 Variovorax Comamonadaceae 0.0 0.0 * Shigella Enterobacteriaceae 0.0 0.0 * * Micromonosporaceae 0.0 0.0 * Micromonospora Micromonosporaceae 0.0 0.0 Thermobifida fusca Thermobifida Nocardiopsaceae 0.0 0.0 Turneriella parva Turneriella Leptospiraceae 0.0 0.0 [Clostridium] sticklandii Peptoclostridium Peptostreptococcaceae 0.0 0.0 Acinetobacter sp. Ooi24 Acinetobacter Moraxellaceae 0.0 0.0 Ochrobactrum sp. SJY1 Ochrobactrum Brucellaceae 0.0 0.0 Carnobacterium sp. WN1359 Carnobacterium Carnobacteriaceae 0.0 0.0 Iamia majanohamensis Iamia Iamiaceae 0.0 0.0 Saccharomonospora viridis Saccharomonospora Pseudonocardiaceae 0.0 0.0 Rhizobium sp. Rhizobium Rhizobiaceae 0.0 0.0 Staphylococcus sp. CDC3 Staphylococcus Staphylococcaceae 0.0 0.0 Shigella sonnei Shigella Enterobacteriaceae 0.0 0.0 Pseudomonas syringae Pseudomonas Pseudomonadaceae 0.0 0.0 Burkholderia vietnamiensis Burkholderia Burkholderiaceae 0.0 0.0 Shigella boydii Shigella Enterobacteriaceae 0.0 0.0 Bacillus weihenstephanensis Bacillus Bacillaceae 0.0 0.0 Erythrobacter litoralis Erythrobacter Erythrobacteraceae 0.0 0.0 Pseudoalteromonas Pseudoalteromonas Pseudoalteromonadaceae 0.0 0.0 haloplanktis Pseudomonas sp. FGI182 Pseudomonas Pseudomonadaceae 0.0 0.0 * Rhizobium Rhizobiaceae 0.0 0.0 * Rickettsia Rickettsiaceae 0.0 0.0 Sphingobium yanoikuyae Sphingobium Sphingomonadaceae 0.0 0.0 Stenotrophomonas Stenotrophomonas Xanthomonadaceae 0.0 0.0 rhizophila * Leuconostoc Leuconostocaceae 0.0 0.0 Aquincola tertiaricarbonis Aquincola 0.0 0.0 Nocardiopsis dassonvillei Nocardiopsis Nocardiopsaceae 0.0 0.0 Carnobacterium Carnobacterium Carnobacteriaceae 0.0 0.0 maltaromaticum * Haemophilus Pasteurellaceae 0.0 0.0 Bordetella parapertussis Bordetella Alcaligenaceae 0.0 0.0 * Dietzia Dietziaceae 0.0 0.0 Shewanella sp. W3-18-1 Shewanella Shewanellaceae 0.0 0.0 Sphingomonas sp. NP5 Sphingomonas Sphingomonadaceae 0.0 0.0 Staphylococcus gallinarum Staphylococcus Staphylococcaceae 0.0 0.0 Micavibrio aeruginosavorus Micavibrio 0.0 0.0 Paracoccus denitrificans Paracoccus Rhodobacteraceae 0.0 0.0 [Cellvibrio] gilvus Cellulomonas Cellulomonadaceae 0.0 0.0 Corynebacterium jeikeium Corynebacterium Corynebacteriaceae 0.0 0.0 * * Staphylococcaceae 0.0 0.0 Meiothermus silvanus Meiothermus Thermaceae 0.0 0.0 Asticcacaulis excentricus Asticcacaulis Caulobacteraceae 0.0 0.0 * Atopobium Coriobacteriaceae 0.0 0.0 Streptococcus constellatus Streptococcus Streptococcaceae 0.0 0.0 Microcystis aeruginosa Microcystis 0.0 0.0 agricultural soil bacterium 0.0 0.0 SC-I-13 [Ruminococcus] obeum Blautia Lachnospiraceae 0.0 0.0 Thermus thermophilus Thermus Thermaceae 0.0 0.0 Shigella flexneri Shigella Enterobacteriaceae 0.0 0.0 * Mycobacterium Mycobacteriaceae 0.0 0.0 Pseudomonas savastanoi Pseudomonas Pseudomonadaceae 0.0 0.0 Staphylococcus capitis Staphylococcus Staphylococcaceae 0.0 0.0 * Cupriavidus Burkholderiaceae 0.0 0.0 Dyadobacter fermentans Dyadobacter Cytophagaceae 0.0 0.0 Dietzia sp. CQ4 Dietzia Dietziaceae 0.0 0.0 * * Methylophilaceae 0.0 0.0 * Neisseria Neisseriaceae 0.0 0.0 Corynebacterium Corynebacterium Corynebacteriaceae 0.0 0.0 aurimucosum * 0.0 0.0 Pseudomonas fulva Pseudomonas Pseudomonadaceae 0.0 0.0 Chromohalobacter salexigens Chromohalobacter Halomonadaceae 0.0 0.0 Brevundimonas diminuta Brevundimonas Caulobacteraceae 0.0 0.0 Streptococcus lutetiensis Streptococcus Streptococcaceae 0.0 0.0 Bordetella petrii Bordetella Alcaligenaceae 0.0 0.0 Erythrobacter sp. JP13.1 Erythrobacter Erythrobacteraceae 0.0 0.0 Methylobacillus glycogenes Methylobacillus Methylophilaceae 0.0 0.0 Candidatus Rhodoluna Candidatus Microbacteriaceae 0.0 0.0 lacicola Rhodoluna Arthrobacter sp. JBH1 Arthrobacter Micrococcaceae 0.0 0.0 Aggregatibacter aphrophilus Aggregatibacter Pasteurellaceae 0.0 0.0 Thauera sp. B4 Thauera Rhodocyclaceae 0.0 0.0 Lysobacter dokdonensis Lysobacter Xanthomonadaceae 0.0 0.0 Clostridiales genomosp. 0.0 0.0 BVAB3 Streptococcus sp. I-G2 Streptococcus Streptococcaceae 0.0 0.0 Pseudomonas mandelii Pseudomonas Pseudomonadaceae 0.0 0.0 Bradyrhizobium sp. S23321 Bradyrhizobium Bradyrhizobiaceae 0.0 0.0 Phenylobacterium zucineum Phenylobacterium Caulobacteraceae 0.0 0.0 Pseudomonas mosselii Pseudomonas Pseudomonadaceae 0.0 0.0 Staphylococcus lugdunensis Staphylococcus Staphylococcaceae 0.0 0.0 Proteus mirabilis Proteus Enterobacteriaceae 0.0 0.0 * * Neisseriaceae 0.0 0.0 Arthrobacter sp. J3-40 Arthrobacter Micrococcaceae 0.0 0.0 * Pantoea Enterobacteriaceae 0.0 0.0 Corynebacterium efficiens Corynebacterium Corynebacteriaceae 0.0 0.0 * Halomonas Halomonadaceae 0.0 0.0 Trueperella pyogenes Trueperella Actinomycetaceae 0.0 0.0 Streptomyces coelicolor Streptomyces Streptomycetaceae 0.0 0.0 Kocuria rhizophila Kocuria Micrococcaceae 0.0 0.0 Bacillus cereus Bacillus Bacillaceae 0.0 0.0 Tannerella forsythia Tannerella Porphyromonadaceae 0.0 0.0 * Alkalibacterium Carnobacteriaceae 0.0 0.0 Atopobium parvulum Atopobium Coriobacteriaceae 0.0 0.0 Serinicoccus profundi Serinicoccus Intrasporangiaceae 0.0 0.0 * Leptotrichia Leptotrichiaceae 0.0 0.0 * * Planococcaceae 0.0 0.0 Planomicrobium Planomicrobium Planococcaceae 0.0 0.0 okeanokoites Jannaschia sp. CCS1 Jannaschia Rhodobacteraceae 0.0 0.0 Paracoccus aestuarii Paracoccus Rhodobacteraceae 0.0 0.0 Rhodobacter blasticus Rhodobacter Rhodobacteraceae 0.0 0.0 Agrobacterium tumefaciens Agrobacterium Rhizobiaceae 0.0 0.0 Shewanella sp. ANA-3 Shewanella Shewanellaceae 0.0 0.0 Pseudomonas cichorii Pseudomonas Pseudomonadaceae 0.0 0.0 Halomonas sp. A3H3 Halomonas Halomonadaceae 0.0 0.0 Serratia liquefaciens Serratia Enterobacteriaceae 0.0 0.0 Sphingopyxis alaskensis Sphingopyxis Sphingomonadaceae 0.0 0.0 * Brevundimonas Caulobacteraceae 0.0 0.0 Deinococcus deserti Deinococcus Deinococcaceae 0.0 0.0 Desulfovibrio vulgaris Desulfovibrio Desulfovibrionaceae 0.0 0.0 Propionibacterium sp. Propionibacterium Propionibacteriaceae 0.0 0.0 NTS31307302 Paracoccus alcaliphilus Paracoccus Rhodobacteraceae 0.0 0.0 Vibrio parahaemolyticus Vibrio Vibrionaceae 0.0 0.0 Candidatus Saccharibacteria 0.0 0.0 oral taxon TM7x Geitlerinema sp. PCC 7407 Geitlerinema 0.0 0.0 * Actinomyces Actinomycetaceae 0.0 0.0 Brevundimonas vesicularis Brevundimonas Caulobacteraceae 0.0 0.0 Acinetobacter sp. YS0810 Acinetobacter Moraxellaceae 0.0 0.0 * Prevotella Prevotellaceae 0.0 0.0 Methyloceanibacter Methyloceanibacter 0.0 0.0 caenitepidi Leuconostoc citreum Leuconostoc Leuconostocaceae 0.0 0.0 * Bacteroides Bacteroidaceae 0.0 0.0 Pseudomonas alcaligenes Pseudomonas Pseudomonadaceae 0.0 0.0 Methylibium petroleiphilum Methylibium 0.0 0.0 Moraxella catarrhalis Moraxella Moraxellaceae 0.0 0.0 Sphingopyxis sp. Kp5.2 Sphingopyxis Sphingomonadaceae 0.0 0.0 Pandoraea apista Pandoraea Burkholderiaceae 0.0 0.0 * Cellulomonas Cellulomonadaceae 0.0 0.0 Olsenella uli Olsenella Coriobacteriaceae 0.0 0.0 Acinetobacter oleivorans Acinetobacter Moraxellaceae 0.0 0.0 Sphingomonas sp. 133 Sphingomonas Sphingomonadaceae 0.0 0.0 Meiothermus taiwanensis Meiothermus Thermaceae 0.0 0.0 Deinococcus geothermalis Deinococcus Deinococcaceae 0.0 0.0 Lactobacillus Lactobacillus Lactobacillaceae 0.0 0.0 sanfranciscensis Acinetobacter sp. LUH5605 Acinetobacter Moraxellaceae 0.0 0.0 Staphylococcus simulans Staphylococcus Staphylococcaceae 0.0 0.0 Arsenophonus nasoniae Arsenophonus Enterobacteriaceae 0.0 0.0 Buchnera aphidicola Buchnera Enterobacteriaceae 0.0 0.0 Weissella koreensis Weissella Leuconostocaceae 0.0 0.0 Psychrobacter sp. G Psychrobacter Moraxellaceae 0.0 0.0 Amycolicicoccus subflavus Amycolicicoccus Mycobacteriaceae 0.0 0.0 Staphylococcus hyicus Staphylococcus Staphylococcaceae 0.0 0.0 Morganella morganii Morganella Enterobacteriaceae 0.0 0.0 Kyrpidia tusciae Kyrpidia Alicyclobacillaceae 0.0 0.0 Ramlibacter tataouinensis Ramlibacter Comamonadaceae 0.0 0.0 Weeksella virosa Weeksella Flavobacteriaceae 0.0 0.0 Acinetobacter junii Acinetobacter Moraxellaceae 0.0 0.0 Acinetobacter sp. 26132 Acinetobacter Moraxellaceae 0.0 0.0 Mycobacterium abscessus Mycobacterium Mycobacteriaceae 0.0 0.0 Neisseria gonorrhoeae Neisseria Neisseriaceae 0.0 0.0 Sphingomonas wittichii Sphingomonas Sphingomonadaceae 0.0 0.0 Bacteroides dorei Bacteroides Bacteroidaceae 0.0 0.0 Corynebacterium Corynebacterium Corynebacteriaceae 0.0 0.0 halotolerans Pelomonas aquatica Pelomonas Comamonadaceae 0.0 0.0 Janibacter sp. TYM3221 Janibacter Intrasporangiaceae 0.0 0.0 * Arthrobacter Micrococcaceae 0.0 0.0 Mycobacterium gordonae Mycobacterium Mycobacteriaceae 0.0 0.0 Pimelobacter simplex Pimelobacter Nocardioidaceae 0.0 0.0 Pseudomonas sp. OM2164 Pseudomonas Pseudomonadaceae 0.0 0.0 Streptomyces albus Streptomyces Streptomycetaceae 0.0 0.0 Halomonas halocynthiae Halomonas Halomonadaceae 0.0 0.0 Nitrosomonas sp. AL212 Nitrosomonas Nitrosomonadaceae 0.0 0.0 Sphingobacterium mizutaii Sphingobacterium Sphingobacteriaceae 0.0 0.0 Vibrio cholerae Vibrio Vibrionaceae 0.0 0.0 Rickettsia felis Rickettsia Rickettsiaceae 0.0 0.0 * * Moraxellaceae 0.0 0.0 Bacteroidetes bacterium 0.0 0.0 WD317 Corynebacterium casei Corynebacterium Corynebacteriaceae 0.0 0.0 Corynebacterium marinum Corynebacterium Corynebacteriaceae 0.0 0.0 * * Flavobacteriaceae 0.0 0.0 Caulobacter segnis Caulobacter Caulobacteraceae 0.0 0.0 Lactobacillus gasseri Lactobacillus Lactobacillaceae 0.0 0.0 * Meiothermus Thermaceae 0.0 0.0 Rhizobium sp. NT-26 Rhizobium Rhizobiaceae 0.0 0.0 Bacillus coagulans Bacillus Bacillaceae 0.0 0.0 * Sphingomonas Sphingomonadaceae 0.0 0.0 * Brevibacterium Brevibacteriaceae 0.0 0.0 Nitrosomonas europaea Nitrosomonas Nitrosomonadaceae 0.0 0.0 Pseudomonas alkylphenolia Pseudomonas Pseudomonadaceae 0.0 0.0 Terrabacter sp. DBF63 Terrabacter Intrasporangiaceae 0.0 0.0 beta proteobacterium CB 0.0 0.0 Moraxella ovis Moraxella Moraxellaceae 0.0 0.0 Shewanella baltica Shewanella Shewanellaceae 0.0 0.0 Mycobacterium gilvum Mycobacterium Mycobacteriaceae 0.0 0.0 * Exiguobacterium 0.0 0.0 * Ochrobactrum Brucellaceae 0.0 0.0 Geodermatophilus obscurus Geodermatophilus Geodermatophilaceae 0.0 0.0 * Devosia Hyphomicrobiaceae 0.0 0.0 Moraxella osloensis Moraxella Moraxellaceae 0.0 0.0 Exiguobacterium sp. 11-28 Exiguobacterium 0.0 0.0 Nocardioides sp. JS614 Nocardioides Nocardioidaceae 0.0 0.0 Nocardioides sp. USM2 Nocardioides Nocardioidaceae 0.0 0.0 Burkholderia gladioli Burkholderia Burkholderiaceae 0.0 0.0 Renibacterium Renibacterium Micrococcaceae 0.0 0.0 salmoninarum Pseudomonas syringae Pseudomonas Pseudomonadaceae 0.0 0.0 group genomosp. 3 Bifidobacterium Bifidobacterium Bifidobacteriaceae 0.0 0.0 pseudolongum toluene-degrading bacterium 0.0 0.0 UCR 021t Corynebacterium imitans Corynebacterium Corynebacteriaceae 0.0 0.0 Corynebacterium callunae Corynebacterium Corynebacteriaceae 0.0 0.0 Bosea sp. WAO Bosea Bradyrhizobiaceae 0.0 0.0 Xanthobacter autotrophicus Xanthobacter Xanthobacteraceae 0.0 0.0 Corynebacterium diphtheriae Corynebacterium Corynebacteriaceae 0.0 0.0 Bacteroides thetaiotaomicron Bacteroides Bacteroidaceae 0.0 0.0 Caulobacter vibrioides Caulobacter Caulobacteraceae 0.0 0.0 * * Pasteurellaceae 0.0 0.0 Finegoldia magna Finegoldia Peptoniphilaceae 0.0 0.0 Anaerococcus prevotii Anaerococcus Peptoniphilaceae 0.0 0.0 Azorhizobium caulinodans Azorhizobium Xanthobacteraceae 0.0 0.0 Erysipelothrix rhusiopathiae Erysipelothrix Erysipelotrichaceae 0.0 0.0 Porphyromonas Porphyromonas Porphyromonadaceae 0.0 0.0 asaccharolytica * Sphingopyxis Sphingomonadaceae 0.0 0.0 Eubacterium rectale Eubacterium Eubacteriaceae 0.0 0.0 Acinetobacter venetianus Acinetobacter Moraxellaceae 0.0 0.0 Variovorax paradoxus Variovorax Comamonadaceae 0.0 0.0 Acinetobacter sp. ED45-25 Acinetobacter Moraxellaceae 0.0 0.0 Bradyrhizobium Bradyrhizobium Bradyrhizobiaceae 0.0 0.0 diazoefficiens Bradyrhizobium japonicum Bradyrhizobium Bradyrhizobiaceae 0.0 0.0 Megamonas hypermegale Megamonas Veillonellaceae 0.0 0.0 * Methylobacillus Methylophilaceae 0.0 0.0 Nocardiopsis alba Nocardiopsis Nocardiopsaceae 0.0 0.0 Modestobacter marinus Modestobacter Geodermatophilaceae 0.0 0.0 Corynebacterium Corynebacterium Corynebacteriaceae 0.0 0.0 doosanense Blastococcus saxobsidens Blastococcus Geodermatophilaceae 0.0 0.0 Anoxybacillus flavithermus Anoxybacillus Bacillaceae 0.0 0.0 Aeromonas caviae Aeromonas Aeromonadaceae 0.0 0.0 Bacillus subtilis Bacillus Bacillaceae 0.0 0.0 Elizabethkingia anophelis Elizabethkingia Flavobacteriaceae 0.0 0.0 Staphylococcus hominis Staphylococcus Staphylococcaceae 0.0 0.0 Ruminococcus bicirculans Ruminococcus Ruminococcaceae 0.0 0.0 Paracoccus marcusii Paracoccus Rhodobacteraceae 0.0 0.0 * Psychrobacter Moraxellaceae 0.0 0.0 Sphingomonas Sphingomonas Sphingomonadaceae 0.0 0.0 sanxanigenens * * Leuconostocaceae 0.0 0.0 * * Burkholderiaceae 0.0 0.0 Bacteroides vulgatus Bacteroides Bacteroidaceae 0.0 0.0 Rhodopseudomonas palustris Rhodopseudomonas Bradyrhizobiaceae 0.0 0.0 Pantoea agglomerans Pantoea Enterobacteriaceae 0.0 0.0 * * Sphingomonadaceae 0.0 0.0 Mycobacterium kansasii Mycobacterium Mycobacteriaceae 0.0 0.0 * Streptomyces Streptomycetaceae 0.0 0.0 Enterococcus faecalis Enterococcus Enterococcaceae 0.0 0.0 Acinetobacter sp. NFM2 Acinetobacter Moraxellaceae 0.0 0.0 Shewanella putrefaciens Shewanella Shewanellaceae 0.0 0.0 Bifidobacterium adolescentis Bifidobacterium Bifidobacteriaceae 0.0 0.0 * Bifidobacterium Bifidobacteriaceae 0.0 0.0 Porphyromonas gingivalis Porphyromonas Porphyromonadaceae 0.0 0.0 Neisseria meningitidis Neisseria Neisseriaceae 0.0 0.0 Rhodococcus pyridinivorans Rhodococcus Nocardiaceae 0.0 0.0 Aeromonas salmonicida Aeromonas Aeromonadaceae 0.0 0.0 Planococcus sp. PAMC Planococcus Planococcaceae 0.0 0.0 21323 Pseudomonas simiae Pseudomonas Pseudomonadaceae 0.0 0.0 Faecalibacterium prausnitzii Faecalibacterium Ruminococcaceae 0.0 0.0 Acinetobacter lwoffii Acinetobacter Moraxellaceae 0.0 0.0 Exiguobacterium sp. N139 Exiguobacterium 0.0 0.0 Streptococcus anginosus Streptococcus Streptococcaceae 0.0 0.0 Thauera sp. MZ1T Thauera Rhodocyclaceae 0.0 0.0 * Shewanella Shewanellaceae 0.0 0.0 * Aeromonas Aeromonadaceae 0.0 0.0 Staphylococcus aureus Staphylococcus Staphylococcaceae 0.0 0.0 Aeromonas hydrophila Aeromonas Aeromonadaceae 0.0 0.0 Aeromonas veronii Aeromonas Aeromonadaceae 0.0 0.0 Bifidobacterium breve Bifidobacterium Bifidobacteriaceae 0.0 0.0 Bacillus megaterium Bacillus Bacillaceae 0.0 0.0

SUPPLEMENTARY TABLE 4 Normalized Reads Per million DNA Patient genus family DNA PC CSF Filobasidiella Tremellaceae 10,696.4^(#) 0.0 Aspergillus Aspergillaceae 1567.5^(#) 0.0 Toxoplasma Sarcocystidae 947.9^(#) 0.0 Filobasidiella Tremellaceae 53.3^(#) 0.0 * Sarcocystidae 27.3^(#) 0.0 Neospora Sarcocystidae 5.1 0.0 Hammondia Sarcocystidae 2.5 0.0 * Aspergillaceae 2.1 0.0 Hammondia Sarcocystidae 2.1 0.0 Aspergillus Aspergillaceae 2.0 0.0 Hammondia Sarcocystidae 1.6 0.0 Aspergillus Aspergillaceae 1.5 0.0 Aspergillus Aspergillaceae 0.7 0.0 Bartheletia Bartheletiaceae 0.7 0.0 Filobasidiella Tremellaceae 0.7 0.0 Aspergillus Aspergillaceae 0.6 0.0 Spirometra Diphyllobothriidae 0.5 0.2 Mucor Mucoraceae 0.5 0.0 Aspergillus Aspergillaceae 0.4 0.0 Setosphaeria Pleosporaceae 0.4 0.0 Caenorhabditis Rhabditidae 0.3 0.0 Aspergillus Aspergillaceae 0.3 0.0 Anisakis Anisakidae 0.2 0.0 Penicillium Aspergillaceae 0.1 0.0 Plectosphaerella Plectosphaerellaceae 0.1 0.0 * Tremellaceae 0.1 0.0 Malassezia Malasseziaceae 0.1 0.0 Gongylonema Gongylonematidae 0.1 0.0 Wallemia 0.1 0.0 * * 0.1 0.0 Candida Debaryomycetaceae 0.1 0.0 Elaeophora Onchocercidae 0.1 0.0 Lichtheimia Lichtheimiaceae 0.0 0.0 Brugia Onchocercidae 0.0 0.0 Parastrongyloides Strongyloididae 0.0 0.0 Sordaria Sordariaceae 0.0 0.1 Bipolaris Pleosporaceae 0.0 0.0 Malassezia Malasseziaceae 0.0 0.0 Cladosporium Cladosporiaceae 0.0 0.0 Strongylocentrotus Strongylocentrotidae 0.0 0.0 Caenorhabditis Rhabditidae 0.0 0.0 Penicillium Aspergillaceae 0.0 0.0 Fusarium Nectriaceae 0.0 0.0 Leptosphaeria Leptosphaeriaceae 0.0 0.0 Chaetomium Chaetomiaceae 0.0 0.0 Myceliophthora Chaetomiaceae 0.0 0.0 Alternaria Pleosporaceae 0.0 0.0 Saccharomyces Saccharomycetaceae 0.0 0.0 Aphanomyces Saprolegniaceae 0.0 0.0 Ophiostoma Ophiostomataceae 0.0 0.0 Rhodotorula 0.0 0.0 Trichosporon 0.0 0.0 Phoma Didymellaceae 0.0 0.0 Erysiphe Erysiphaceae 0.0 0.0 Leptosphaeria Leptosphaeriaceae 0.0 0.0 Zymoseptoria Mycosphaerellaceae 0.0 0.0 Debaryomyces Debaryomycetaceae 0.0 0.0 Alternaria Pleosporaceae 0.0 0.0 Phoma Didymellaceae 0.0 0.0 Umbilicaria Umbilicariaceae 0.0 0.0 Penicillium Aspergillaceae 0.0 0.0 Wuchereria Onchocercidae 0.0 0.0 Cladosporium Cladosporiaceae 0.0 0.0 Saccharomyces Saccharomycetaceae 0.0 0.0 Fusarium Nectriaceae 0.0 0.0 Dasybranchus Capitellidae 0.0 0.0 Neofusicoccum Botryosphaeriaceae 0.0 0.0 Cladosporium Cladosporiaceae 0.0 0.0 Penicillium Aspergillaceae 0.0 0.0 Ustilago Ustilaginaceae 0.0 0.0 Mucor Mucoraceae 0.0 0.0 Albugo Albuginaceae 0.0 0.0 

What is claimed is:
 1. A method for identifying one or more pathogens in a sample of biological material, the method comprising, at a computer system: receiving a plurality of sequence reads obtained from a sequencing of DNA molecules from the sample of biological material, the sample including DNA molecules from a plurality of organisms; using a first alignment technique to align the plurality of sequence reads to a plurality of classified reference genomes in a database, thereby obtaining alignment results that include, for each of at least a portion of the sequence reads, at least one matching reference genome to which the sequence read aligns, thereby obtaining matching reference genomes, and wherein the alignment results include classification information for each of the matching reference genomes, the classification information including a plurality of taxonomy identifiers; filtering the alignment results that are initial alignments results by: for each matching reference genome of a set of matching reference genomes: based on the initial alignment results, identifying an optimally-aligning sequence read that aligns to the matching reference genome with an optimal alignment score that exceeds or is equal to alignments scores of other sequence reads that align to the matching reference genome, thereby obtaining optimally-aligning sequence reads; for each of the optimally-aligning sequence reads: applying a second alignment technique for the optimally-aligning sequence read to the classified reference genomes to obtain a plurality of new alignment scores; determining whether any of the new alignment scores exceeds the optimal alignment score for the optimally-aligning sequence read; when one new alignment score exceeds the optimal alignment score for the optimally-aligning sequence read, removing the matching reference genome from the set of matching reference genomes; identifying a set of the sequence reads that aligned to two or more of the classified reference genomes with at a least a minimum alignment threshold; for each sequence read of the set of the sequence reads: identifying two or more matching reference genomes of the classified reference genomes to which the sequence read aligns with at least the minimum alignment threshold; assigning a taxonomy identifier from the classification information to each of the two or more of the classified reference genomes, thereby obtaining assigned taxonomy identifiers, each taxonomy identifier including at least two levels of classification, the at least two levels having a hierarchy such that there is a lower level and at least one higher level; comparing the assigned taxonomy identifiers of the two or more classified reference genomes at each of the at least two levels of classification; removing each level of the at least two levels from the assigned taxonomy identifiers that do not match between the two or more classified reference genomes; and assigning to the sequence read a lowest level of the at least two levels of the assigned taxonomy identifiers that is shared between the two or more classified reference genomes; and providing an identification of one or more taxonomy identifiers corresponding to one or more candidate pathogens based on numbers of corresponding sequence reads assigned to each of the plurality of taxonomy identifiers.
 2. The method of claim 1, further comprising: when none of the at least two levels of the taxonomy identifiers match, associating an unassigned state with the sequence read.
 3. The method of claim 1, wherein providing the identification of the one or more taxonomy identifiers corresponding to the one or more candidate pathogens includes providing an amount of corresponding sequence reads.
 4. The method of claim 1, further comprising: using sequence reads having an assigned taxonomy identifier to identify the one or more taxonomy identifiers corresponding to the one or more candidate pathogens by: for each of a set of taxonomy identifiers: determining a total score based on at least a coverage of a reference genome corresponding to the taxonomy identifier, thereby obtaining total scores for the set of taxonomy identifiers; ranking the total scores; and identifying the one or more taxonomy identifiers exceeding a score threshold.
 5. The method of claim 4, wherein the total score is determined further based on: a length of the reference genome corresponding to the taxonomy identifier, and an identity score corresponding to a fractional sum over genomic positions in the reference genome, the fractional sum of a fraction of sequence reads having an exact match at a genomic position relative to all sequence reads covering the genomic position.
 6. The method of claim 5, wherein the total score for a taxonomy identifier is determined as (the length+the identity score)*a percent identity score, wherein the percent identity score corresponds to the identity score divided by a number of genomic positions covered by a sequence read aligned to the reference genome.
 7. The method of claim 1, wherein the sample of biological material is from a host organism, the method further comprising: performing a clinical intervention for the host organism based on the identification of the one or more taxonomy identifiers corresponding to the one or more candidate pathogens.
 8. The method of claim 1, further comprising: preprocessing the plurality of sequence reads to perform at least one of: remove sequence reads identified as low quality sequence reads, and remove sequence reads identified as low complexity sequence reads.
 9. The method of claim 1, wherein the classified reference genomes in the database include reference genomes of microorganisms, the method further comprising: prior to aligning the plurality of sequence reads to the classified reference genomes in the database: aligning the plurality of sequence reads to a reference human genome; and removing a portion of the plurality of sequence reads that align to the reference human genome.
 10. The method of claim 9, wherein the classified reference genomes in the database include bacterial genomes, fungal genomes, and viral genomes.
 11. The method of claim 1, further comprising: sequencing DNA molecules from the sample of biological material to obtain at least one million sequence reads, wherein the alignment technique aligns the at least one million sequence reads.
 12. A method for identifying pathogens in a sample of biological material, the method comprising, at a computer system: receiving a plurality of sequence reads obtained from a sequencing of DNA molecules from the sample of biological material, the sample including DNA molecules from a plurality of organisms; using a first alignment technique to align the plurality of sequence reads to a plurality of reference genomes in a database, thereby obtaining initial alignment results that include, for each of at least a portion of the sequence reads, a matching reference genome to which the sequence read aligns; for each of a set of matching reference genomes: based on the initial alignment results, identifying an optimally-aligning sequence read that aligns to the matching reference genome with an optimal alignment score that exceeds or is equal to alignments scores of other sequence reads that align to the matching reference genome, thereby obtaining optimally-aligning sequence reads; for each of the optimally-aligning sequence reads: applying a second alignment technique for the optimally-aligning sequence read to the plurality of classified reference genomes to obtain a plurality of new alignment scores; determining whether any of the new alignment scores exceeds the optimal alignment score for the optimally-aligning sequence read; based on one new alignment score exceeding the optimal alignment score for the optimally-aligning sequence read, removing the matching reference genome from the first alignment from the set of matching reference genomes; and providing the set of matching reference genomes.
 13. The method of claim 12, wherein the matching reference genome from the first alignment is removed from the set of matching reference genomes when the matching reference genome shares a same taxonomic level as a reference genome corresponding to the one new alignment score.
 14. The method of claim 12, further comprising: providing the initial alignment results corresponding to the set of matching reference genomes.
 15. The method of claim 12, wherein the initial alignment technique uses a global alignment algorithm and the second alignment technique uses a local alignment algorithm. 